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FOREWORD 


This report entitled "Development of Response Models for the Earth Radiation 
Budget Experiment (ERBE) Sensors" consists of the following four parts. 

This is Part X, NASA CR— 178292, entitled "Dynamic Models and Computer 
Simulations for the ERBE Nonscanner, Scanner and Solar Monitor Sensors". 

Part II, NASA CR-178293, is entitled "Analysis of the ERBE Integrating 
Sphere Ground Calibration". 

HI* NASA CR — 178294, is entitled "ERBE Scanner Measurement Accuracy 
Analysis Due to Reduced Housekeeping Data". 

NASA CR-178295 , is entitled "Preliminary Nonscanner Models and 
Count Conversion Algorithms". 
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1. INTRODUCTION 


The radiometric measurement system for the Earth Radiation Budget Experi- 
ment (ERBE) is designed to measure the earth's radiation fluxes which result 
from the earth and its atmosphere's emmitted energy as well as from reflected 
energy, i.e., the earth's albedo. A general overview of the concept behind the 
three-satellite ERBE radiometric system and its instruments has been described 
in many internal NASA documents, such as the Science Team Minutes and Contractor 
Status Reports, as well as in journal publications, (e.g.. Ref. 1). 

The radiometer package on each ERBE satellite consists of eight instruments 
distinguished by their mode of operation, f ield-of-view and spectral response. 
Four of the eight instruments are nonscanner (NS) type radiometers, two with 
a Medium Field-of-View (MFOV) and two a Wide Field-of-View (WFOV) . Two of the 
nonscanner radiometers, one MFOV and one WFOV, have a suprasil filter which 
eliminates the instrument's response to longwave radiative effects so that they 
only remain sensitive to shortwave radiative effects. Shortwave (SW) radiative 
effects are defined here as having wavelengths between 0.2 and 5 pm and long- 
wave radiative effects as having wavelengths between 5 and 50 pm. Since radio- 
meters without filters measure total radiation, they are designated "Total" (T) . 
The four nonscanner radiometers are therefore the Medium Field-of-View Total 
(MFOVT) , Medium Field-of-View Shortwave (MFOVSW) , Wide Field-of-View Total 
(WFOVT) and Wide Field-of-View Shortwave (WFOVSW) . Three of the remaining in- 
struments or "channels" are Narrow Field-of-View (NFOV) scanning radiometers, 
two of which have filters. One scanner radiometer has a filter which eliminates 
longwave radiative effects and another has a filter which eliminates shortwave 
radiative effects. Hence the three scanner radiometers are the Narrow Field-of- 
View Shortwave (NFOVSW) , the Narrow Field-of-View Longwave (NFOVLW) and the 



Narrow Field-of-View Total (NFOVT) . The eighth instrument in the ERBE radio- 
metric package is the solar monitor radiometer. 

After the ERBE radiometers perform their earth, space and/or sun view and 
calibration measurements, the measurement results undergo an analog- to-digital 
(A/D) conversion and the data now in counts is transmitted to the earth through 
the spacecraft telemetry system. Counts for the nonscanners are measured 
every 0.8 seconds. Heat sink temperatures are read every 8 seconds and other 
housekeeping (HK) temperatures every 16 seconds. Therefore, in 16 seconds 
the data record for one nonscanner consists of 20 radiometric count readings, 

2 heat sink temperature readings and 1 reading for each of the other HK 
temperatures. 

The scanners perform a measurement cycle once cvciy h seconds. During a 
normal, earth-viewing scan each scanner performs 8 space-look measurements, 62 
earth-view measurements and 4 calibration measurements. The Longwave and Total 
NFOV scanners utilize the Internal Blackbody (IBB) for calibration measurements, 
whereas the Shortwave NFOV scanner utilizes the Shortwave Internal Calibration 
Source (SWICS). During a Mirror Attenuator Mosaic (MAM) scan cycle, the NFOVT 
and NFOVSW scanners measure solar flux performing 8 space-look measurements, 

46 MAM measurements and 4 calibration measurements. Housekeeping (HK) data is 
sampled once per scan. 

The solar monitor uses the same type of Active Cavity Radiometer (ACR) that 
is used by the nonscanners. The three major differences between these ACR's 
are as follows: 

a) the solar monitor has a shutter to cut off the solar flux every 
32 seconds, 

b) the configuration of the solar monitor housing is totally different 
and 
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c) the solar monitor has no heat sink heater. 

Solar Monitor count measurements are made every 0.033 seconds and HK measurements 
are made every 4 seconds. 

The radiometric instruments described above will have inherently different 
behavior and responses to arbitrary input signals since no instrument can be 
exactly duplicated or subjected to exactly the same calibration or operating 
environment. Even when viewing the problems of sensor performance from this 
macroscopic point of view, there are many things that have to he dealt with and 
accounted for. For example, environmental contamination and aging effects in- 
fluence the material degradation of the sensors. Sensor behavior may also be 
questioned when the thermal load and leakage exceed the design level. One ex- 
ample of abnormal thermal load is in the sun-blip case, where the sensor is 
subjected to a sunrise and sunset during its orbit. 

In attempting to account for such problems found in the sensor's radiometric 
performance and capability, it is therefore necessary to develop a radiometer 
sensor model. A sensor model should describe both a transient and steady-state 
sensor response for a given radiative input. It must provide sufficiently simple 
analytical expressions for parameters of particular interest and describe the 
model's behavior, e.g., the model's sensitivity to changes in the modeling assump- 
tions. The model must also be able to incorporate changes resulting from in- 
strument design requirements during its construction as well as instrument aging, 
i.e., sensor material degradation and unexpected thermal load and leakage during 
its long-term orbital operation. The products obtained from the sensor model and 
simulation must be in a form that is useful to the users. The sensors' response 
must therefore be validated by updating model parameters and converted into data 
with useful engineering units. 
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This report entitled "Development of Response Models for the Earth Radiation 
Budget Experiment (ERBE) Sensors" consists of the following four parts. 

Part I, NASA CR-178292, is entitled "Dynamic Models and Computer Simulations 
for the ERBE Nonscanner, Scanner and Solar Monitor Sensors". 

Part II, NASA CR-178293, is entitled "Analysis of the ERBE Integrating 

Sphere Ground Calibration". 

Part III, NASA CR-178294, is entitled "ERBE Scanner Measurement Accuracy 

Analysis Due to Reduced Housekeeping Data". 

Part IV, NASA CR-178295, is entitled "Preliminary Nonscanner Models and 

Count Conversion Algorithms". 
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2 . SENSOR MODELING 


The following ERBE radiometer sensor models were developed by Information & 
Control Systems, Incorporated (ICS) to consider thermal and radiative exchange 
effects. These models consider the surface specularity and spectral dependence 
of a filter as well as the thermal conductive and radiative interactions among 
the nodes in an enclosure of a radiometer. An enclosure is a domain in which 
2 d instrument's geometry is divided into appropriate elements. These enclosures 
in an instrument's model are set up using real or fictitious boundaries to ef- 
fectively alleviate the computational load for defining configuration factors. 

The instrument's model also includes the electronic thermal control sytems which 
are coupled with the energy equations. A general form of the energy equations 
can be used for all the sensors, despite differences in the mode of operation, 
field-of-view and spectral response, as long as the radiative enclosures and 
thermal nodes in the sytem's sensor configuration comply with the first law of 
thermodynamics . 

A sensor model can be tested and evaluated by comparing its behavior with 
the actual results obtained from the sensor's ground and inflight calibrations. 

The sensor model includes additional uncertainty since it incorporates outputs 
from the calibration source models. Once the sensor and calibration source models 
produce sufficiently accurate results, they can be used to define the parameters 
for the count conversion procedure. 

The characteristics of a sensor's simulation model are such that it has 
a dynamic response to the time dependence of an arbitrary input signal. The 
model also has a spectral response which is due to the spectral dependence of 
a sensor's filter and its exposed surfaces. Thermal effects for the model can 
be determined from the conductive and radiative properties of the sensor. The 
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model's measurement error can therefore be determined after its time response, 
spectral response and thermal effects have been obtained. The model can simu- 
late the sensor’s response to the temporal incident radiation received through 
the field-of-view limiter in both the transient and steady-states. 

The purpose of modeling and simulating is therefore 1) to understand sensor 
performance, 2) to improve measurement accuracy by updating the model parameters 
and 3) to provide appropriate constants necessary for the count conversion 
algorithms. 

2.1 General Description of the ERBE Sensor Models 

2,1.1 Models 

A sensor's geometry can be appropriately divided by considering its ther- 
mal and radiative interactions as a consequence of its design. Figures 1 
through 6 provide diagrams of how sensor elements have been geometrically 
partitioned for the nonscanner, scanner and solar monitor. Each element has 
a time resonse to the energy flowing to and from its neighboring elements. 

One might assume that maximizing the number of sensor nodes would maximize 
the accuracy of the model’s prediction of energy from a target source. However, 
increasing the number of nodes would introduce additional sources of error 
since the large set of equations required would incorporate additional para- 
meter approximations. Such equation approximations would need to be made to 
calculate constants such as thermal resistance, exchange factors and surface 
spectral and specular characteristics. There is an accuracy tradeoff between 
increasing the number of nodes and approximating the sensor's size, material 
properties, geometry, enclosure surface characteristics, temporal incident 
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energy flux intensity and electronic control system. If the temperature 
gradients within a node and between nodes are small enough, a model approx- 
imation with only a few nodes would not jeopardize the accuracy required to 
perform the parameter estimates and error analysis. Factors which would de- 
crease these temperature gradients within a sensor include uniform thermal 
properties, simple and symmetrical geometry, uniformity in directional and 
spectral properties of an enclosure surface and small temporal change in in- 
cident energy flux. 

A nonscanner sensor has 9 nodes if it has a suprasil filter dome (Figures 
2 and 4) and 8 nodes if it does not (Figures 1 and 3). The sensor's geometry 
is divided into nodes based on material homogeneity, overall sensor size, thermal 
conductivity of materials, thermal resistance between materials and surface 
radiative properties. The 9 nodes of the MFOVSW model (Figure 2) are there- 
fore, as numbered: 1) the active cavity with a heating and sensor coil 2) the 
copper heat sink with a heater, 3) the reference cavity with a sensor coil, 4) 
the aluminum heat sink, which is partly covered by a heater, 5) the base .plate 
or substrate, 6) the f ield-of-view limiter, 7) the aluminum heat sink which 
houses the reference cabity, which is partly covered by a heater, 8) the pre- 
cision aperture and 9) the filter dome for shortwave channels. 

Similarly, the scanner model (Figure 5) is divided into 12 nodes for the 
total channel, 13 nodes for the longwave channel and 14 nodes for the short- 
wave channel. The solar monitor is divided into 19 nodes (Figure 6). 

In addition to their nodes, these sensor models incorporate fictitious sur- 
faces to define and simplify enclosure geometries. Since the specular prop- 
erty of each surface in an enclosure is taken into consideration, the simple 
geometry facilitates the exchange factor calculation, alleviating the compu- 
tational work load for determining the exchange factors. A net radiation 
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analysis is performed for each enclosure which has surfaces with specular and 
spectral characteristics. This analysis appropriately deals with any partial 
specular and diffuse characteristics of an enclosure surface (Ref. 2). The 
model can also easily include transient response and thermal conduction. This 
net radiation method utilizes the waveband approximation to an enclosure's 
surface which has a specific spectral dependence. 


2.1.2 Thermal Conductive and Radiative Int eractions 

A volumetrically defined node, or control volume, complies with the first 
law of thermodynamics as in the following equation 


de 

dt 


A q" + A 4> 
c c r 


( 1 ) 


The first and second terms of the right hand side of Equation 1 represent the 
conductive transfer through contacting surfaces and radiative transfer from 
surrounding elements of an enclosure, respectively. The rate of energy change 
in a single node, therefore, depends on the difference between the ingoing and 
outgoing energy fluxes through a control surface or "boundary". This rate of 
energy gain or loss is composed of the rate of energy stored in a node as well 
as the rate of energy source or sink as indicated in the following equation 


£ ■ ^ £ - * 

nere V is the volume of a node, y the density of a node, 
the temperature and Q the rate of energy, or heat, at a 
Combining Equations 1 and 2 yields 

v yC - A c ^ + A r ♦ + ^ 


( 2 ) 


C the specific heat, 
source or sink. 


(3) 
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Equation 3 represents an energy balance in a single node. In the geometrically 
partitioned model, the conduction term can be written in the following form 



where A is the contacting area of i element to j element, qj; the heat flux 
i-j 

by conduction through the contacting area between i and j elements. is the 

thermal resistance between 1 and J elements. The thermal resistance between nodes 
is defined by the following realationship 


'i-j 




k. A 


i i-j 


'i-j "i-j 




k j A j-i 


(5) 


where is the distance from the thermal center of a node to the center of the 

contacting surface, k the thermal conductivity and h the thermal conductance of 
the contact between i and j elements. The thermal center is normally coincident 
with the geometrical center if a node material has homogeneous and isotropic ther 
mal properties. The contact resistance (1/hA) is mainly dependent on the surface 
roughness and the contacting pressure. The well-machined surface has a regular 
and uniform roughness. In such a case, the contact resistance can be determined 
by the following empirical relations (Ref. 3) as functions of the surface rough- 
ness and the thermal conductivities of two contacting materials. 
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and thermal 


Let surfaces contact each other with surface roughness 6^ and ^ 
conductivities k and k 2 , respectively. Assume the contacting area per periodic 
2 1 


and the 


contact is Tra. Then the number of contacts per unit area n = ^ 

, » • A 

total number of contacts for two surfaces with total contacting area A is ^ 


Conductance can therefore be calculated: 


h 




where 


( 6 ) 



6 = 6 1 + 5 2 2<S 1 ° r 26 2^ 

n = 2.136 /r T , where n is a parameter for roughness 

_ _ c where c is the half-width of the surface cavity 

“ a + c ’ 

n ts — i-- = number of contacts per unit area 

k l ,k 2 = thermal conductivity of nodes 1 and 2 

<5^ >( $ 2 = root mean square surface roughness of nodes 1 and 2 


A - total node interface area 


a = radius of a contacting point 

For the radiative interchanges in an enclosure including both specular 
and diffuse surfaces, a radiative property such as a reflectance, p, can be 
divided into a diffuse and a specular component as well as shortwave and long 

wave bands. 

That is, for shortwave bands: 


i 
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and for longwave bands: 

d s 

P * = P A + P «, 

where the superscripts, d and s, signify the diffuse and specular, and the sub- 
scripts the short and longwave bands. For a semi-transparent or an absorbing 
medium, the radiative properties have the following relationship among the 
reflectance p, emittance or absorptance c, and transmittance, T. 

e +o +T =1 for shortwave 
s s s 

e. + p 0 + T. = 1 for longwave 

The energy fluxes incident upon and leaving a typical semi-transparent 

th 

surface of the enclosure are described in Figure 8. Consider the i inside 
surface area, A^, of the enclosure which neighbors with other enclosures. The 
quantities and M* are the rates of incoming and outgoing radiant energy 

while considering diffuse and specular reflectance of the surface per unit in- 
side area, respectively. The quantity M* does not include the rate of trans- 
mitted incident radiant energy from the other neighboring enclosures. The 
quantity E ± is the rate of incoming radiant energy to the non-opaque boundary 
surface from the neighboring enclosure. The quanitity <J>^ is the energy flux 
supplied by some external means to the surface to make up for the net radiative 
loss and thereby maintain the specified surface temperature. A heat balance 
at the surface provides the relation 

4>i ■ E i " 

A second equation results from the fact that the energy flux leaving the 
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surface having diffuse and specular characteristics is composed of directly 
emitted plus reflected energy. This gives 

M i = £ i 0 T i + P i E i + P i E i (8) 

Equation 8 is called an opaque radiosity. The total radiosity, considering 
the rate of transmitted radiant energy which is originally coming from the 
neighboring enclosure, is composed of the opaque radiosity M* plus the energy 
flux transmitted through the surface with transmissivity, T This gives 


M i = M* + T ± E ± 


(9) 


The diffuse total radiosity considering only a diffuse factor in Eq. 9 
and excluding the specular reflection term is described as the following 


m: 


£ i 0 T i + P i E i + 


i i 


( 10 ) 


The contribution to from the other diffusely reflecting surfaces in 
an enclosure can be represented in accordance with the relationship defined by 
the configuration factors in an enclosure with diffuse surfaces, except that 
now the configuration factors are replaced by the exchange factors, or 


where F is the exchange factor defined between i element and j element through 

ij 

j's images on the element i. Equations 7 through 11 can be represensted by the 
short and longwaves, respectively. 

To obtain the heat balance (Eq. 7) as functions of surface temperatures in 
an enclosure and the irradiance, or incident radiation, from the target source at 
the f ield-of-view limiter of the radiometer models, first Eq. 9 is substituted 

into Eq. 7, yielding 
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♦i - (1 - pj - Pi> E i - e i “ T i E i 

in which p d + p® - P ± where p i is the t0tal r P flectance of the surface and 
ft = 0 T 4 . 

Thus 


p 4 ) - c t a i ■ t i ' 

Then, combining Eqs. 10 and 11 and rearranging, yields the following equation 


E i * f <{ ij ' hi pd i rl ? i J <e J "l + T J V’ 

For brevity, we introduce new variables from Eq. 13, that is 

*ij - «« - hi op ' 1 hi v 
c u ■ <6 ij ■ f « pd> 1 F « v 

Thus, Eq. 12 with the above constants can be written 


(13) 


(14) 

(15) 


Equation 16 has an important role to couple an enclosure with the neighboring 
enclosures through a semi-transparent or transparent boundary surface. 

The above relation (Eq. 16) can be used for the cases when an enclosure 
is partly or totally surrounded by the neighboring enclosures and also when 
part or all of the surfaces of the enclosure are non-opaque, i.e., semi-trans- 
parent or transparent. To include all the configurations of the enclosures and 
to couple them together through interactive boundary surfaces, the above rela- 
tion (Eq. 16) can be generalized In the matrix form, that is, 


E = A ft + C E 


(17) 
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where A and C denote Eqs. 14 and 15, respectively. The bars above the notation 

A and C are to distinguish Eq. 14 and 15 from the notation for the area and 

the conduction constant of governing equation to be discussed in the Section 2.1.4. 


2.1.3 Coupling Enclosures 

In the above Eq. 17, E can be directly related to E by defining the geo- 
metrical coupling factor. The relationship between E and E is apparently de- 
pendent on the geometrical formation of enclosures through the interactive 
(non-opaque) surfaces as in the following example. Consider the system shown 
in Fig. 9 which consists of nine enclosures and interacts at the boundaries with 
the incident radiation from the environment. The incoming radiation E can be 
written 



(18) 


\jhere the superscript m denotes the number of enclosures in the prescribed 
system. Each enclosure of the system has a certain number of interacting sur- 
faces. In such a case, the incident radiation upon the surfaces of the i 
enclosure can be described by 



and 1 £ i £ m 


( 19 ) 
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where the subscript n denotes the number of a surface in an enclosure of the 
system. As shown in Fig. 9, eJ and E* define the incident radiation on an in- 
teractive boundary surface between the first and second enclosures. The sub- 
scripts 4 and 2 here, respectively denote the surface numbers as seen or counted 
from each enclosure. These surfaces designated by 4 and 2 in the enclosure 1 and 
2, respectively, are supposed to be the interactive boundary surface co-owned by 

enclosures 1 and 2. On the contrary, incident radiation E coming from the op- 

~ 2 

posite, or neighboring, enclosure upon the boundary surface is defined by E 2 and 
lj, respectively. Accordingly, observing enclosures 1 and 2 with respect to the 
incident radiation on the interactive boundary surface, we find the following 
coincidences : 

1 ~2 , _ 2 _ «1 ( 20 ) 

E 4 = E 2 and E 2 - 1 4 

Equation 20 tells us that the incident radiation E*, upon the interactive 
boundary surface of the first enclosure is virtually regarded as the opposite of 
incident radiation E 2 , on the back surface of the same boundary, whether the 
second enclosure is considered to be formulated or vice versa. This coincidence 
on each interactive boundary surface is uniquely determined by the formation 

of the enclosure in a given system. 

To generalize the above relationship, Eq. 20, for the enclosure which 
has interactive boundaries and also system boundaries, let us consider the 
first enclosure shown in Figure 9. The first enclosure has not only two inter- 
active boundaries at 3rd and 4th surfaces, but also two system boundaries at 
1st and 2nd surfaces (Fig. 9). Extending the above relationship (Eq. 20) to 
include the system boundaries and describing it in a matrix form, we have 
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( 21 ) 
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Likewise, the other enclosures have the same expression as depicted for the 
first enclosure. Thus, introducing a new variable -D for the enclosure coupling 
factor, the relationship (Eq. 21) can be written 

( 22 ) 

£ = DE + FS 

where D is a matrix representing the enclosure coupling factors through the 
interactive boundaries, F is a matrix representing the system boundary conditions 
imposed by surroundings, and S is a column vector of the heat fluxes arriving 
at the system boundaries. Now substituting Eq. 22 into Eq. 17 and collecting 
the same variables, we have, with I being an Identity matrix 

E = (I - CD)' 1 fil + (I - CD)” 1 CFS (23) 

Equation 23 illustrates that the incident radiation E upon a surface of the 
enclosure accounts for all the effects of radiation from the surfaces with 
different temperatures not only in its own enclosure, but also in neighboring 

enclosures, and the heat flux at the system boundaries. 

The energy balance at a node or a surface in the enclosure due to the 
incident radiation is then described by inserting Eqs. 22 and 23 into Eq. 12. 
Then, by eliminating E and E from Eq. 12, the radiative energy balance is 
obtained as only functions of the temperatures of enclosure surfaces and the 
irradiance from the target source or any other kinds of heat flux which inter- 
acts with the system boundaries. That is. 
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( 24 ) 


<j) - [(I - P - TD)(I - CD) -1 A - e ft] 

+ [(I - p - TD)(I - CD)" 1 C - T Fs] 

When the radiative energy balance described by Eq. 24 is applied to the 
enclosures obtained by properly dividing the geometry of the system, the 
number of the radiative energy balance equations, <}>, obtained from Eq. 24 is 
more than the number of nodal points because consideration is given twice to 
the interactive boundary between two enclosures. Accordingly, in the multiple 
enclosure system the radiative energy balance is calculated for each surface 
of enclosures. Then, for the overlapped boundary surface, the radiative energy 
balances from both enclosures are added together. The radiative energy balance 
(Eq. 24) is then plugged into the energy equation (Eq. 3) of the system. 

The enclosure coupling technique developed here has many potential 
applications and advantages over the conventional methods currently used. 

Among the advantages, this technique can be used for any kind of geometry 
by only setting up enclosures via dividing the geometry and defining the D 
and F matrices in Eq. 22. The enclosure space with irregular geometry can be 
properly divided into many smaller regular sub-spaces to adopt the enclosure 
coupling technique. For an enclosure with complicated geometry, the division 
of the space geometry into well-defined sub-spaces facilitates the computation 
of the configuration factor or the exchange factor, especially for an enclosure 
with specular surfaces which is totally dependent on the geometry. 

This technique can also be extended to treat radiation interactions in- 
cluding the scattering and absorption with optically thin or thick participating 
media in an enclosure. Even though there are temperature and density gradients 
in participating media, by dividing the enclosure space into sub-spaces which 
are small enough to assume the sub-spaces to be isothermal and isotropic, the 
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enclosure theory can be directly applied for the system without loss of gener- 


ality. 

In short, the advantages and potential applications of the enclosure coupling 
technique are listed below. 

ADVANTAGES : 

1. It minimizes the effort required to define the exchange factors 
among the surfaces in an enclosure with complex geometry, and so 
may be able to decrease a degree of possible approximation in the 
computation of exchange factors. 

2. It may take the wavelength characteristics of materials or participating 
media into account in a simple manner. 

3. It can be easily coupled with the thermal model of the system so that 
the material and thermal properties are considered in the coupled 
dynamic system without approximations and assumptions such as those 
in the boundary definitions. 

4. It enables the system to deal with any imposed initial conditions 
and prescribed boundary conditions. 

5. It may enhance the accuracy in the computation of a geometrically 
partitioned conductive and radiative model. 

6. As a consequence, this technique can avoid the length of time for 
computation and a cumbersome manipulation which is often revealed in 
the other methods such as ray tracing or Monte Carlo methods. 

APPLICATIONS : 

1. It can be used for any enclosure geometry. 

2. It can be used for either the enclosure with a vacuum space or the 
enclosure with participating media (solid, liquid or gas). 
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3. It can be applied to any optical system's anslvses and design regard- 
less of the directional, spectral, and/or temperature dependences of 
its enclosure spaces. 

4. Applications of the enclosure coupling technique to optical systems 
may vary from a typical radiometer to any waveguide system. 

2.1.4 Governing Equation 

The governing equation which describes the performance of a radiometer that 
has a dynamic response to an incoming radiance can be obtained by substituting 
Eqs. 4 and 24 into Eq. 3 and collecting terms. Then we can separate the in- 
coming radiation flux from the radiation balance term. The equation can be 
written in a general form to describe the energy balance in a node as follows 


AT + DT + BE + 


V 


+ b + e 


where 



e 


R 

M 


_1_ 

MR 

^ [(I - p - D)(I - CD )" 1 A - e] 

^ [(I - p - D) (I - CD )" 1 C - T f] 

, Ac. 

M 

: a bias term 
= a noise term 

= the thermal resistance between nodes (see Eq. 

vyc 


5) 


(25) 
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A = the area for the radiative interaction 
r 

A * the node area covered by the heat sink heater 
c 

Equation 25 represents a general form that consists of a set of equations 
for a multinode system. Considering only the thermal model for the nonscanner 
total radiometer which is divided into 8 nodes, we have a set of 8 equations 
which represent the energy balance in each node separately. The inclusion of the 
electronic thermal control system equations will increase the total system by 
the number of the transfer functions. The constants of Eq. 25 for the nonscanner 
the scanner, and the solar monitor are described in Appendices C, D, and E. 

2.1.5 Electronic Control System Model for the Nonscanner 

The block diagram shown in Fig. 10 illustrates the nonscanner ACR control 
system model. The control electronics transfer function is of the form (cf. 

ERBE Instruments, Monthly Status Report - Vol . 1, Nov. 1980, unpublished) 


K^CtjS + 1) (t 2 s + 1) 
F E^ = s(t 3 s + 1) (t^S + 1) 


(26) 


To obtain the differential equations describing the electronics transfer 
function, we define the state variables T 1Q , T^, T 12> as follows 


r W 1 

T 10 = V B K a h [W + R r (T 3 ) ' ®eJ + K E 


(27) 


Note B - — is an input constant, where and R^ are the bridge 

E R 3 + R/j 

resistances. 


w ■ A i + B i T j + c i T r 

hi ' t ( ' T n + ho + VlO> 


i = a, r; j = 1,3 


(28) 

(29) 
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FIGURE 10. Active Cavity Radiometer (ACR) Control System Model 



Substituting Eq. 27 into Eq. 29 gives 

t <- T n +T io> +Tl 
3 


*11 - T„ <- T U + V + ^ V B K a K E _ R a (l') l R r (T 3 )_ + t 3 K E n E (30) 


T = — (-T,„ + T 31 + T_ T .) 
12 T v 12 11 2 11 


(31) 


For full expression, substitute T u in Eq. 30 intoEq. 31. 

Note that T 12 is now the output of the electronics, so that the voltage 

output is 

(32) 

v " T 12 + n v * 

To convert to counts, C, use the following relation: 


C - V - K c T 12 + K c n v . 


K = 819.1 counts/volt, 
c 


The heat flow provided to the active cavity to maintain its temperature 
at a constant level is given by 


Q x " V 2 " K a T 12 + 2K a T 12 n v + K a n v 


(33) 


This expression for ^ brings to light the major effects that the noise 
terms introduce into the operation of the radiometer. Note that due to the 
use of square feedback, the linear noise terms acquire a gain of 2K & T^. 
while a quadratic noise term (having non-zero mean) is introduced. Coils 
which are wrapped around the active and reference cones generate heat and are 

expressed by the following terms* 

2 M*l> 

Qj i " B [ R a (T l ) + WJ 

2 'riV 


3 J 3 ' V » [R a (T i> + R r ( T 3 )J' ! 
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To obtain the closed-loop system, the expression Eq. 33 is simply 
substituted into the equation governing the active cavity temperature. 


. k „ 2K T.. n 

*1 ■ A 12 (T 2 - V + l D ik T k + B 1 E + ST T 11 + V V 


K n 2 
a v 


V T 1 > 


+ «i +Vb 

The block diagram shown in Fig. 11 illustrates the non-scanner ACR heat 
sink heater control system model. 

The control electronics transfer function can be approximated in the 

form 


( 34 ) 


F HE (S) “ ’SiE 


<v + x) 


(35) 


The, differential equations describing the heat sink heater controller can be 
expressed as 


*13 " \ 3 *SiE (I H ” T 2 ) + ‘SlE n H 

V H = T 13 + T H SiE \ 3 (X H " T 2 ) + T H ‘SlE °H 

«2 ■ V H 


(36) 

(37) 

(38) 


In Eq. 25 the heat source term, Q, is substituted by Eq. 38 for the aluminum 
and copper heat sinks. 


2.1.6 Simulation of the Nonscanner Model 

The simulation of the ERBE nonscanner sensors is divided into two major 
parts; the first determines the coefficients of the system, the second solves 
the systems over time and plots the results. Figure 12 is a flow diagram 
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FIGURE 12. Flow Diagram Indicating File and Program Usage in Simulation 







indicating file and program usage. The equation describing the change of temp- 
erature in a given node at time i is 

T(i) - AT(i) + DT 4 (i) + BE + B Q Q + B n n (39) 

Where A is the conductance coefficient matrix, D is the radiative coefficient 
matrix, B is the irradiance source vector and B^ is a vector of coefficients 
which describe heating contributions from the control system. T(i) represents 
either the temperature change for a thermal node or a voltage change for an 
electronic node. The general method of solution is a fourth order Runge-Kutta 
algorithm, all of which is discussed in great detail in the following sections. 

2. 1.6. A Parameter and Coefficient Determination 

The coefficient matrices and vectors are determined in the first program 
in the sequence ERBEPRM. The data inputted to this is determined from material 
specifications and engineering drawings supplied by the manufacturer. The code 
is partitioned into 5 major sections, not including the program input/output. 

The first of these is the subroutine "COMPUTA". This routine takes the heat 
capacities and node to node thermal resistances and returns the conductance 
coefficient matrix. The program then enters the routine ADDELC which is designed 
to modify the conductance matrix in such a way as to allow the ACR heater con- 
trols to use the same equations in the simulation as the node temperatures. 

The output from this routine is an enlarged conductance matrix and several elec- 
tronic vectors. Once these two routines have been computed COMPUTH calculates 
the radiative exchange coefficient between surfaces for each cavity. Due to 
the linking of cavities a and b, as defined in Appendix C, dome filter or 
fictitious surfaces, the radiative exchange coefficients for these cavities 
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are stored for later use. The results of COMPUTH are then added together to 
give the total radiative exchange coefficient matrix in COMPUTD, COMPUTB and 
CMPUTBP. The results are then issued to two files, one formatted for in- 
spection and one binary for use by the simulation program. 


2.1.6.B Instrument Simulation 

The program consists of four logical sections, an Input, the Integration, 
the intermediate Output and the Plotting section. The actual solution to the 
system of equations that describes each instrument is accomplished by the pro- 
gram ERBESIM. The program steps through a determined time interval using an 
evenly spaced discrete time step and solves each system using a Fourth Order 
Runge-Kutta scheme to evaluate the system at the next step. Each system is 
sampled at another rate controlled by the user until 20 samplings are taken 
or the end of simulated time is reached. The results are then stored on mag- 
netic medium for use by the plotting routines. The final product is the sen- 
sor's resposne described by plots and tables of the modeled response of the 
instrument, i.e., time-dependent temperatures of the nodes, instrument response 
(in counts), electronic voltages (in volts) and input source irradiance (in 

W/m 2 ). 

The first section of the program reads the coefficient matrices and vectors 
intitial conditions from binary files, and the simulation control variables 
from a namelist file. The Namelist file contains variables for controlling 
simulation time step size, sampling rates and the combination of instruments 
and radiation profiles, etc. which are to be simulated. The simulation is 


35 



capable of allowing the user to choose either a sawtooth, Fourier series, step 

function, constant or ramp radiance profile. 

The Runge-Kutta method used was developed by Fehlberg and solves the system 

5 

Y = Y + h .E n A. K.., i = 1, 2, 3, ... number of nodes 

i,m+l i,m j=0 J Ji 

where 

i 

K jl ' t(t m + C j h> y m + h 1-0 V V 

where f is the function describing T(i) from Eq. 39 as functions of time and 
temperature, a., b , and C. are constants determined by Fehlberg for solving 
first order non-stiff systems of ordinary differential equations. 

2,1.7 Count Conversion for Nonscanner 
2. 1.7. A Nonscanner Total Channel 

In the active cavity radiometers, the electronics control the amount of 
heat added to the active cavity to maintain the temperature of the cavity 
about half a degree above the reference cavity. The control of the heat input 
to the active cavity is inversely proportional to the incident radiation en- 
tering through the field of view and precision aperture from a given source. 

In other words, the control system provides more heat input for the weak incident 
radiation and less heat input for the strong incident radiation. Thus, the 
voltage driving the active cavity heater provides the sensor output proportional 
to the incident radiation. The voltage reading, after A/D conversion, provides 
the sensor output in counts. The sensor output in counts is then converted into 
engineering units. The algorithms for the count conversion are of the gain/offset 
type, and are derived based on the sampled instrument data, e.g., sensor response 
(in counts), HK temperatures, etc. and the instruments radiative and conductive 
interactions. 
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The following is the count conversion equation for NS total channels. 


E(t - T) A y V (t) + A^T^t) - T H J + A a |t a (c) - T a J 


+ ^[V 0 - t f ] + “ 


The constants in Eq. 40 are defined as follows 


B K 

A = - - a - 
v B, 


a _ JL 2 n 

*H = ~ B. kM.F lk k 

1 c 

. °18 3 

A A ■ - 4 ‘bT T 8 
1 o 


A A 16 T 3 

a f " - 4 IT T 6 

1 O 


B _ hi b _ £ ^ t 4 - ^ 

B ■ B, b E k=l B. k 4R_ M. B. 

1 1 O D 1 1 


where 

V = instrument output (in volts) 

T u = measured heat sink temperature (°K) at time t 

T^ * measured precision aperture temperature (°K) at time t 

T„ = measured FOV limiter temperature (°K) at time t 
r 


( 40 ) 

( 41 ) 

( 42 ) 

( 43 ) 

( 44 ) 

( 45 ) 


The remaining quantities, i.e., B , K , B, , T, D , , A „ and b are con- 

sal ik 1/ L 

stants which are specified for each instrument. These constants for the count 
donversion equations are determined for the simulation study of the instrument 
model using ground and inflight calibration data. These parameters or constants 


37 



may be altered at various points in time, usually after a flight calibration, 
in order to reflect any changes in the instrument characteristics. 

2 . 1.7.B Nonscanner Filtered Channel for Shortwave 

For the filtered channels, the count conversion has a set of two equations 
for short and longwaves as follows: 



A 

E 

+ a, „ AT.. = 

6. 

(46) 

11 

s 

12 9 

1 


21 

e 

s 

+ a 22 AT 9 = 

*2 

(47) 


The coefficients of the terms in the above equations are given by the 
following equations. 


“u ■ B ls • B u 
“12 ' 4D 1,9 T 9„ 


a 21 B 92 ” B 9 s 


“22 " *9,4 + *9,8 ‘ 4D 9,9 T 9 


n. 




3 V B g r 

^1 * (A 12 b E - jJjj. D il ARg 4 ^k^8,6 ik k Q i6r 2 

* ( T H ' T H } “ 4D 18 T 8 (T A ' V ' 4D 16 T 6 o (T F “ V 


(48) 

(49) 

(50) 

(51) 


- B K v - B E 
a a s ie 


(52) 


61 - "Jl °9k T k o + *98 (I A ' T 9 o ) + *d8 (T H - T 9 o ) + 4< W4?8,9 D 9k T k o > 
* (T H ■ T H ) + 4D 98 T 8. (T A ' V ' 4D 16 T 6 o (T F " T F () ) + B U E 


(53) 
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The filtered channel count converstion is obtained from the solution of 
two simultaneous equations (Eqs. 46 and 47). Thus the solution to the equations 
is given by 


a ll a 22 " a 12 a 21 


( a 22 ^1 “ “l2 ^2^ 


(54) 


T = ~ 

9 a 


11 a 22 " a l2 a 21 


(- a 21 + a lx B 2 ) 


(55) 


The longwave incident radiation estimate is then defined by subtracting 
the shortwave incident radiation from the total incident radiation estimate. 
That is, 



The constants in Eqs. 46 - 56 are 

v ■ shortwave channel output (in volts) 

s 

2 

E » shortwave incident radiation estimate (W/m ) 
s 

/s 2 

= longwave incident radiation estimate (W/m ) 

2 

E = total incident radiation estimate (W/m ) 


The temperatures of the heat sink, aperture and FOV limiter, T R , T^, and T p , 
are measured as before in total NS channel. 
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2.1.8 Sola r Monitor M o del Analysis 

The solar monitor is different from the nonscanner in size, shape and 
design feature as shown in Figure 6. It has a barrel structure - 9.5 cm long 
between the primary and secondary apertures. The tail of the core body of the 
instrument is attached to the housing and includes the primary and secondary 
apertures, barrel, heat sink and active and reference cavities. This instru- 
ment has a shutter located in front of the secondary aperture. The shutter 
operates on a 64 second -on and off" cycle in which it is open for 32 seconds 

and then closed for 32 seconds. 

Since this instrument deals with a strong intensity of the incident 
radiant energy and does not have the heat sink temperature controller, the 
heat sink has a larger heat capacity for damping out the effect of large in- 


put temperature gradients. 

However, it is desirable to check whether there is any transient effect 
due to the periodical on and off mode of the shutter, the strong intensity of 
solar radiation, and the excessively long barrel geometry. In other words, 
it is desirable to check whether the steady-state approximation is appropriate 
to develop the count conversion algorithm for the solar monitor. Thus, we 
introduce two simplified geometries here to check how fast the temperature 
reaches a stable condition, as every 32 seconds the solar radiation is 

viewed. 

First, let us assume that the shutter has a rectangular shape of thin 
aluminum plate and the shutter is partly exposed to the incident solar radiation 
for a 32 second interval of the full 64 second cycle and has a radiation inter- 
action with the isothermal neighboring components at any time. Figure 13 shows 
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the shutter model. 



FIGURE 13. The Shutter Model 


It is assumed that the left end of the shutter is insulated and the right 
end is attached to a temperature reservoir which has a large heat capacity. 
Then the energy balance equation, the initial and boundary conditions govern- 
ing the above model can be written as 


G.E . 



k D 2 T ^ q" 
pC 3x2 5 


(57) 


I.C . T(T - 0) - T q 


B.C • T x (x = 0, T) * 0 
V* - L, T) - T b 


The above equation with Neumann and Dirichlet boundary conditions can be 
easily solved by introducing a new variable and separating the variables. 
Thus the solution to the equation is described by 


(-l) n -aXn 2 x 


T(x, T> - T + 2(T o -T b ) n J 0 e 


COS(Xnx) 


' 4 n=0 OS L)* 6 ’ aXn T C0S(Xnx) 
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where T is an initial temperature set by 273.16 °K, L the total length of the 

o 

shutter (6.5 cm), 6 the thickness of the shutter (0.2 cm), k the thermal con- 
ductivity (1.557 w/cm°K), i, the exposure length to solar radiation (2.5 cm), 

p the density of the shutter (2.713 g/cm 3 ) , c the specific heat of the 

2 

shutter (0.887 J/g°K). a the thermal wave velocity (0.667 cm /sec) 


. tfL+m „ . 0. 1. 2. ... 

n 2L 

2 

q" is the solar flux (0.142 w/cm ) x 0.2 
s 

4 4 

q" = Fee (T - T ] 

H cx> 00 

q" = q" u(£-x) u(32-t) + q'^ u(£-x) u(x-64) u(96-t) ... + 2 q 
F the configuration factor between the shutter and the neighboring 
components (0.85) 

e the emissivity of the shutter surface (0.8) 
a - the Stephan-Boltzmann constant 
T = the ambient temperature (use = T b ) 

OO 

T — the reservoir temperature set by 273 K, 
b 

X — current shutter temperature, but use T(t 1) 
u(x) = the unit step function. 


The graphical description of the solution is shown in Figures 14 and 15. 
Figure 14 shows that the shutter temperature not only responds to the on and off 
mode of the solar flux, but the peak temperatures as well as the minimum temp- 
eratures of each cycle (64 sec) also stay at the same ranges. That is, the time 
averaged temperature of the shutter does not change sifnif icantly . And Figure 
14 also shows that the temperatures before reaching minimum and maximum points 
become stable. Figure 15 shows the temperatures along the length of the shutter 

after 100 seconds passed. 

Secondly, the long barrel which Is placed In the housing and Is between 
the primer, end seconder, aperture, is considered by the following assumptions. 

1) no direct exposure to the solar flux 
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TErtPERATURE ( DEG K) 


273 



FIGURE 15. SHUTTER TEMPERATURE AT 100 SECONDS WITH 20% ABSORPTION 
OF SOLAR FLUX Qs = 0.142 W/CM . 
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2) radially uniform temperature since the wall thickness is thin 

3) no circumferential temperature gradient 

4) radiation exchanges only at inside and outside walls 
The configuration of the model is shown in Figure 16. 



FIGURE 16. A SIMPLIFIED CONFIGURATION OF SOLAR MONITOR BARREL 

The equations and boundary conditions governing the thermal conditions 
imposed upon the above geometry can be written as 

3T , 3 2 T q" 

G.E. pc ^ = k 3^2 + R o - Ri 

I.C. T (T=0) = T, 


(59) 


B.c. T (x=0) * (set to be constant) 


T(x=L) = T. 


Then the solution to the equation is obtained as 

2k(R o -R 1 ) q” 


T “ T b + 2k(R 


- 2 


q n . lL { h - (*) 2 1 - 

wU L J 

f to . 

n*0 l(LXn)^ 


• M T 

q L 


a - 


2k(R o -R.) q£ 1 

— 5^ axil yz 


b 


e ~aXn t Cog (X nx ) 


II 


(60) 
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where 


q" = qV + q^ 


" = F, .e.a[(T A (x, t-5) - T (x, t-1))] 


1-11 


qV = q" 

M i 1 s 


'2 1-2 


Q" = F e a[ T A (in x) - T 4 (X=0, t-l)J[u(64-T) u(x-32) + u(128-t) u(t-96) + 

^oo 1-3 3 ~ ,,r 


.] 


av 


Tb 
L = 
k - 
a = 
Ro- 
An 


= 273.16 °K 
10 cm 

1.557 w/cm °K 

0.647 cm 2 /sec 

Ri = 0.5 cm 
(2n+l)_ 


til 


1- 


£ - ir 

2L _ 

0.146 £u(32-t) + u ( t-64) u(96-t) + ...J 


= 0.75 


1- 


1 

= 0.8 
2 = 0-4 
= 0.08 


1- 


= 0.4 


3 

= 0.5 


The above values are used in both the shutter and barrel cases. These estimates 
were selected to best represent the solor monitor based on the material and 
thermal properties of the shutter and barrel. The results of the calculations 
described in Figures 17 and 18 demonstrate that the temperatures do not change 
very much with respect to time or in the axial direction of the barrel since 
they exhibit only about a half degree of temperature change in both cases. 
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THETft TEMP ( DEG K) 



FIGURE 17. BARREL TEMPERATURE WRT TIME AT X*3 CM FROM THE SECONDARY 
APERTURE. 


THETA TEMP C DEG K) 



FIGURE 18. BARREL TEMPERATURE ALONG THE AXIAL DIRECTION AT THE TIME OF 200 
SECONDS. 



Therefore, these results show that the steady-state approximation for the 
count conversion algorithm does not contain any crucial errors. Therefore, 
the open and closed modes of the shutter were considered in the algorithm 
development. Namely, the output responds to the direct radiation which falls 
into the active cavity according to the shutter's on and off modes. Even 
though the instrument has a small time constant and hence it is stabilized 
quickly, within the 32 seconds of a half cycle, the output data which is used 
for the conversion algorithm has to be picked up at the near-end of a half 
cycle just before the shutter closes the aperture. 
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2.1.9 Solor Monitor Count Conversion Algorithm 


The approximations made here including the steady-state approximation are 
the following: 

1. The nodal points are selected to be adjacent to where the temperature 
probes are so that the temperatures measured by the probes may be 
regarded as the temperature of the nodes. 

9 a. _ t + b where b is the electronical bias, 
l 2 E E 

3. The temperature of the reference cavity is the same as that of the 

heat sink, that is, T 2 = Accordingly, T 1q - T 3o + b £ where the 

subscript o signifies the initial temperature. 

4 . t^ - T lo can be replaced by T^ - T 3q . 

With the above approximations, the equation for count conversion was derived as 

E(t-T) - A v V 2 (t) + Aj [T 3 (t) - T Jo ] + Aj [T 5 <t) - tJ 

+ *7 [l 7 ( t) - T 7o] + “ (<>1) 

where T 3 , T 5> and T ? are temperatures in °K which are measured in the instrument. 
These temperatures are measured and telemetered every 16 seconds. If these 
temperatures are measured twice during the open period of the shutter, then the 
second measured temperatures can be used for the count conversion algorithm. 
Otherwise, as the temperatures are measured arbitrarily without regard to the 
shutter on-and-off mode, the temperature measured at or near the time before the 
shutter closure must be used. As well as the temperatures measured, the count 
has to be picked up at every 40th data point or near the time before the shutter 
closure. V(t) in the equation is the instrument output in volts and is tele- 
metered at every 0.8 second intervals. The constants in the equations are 
obtained from the thermal control electronics and conductive and radiative 
heat transfer, all of which interact with the active cavity and the nodes 
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which have temperature probes. In the equation, 


B.K, 

Bl 


a 3 - - ^ [D^iTio + D 1 _ 2 T^ 0 + D^T^ + DWl| 0 + 1 >i-6 T 6o 
+ D l-18 T 18o + D 1-19 T 19 o] 
a 5 “ -^[ d 1-5 t L + D 1-6 T 6 o] 

A 7 “ TJ- [di_ 7 t 5 0 + D^gT J 0 + Hi-9^96 + D l-ll T ?lo + D 1-1A T 14 o 
+ d 1-15 t 156 + D 1-17 T 17 o] 


Ayi . _ E 0^ 4 _ BjKjV 

b E ^k Bl ko Bj 


B l‘ 


( 62 ) 

(63) 

(64) 


(65) 

( 66 ) 


the quantities such as B a> B^, K & , Kj, V, B^, T, , A^ a °d bg are the 

o 

constants to be specified for the instrument. 


2.1.10 Scanner Control Electronics 

The control electronics for the scanner have a different feature from 
that of the nonscanner. The difference between them is attributed to the 
type of sensor in the scanner. The scanner uses thermister blolometers for 
both active and reference flakes. When the temperature of the thermister 
varies due to the energy exchange, the current through the thermister varies 
with the change of the temperature-dependent thermister resistance. In the 
scanner instrument, the active flake that is attached to the heat sink in- 
directly responds to the incident radiation from the source through its op- 
tical arrangement. The reference flake is also attached to the heat sink 
and is placed where the space is totally sealed by the surrounding materials 
such as the heat sink and housing. Thus it is assumed that the reference flake's 
temperarture is equal to that of the heat sink and its housing. Accordingly, 
the observed difference in resistances between the active and reference flakes, 
is proportional to the incoming radiation. 
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Figure 19 is the block diagram for the scanner control electronics and 
their transfer functions. 

The thermister resistance as a function of temperature is given by 


Rj.(T) 


W EXP 6 i ( ^ 



i - 1. 3 


(67) 


Using the polynominal expansion of the exponential in Eq. 67» the difference 
in resistances between the active and reference flakes is approximated as 


W 


r 3 (t 3 ) 


R lo 6 l T + !3aii 

— f2— T i 


$1 ~ Ro Bq 


t 3 + 


dO. 


( 68 ) 


The constant in the bridge is defined by 


2v B k(l + k) (69) 

K B = R (1 + 2k ) 2 
o 

The transfer function of the bridge is transformed into a differential 
equation as shown in the following 

T 1 *1 + x i " k b (r i " V (70) 

where T is a time constant which, for the case of the scanner is less than 
10 ms. 

For the preamplifier, we end up with the equation 


t 2 * 2 + * 2 - Xj + T 3 - K c C B + 

Thus, the following equations describe 
the differential forms of the transfer 
ized forms of the amplifier and filter 

T = AT + DT 4 + B q Qg + BE 


(E os + n) + T 3 (R os + A) (?1) 

the energy balance equation (Eq. 72) and 
functions, (Eq. 73) which include general- 
and a model of each scanner channel 

(72) 
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A x + 
x 


B T + 


B (E + 
os os 


n) 


+ B C_ + b 
c B 


(7 3) 


where T is a 14 component vector consisting of node temperatures, X Is a 6 
component vector consisting of electrical signal processing variables, Q s is 
the heater heat flow rate, E^, (:„, and b are electrical and flake related vari- 
ables, n is the total noise as viewed at the preamplifier input, and E is the 
incoming radiant incidence vector. The measured quantities which are telemetered 
to ground are the raw counts, x,., and the heat sink temperature T ? , measured by 


a thermister. 



where v is the measurement noise including quantization error due to A/D 
conversion. 

2.1.11 Scanner Steady-State Analysis 

Two types of transient behavior will be noticeable in the instruments. 

One type occurs when the heater Qg is activated, i.e., when the instrument 
is turned on. The transient behavior in this case is rather slow, as the 
heat sink time constant is large and the conductivity of the filters (for 
SW and LW instruments) is low. Thus, the time period necessary 
various node temperatures to reach their relatively steady operating temper- 
atures, particularly the filters, will be long. The precise period can be 

determined by simulation or measurement. 

The second and more important transient behavior is due to quick 
variations in the input (or observed) radiation incidence, E, after the opera 
ting condition has been reached. The time constants, in this case, are 
determined largely by the electronic signal processing parameters and the 
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flake characteristics. After the operating condition has been reached, 
the temperature variations of most nodes should be negligible. The active 
flake temperature then varies with the incoming radiation. The flake 
resistance and other electrical variables also vary with the incoming radiation 
and can produce transient behavior which can be observed if the target level 
is changed suddenly, as a step function. On the other hand, if the incoming 
radiation varies slowly relative to the system time constants, then the system 
will operate in quasi-steady state condition. A preliminary analysis indicates 
that, for the parameters selected by the designer, the applicable time constants 
are in the vicinity of 10 ms. Thus, if the input radiation does not contain 
frequencies above 10 Hz, the transient will be negligible, while frequencies 
above 17 Hz, will^be attenauted. Whereas the earth scan may have low content 
in the high frequencies (which are of less interest to ERBE) , the transition 
from space to Earth may produce a necessarily short period of transient behavior. 
It may be of interest to analyze the magnitude and duration of such transients, 

which probably last less than 50 ms. 

Since the quasi-steady state condition is of interest to this analysis, 

let us first consider one scan period. The scanner looks at space for a period 
"sufficient" for the transient effect to die down and obtains an average of the 
last few values (to reduce noise and errors) as the "space clamp". Then the 
scanner sees the Earth and atmosphere, then looks at two black bodies at ambient 
but known, temperatures as well as the SWICS in the off condition. The kAt 
correspond to the kth space clamp, so that T(kAt) and x(kAt) are the temperature 
and electronic vectors at the space clamp measurement. Now let 

Ax(t) = x(t) - x(kAt) , AT(t) = T(t) - T(kAt), kAt < t < (k+l)At. (75) 
From Eq. 72 and Eq. 73 
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(76) 


AT * AAT + 4D T^ 'AT + Bq AQg + BE 

Ax * A x Ax + B t AT + Ad + B Qs An, 

where it has been assumed that the thermal and radiative properties are in- 
cluded in A, D, Bq and B and the time constants for the electronics in A x , B Qg 

and B have remained unchanged over the scan period of 4 seconds. Also we can 
c 

define Ad, the drift component of the voltage as. 

Ad = AB t T(kAt) + B og AE og + Ab + B c AC g (78) 

Considering the amplifier and Bessel filter to be in a quasi-steady 
state condition, it can be seen that 


Ax, (t) = K f Ax (t) + K ^(t) (79) 

6 t o o 

where A^(t) is the part of Ax 2 (t) with lower frequency content (below 10 Hz) 
consisting of an averaged part of Ax 2 (t). which may be considered the signal 
part, and ^(t) is the part of Ax fi (t) consisting largely of noise and irradiance 
variations over short distances at the top of the atmosphere with zero mean. 

As the filter attenuates high frequencies, Ax fi (t) will have a low amplitude. 

The preamp equations are expressed as 


Ax 2 + 


Ax, 


Ax l + T 3 


Ax x + 


AE + An + T,(AE + An), - K AC 
os 3 OS c 


B 


(80) 


Now assume that 


T 2 E 2 - AE os « 0 


(81) 


then 


Ax 2 - Ax x + t 3 Ax x + AE_ 


os 


K AC + An 
c B 


(82) 
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where Ax^ is the low frequency portion of Ax^, and where it has been assumed 

that AE and AC„ are baseband signals. Note that 
os B 


Ax£ = Ax£ + Ax£, 


(83) 


where Ax 3 contains the higher frequency components of AX 2 and causes the term 

Ax,. The bridge equations are now expressed as 
6 

T 1 A *1 + A *1 = T 1 B T u AT 1 + T 1 B T 13 AT 3 + Ab l + T 1 AB T n T 1 


+ x.AB T + T AB AT + T AB AT 
1 L 13 11 x x 13 


(84) 


It is assumed that the only variable on the RHS containing higher 
frequency components is T^, whereas the remaining variables are slow moving 
drift components. 


Ax l - T 1 (B T„ AT 1 + V AT 3 - ^ 


11 


13 


(85) 


Combining Eqs. 85, 82, and 79, 


’ T 1 (B T u W 1 + B T 13 iT 3> + (T 3 - T l ) ‘ T 2 x 2 


+ E - K AC„ + An 
os c B 


( 86 ) 


ix 6 ■ K f t 1 <B T„ 4T 1 + B T„ “s’ 

o 11 id 

• • 

+ K f r (t 3 - t 3 ) - t 2 x 2 

o L 


+ AE - K AC„ + Si + Ax, I 
os c B 6 J 


(87) 


Now consider the thermal equation 


N 


AT 1 " A 12 (AT 2 - A V + 4 V.h °ik T k AT k + B 1 E 
B 1 E = A 12 (A ^1 " AT 2 J ‘ 4 k=l D lk T k AT k + ^1 


( 88 ) 

(89) 
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1 


(90) 


AT 


1 Kj ®x 

f o 11 


®Ti3 

.Ax, - ^ AT 0 - 


v 6 B 3 T- B 

T 11 1 11 


«J - (T 3 - V ^ - t 2 i 2 + iE os - K c 1C B + 45 + 45 6 
Substituting Eq. 88, 90. and 91 into Eq. 89. we have 


(91) 


B 1 E 


<* 12 - 4 D n T i } 

K f T, B 

f o 11 


N 


Ax, — A, « AT« — 


13 


12 “*2 B 


( A 12 - 4 D H T 1 )AT 3 


11 


_ i AD., T. AT. + AT. r~ 

k=2 ik k k 1 


(A., - A D u V 


11 


(92) 


2.1.12 Scanner Count Conversion 

The scanner count conversion algorithm below is based on steady state 
approximation. The structure (i.e. , the form of the equations) of the algorithm 
is the same for each of the three scanner channels. The parameters have 

different values for each channel- 

The average of scan points during space clamp is given below 

1 n , - (93) 

-< t k>-ni£l V<t ki ) - n ' 8 

The mean variance of the counts compared with the space clamp during a scan 
cycle is obtained by the equation 

<4 ■ i Ji [ v(t k ) ‘ " (c ^ 2 

The scanner count conversion equation is then defined by the following 
equations 

_ -i t ~ ^ 

E(t - T) = A v [v(t) - v(t k )J + A h [t h (0 - T H (t)J + A, 6 k At 


( 94 ) 


fc k-l ‘ ‘ C k 


( 95 ) 
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The coefficients in Eq. 95 are the following 


A 12 " A °11 T 1 


Ja. 


A v'KK t T x B VO 

o 11 




' - { 4 [ D U T l + ° 


12 T 2 + D 13 T 3 + D 14 i + D 15 T 5 + D 16 T 6 


+ D 18 T 8 + D 113 T 13 


A = c 


A 12 ~ 4 D 11 T 1 
KT B 
1 11 


] - + I 

J 11 


Vb o 

-A- = c K f A 

V° £ o V 


{ k ' kT [ I(t k> - '' (t k-i ) ] - T i 


(B T + B t 
A 11 A 13 


> [W - T H (t k-l ) ] 


V ( \i Tl 

o 


+ ®t t ) r v 

13 3 A |_ B 

K d[w - V d (t k-1>] 


(t k> - V 8 (t k 


.>} 


C k " C k-1 + Lt 


( 96 ) 


(97) 


(98) 


(99) 

(ICO) 


In the last term of Eq. 95, (t-t R )/At can be replaced by the number of measure- 
ments in the scan mode as follows: 


t - tk 
At 


= 0.0135 (J-8) 


(101) 


where J is the number of measurements in the scan mode. The constants in the 
above equation are the following: 

v(t R ) = space clamp value at time t R 

v(t ) = instrument output (in volts) when viewing space at t R 
ki 

o 2 = noise variance estimate during space look 
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v (t) = instrument output sample at time t 

T (t) = heat sink temperature measurement at t or most recent value ( K) 

H 

At = total scan period (4 sec.) 

t = time at end of space look (sec.) 
k 

t = sampling instant (sec.) 

V (t) = detector bridge bias voltage measurement at time t or most recent 
B 

value (v) 

v (t ) = drift balance DAC voltage measurement at time t or most recent 
d k 

value (v) 

6 = estimate of unaccounted drift during k th scan period (V) 

k 

T — average time lag (sec.) 


k = post amplification gain 
o 

Note thet to compute (y it is necessary to have the following space clamp 
value; i.e.. vU k >. If the term with 6, in Eg. 95 is neglected, then the count 
conversion may be done without using Wt k >. The variance d* la computed for 
diagnostic purposes and should be output as described below. The converted 
value E ± (t - T), 1 - s,E, T (for shortwave, longwave, and total, respectively) 

can be interpreted as follows 


E.(t) = / T t (X) E x (t) dX + e A (t) 


( 102 ) 


where T^X) is the transmittance of the 1 th channel, E x (t) is the spectral 
irradiance at the instrument due to radiation from the footprint at time t, 
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and e^(t) is the error in the measurement combined with the count conversion 
process, and can be estimated through analysis of calibration data and model 
results . 


2.1.13 Accuracy and Diagnostic Checks 

During each scan, the instruments view two black bodies at ambient, but 
known temperatures. These data can be used to determine if the basic operation 
of the instruments have changed as well as provide an approximate accuracy check 
of the combined error e^(t) in Eq. 101. 

1. Store a table of elements 

(Vg s (Tg) , VpjCT^), v nT (T n ), E Ro (T R ), E r 9 (T b ), E rt (T r )) where T„ “ T 


'BJT B' 


BT B 


Bs B 


-B£ v B' 


BT B 


B 


Bj‘ 

th 


1 £ 1 £ N. The values v (T„) represent the steady state output of the i 

BI B 

channel when viewing the black body at temperature T g , and may be obtained in 

calibration tests. The values E n .(T„) represent the output of count conversion 

Bl B 

when viewing the black body at temperature T fi , given by 


E Bi (T B > = tz t <» e bx (t b ) dA 


(103) 


where E„, (T„) is the spectral irradiance from the black body at temperature T . 
BA B o 


2. Compute E 0 . using the measured black body temperature T R (t) as 
Bi B 


E Bl (T B <t); ' ‘ E B1 + 


T B (t) ~ T B1 
T B2 “ T B1 


[ 


E Bl<V * W. 


(104) 


where T fil £ T fi (t) < Tg 2 

3. Compute E^t) using the count conversion algorithm with the corresponding 
value v^(t) . 


4. 

Plot 

V ± (t) 

and 

v D .(T,,(t)) on same axes. 

ul D 

9. 

Plot 



®Bi^ T B <fe ^ * MMB * onmo. 
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( 105 ) 


6. Compute and plot o^» and 


'v. ,k+l 


a v k + (1 ' a)(v l (t) " V Bi (T B (t))) 


” i = a e ±k + (1 - a)(E i (t) - E Bi (T fi (t))) 


°x^,k+l 3 °v 


0 i,k+l = 3 ~ik 


J k + (1 " a) [ (v i (t) ' V Bi (T B (t)) " 5 v 1 ,fcfJ i 

o], + (1 - a) [<^(0 - E Bi (T B (t)) - i 1>k+1 ] 2 


( 106 ) 

( 107 ) 

( 108 ) 
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3 . 


INSTRUMENT PERFORMANCE AND COUNT CONVERSION TEST 


The instrument models developed for the non-scanner, the scanner, and 
the solar monitor are based upon the theoretical approach with a certain 
degree of approximation. For instance, the energy equations for the instru- 
ments were formulated with the assumptions that the node (or control volume) 
has uniform temperature and homogeneous thermal and radiative properties. 
However, this lumped geometry of the node, in reality, cannot hold the 
uniform temperature through the node geometry unless the temperatures of the 
neighboi nodes are the same, or the node boundaries are perfectly insulated, 
or the thermal conductivity of the node material is infinity. Otherwise, 
if there is any heat flux allowed into or from the node, a temperature grad- 
ient from the center of node to the boundary will exist. 

The radiative properties that are used in the ERBE radiometer instrument 
models may be the values picked out of the handbook or measured and averaged 
values over a certain band of wavelength. No matter where these values were 
obtained, they would introduce some errors when they are used for the compu- 
tation of the radiation exchange between the surfaces. The spectral distri- 
bution of the emissive energy flux from a source of radiation such as MRBB 
may lie on the limited range of wavelength in which the averaged radiative 
properties used may be far above or below the real values within the spectral 
energy distribution. Therefore, the averaged values may not represent 
reality in instrument performance. Also, during long term operation in space 
these properties will be subject to changes due to various reasons such as 
thermal and environmental contamination, and aging effect. 

The electronic systems are used in the radiometer to measure the temper- 
ature difference between the active and reference cavities or flakes, and to 
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digitize the signal, then to transmit it through a telemetry system. The 
parts such as resistors of this electronic system are sensitive to a wide 
range of temperature changes. For instance, if the resistance of the bridge 
slightly changes from the reference point, then the error in irradiance 
measurement due to this temperature change may be in error by a few watts 
per square meter. Although the instrument models are well-conceptualized and 
designed, the real behavior of the radiometers are, as expected, different 
from what is was designed to be. 

Other factors to be considered, which influence the sensor performance 
and measurement accuracy, are the inflight changes from the states that the 
sensor normally operates in and the calibration source variability in its 
operation mode or the changes in source characteristics. An environmental 
change was expected in the long-term operation of the instrument due to the 
orbital changes which caused cyclic variations in the heat load of the in- 
strument. Another change, was due to the contamination or degradation of 
then senosr surfaces due to radiative interactions with other surfaces. 

The internal calibration sources such as MRBB and SWICS are also sub- 
ject to change. The monitoring of the change in intensity of these sources 
by possible variation in power supply and shift of the spectral distribu- 
tion in the source output is possible by using the SiPD, although it is also 
subject to change, and intercomparing calibration sources. 

To account for the possible errors due to the reasons discussed above 
a procedure was developed to detect those errors. The procedure uses the 
algorithm shown in Figs. 20 and 21 for that purpose. This algorithm is 
capable of evaluating the sensor model, adjusting the simulation model para- 
meters for minimum errors, and establishing confidence intervals for the 
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FIGURE 21 





converted data reflecting sensor characteristics and calibration sources. 

The sensor simulation model responses, using known calibration sources 
as inputs to the model, are checked against the instrument calibration out- 
put. Any differences in the slope, offset and curvature between the simula- 
tion and calibration results were interpreted as uncertainties of the model 
simulating the instrument. These differences, however, generally reflect 
the problems discussed earlier. In the algorithm, these differences are 
narrowed by checking and adjusting the sensor model, source, and environmental 
parameters to minimize the errors. However, the estimated parameters must 
stay within the tolerance that can be allowed by the approximations made in 
the modeling or by the nature of reasons described earlier. 

The algorithm consists of three parts (see Fig. 20): 

The first part is a routine to adjust the sensor model parameters by com- 
paring sensor simulation outputs such as counts and HK temperatures to the 
inflight and ground calibrations obtained for the known sources. 

The second part is a procedure to check the errors generated by simulating 
a sensor by comparing them to the calibration data for the known sources. It 
checks the error generated from the sensor simulation target source, or envrion- 
ment at different time intervals for the same known irradiance source as in 
the inflight calibration. 

The third part consists of routines to search for any changes in the target, 
environment and sensor performance due to the degradation of the instrument 
surface, vibration, contaminations, and thermal and aging effects on the radia- 
tion interactive surfaces. 

To do the above job, the algorithm has been designed to have the following 
capabilities : 


67 



1. to adjust the sensor model parameters in comparison with the 
calibration data from known sources to minimize the error be- 
tween the sensor model and the calibration data, 

2. to periodically check any changes in sensor, target and environment 
parameters between inflight calibrations by comparing the calibration 
data for the same sources at different time intervals (e.g., 1 or 2 
week periods) , 

3. to re-define the sensor, target and environment parameters, if 
they have been changed, 

4. to test and check any differences between calibration data and 
responses from the simulated sensor model for the same known 
sources, 

5. to repeat the routines 1 through 4 until the errors produced by 
step 4 are minimized, 

6. to generate the count conversion parameters. 


3.1 Description of Routines 

3.1.1 Routine I 

The calibration sources identified during either ground or inflight cali- 
brations are numbered by j =1, N at time t ± and t i+1 where i - 1, n. The sensor 
simulation model generates temperatures and sensor responses corresponding to the 
calibration sources at a time t^. These temperatures and sensor responses, which 
are obtained from the model simulation, are compared with measured data to com- 
pute errors of the simulation model. Once these errors fall within the pre- 
determined error bounds (after repetitive refinement of simulation parameters), 
then these errors will be used as criteria to determine Lhe measurement errors 
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against known targets on the ground. If not, the parameters for both the tar- 
get source and environment effects will be checked; then in turn the sensor 
parameters will be checked and adjusted, in case when the computed parameters 
are not within the acceptable range based on the materials properties, e.g. its 
optical properties. Such a routine will inlcude a decision block to check targets 
and environmental effects to avoid unnecessary adjustments of sensor parameters 
when the errors may be caused by and come from the target and environment para- 
meters. 

3.1.2 Routine II 

If the errors are within bounds, the sensor responses , j - 1, N at dif- 
ferent times t^ and t^ + ^ will be sorted and evaluated in Routine I to check the 
errors from the sensor simulation model for the targets at different times t. 


3.1.3 Routine III 

After finishing Routine II, using all the target sources at different times, 
t and tj_^, Routine III checks the sensor respones of the same target at a dif- 
ferent time and computes the errors. If the error of the same target at a dif- 
ferent time falls within the error bounds, the Routine will check the next tar- 
get. If not, the history of target source variation will be studied, and at 
the same time the sensor and environmental conditions will be rechecked to deter- 
mine the source of errors. This evaluation in Routine I is performed whenever the 
error signals the flag. The process continues until the error source is recog- 
nized and its correction of error is made by adjusting the related parameters in 
the simulation model. When this Routine finishes satisfactorily evaluating all 
of the target sources, the program will print all the decisions, errors, sensor 
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and calibration source HK data, sensor parameters and counts and will supply the 
final results to the count conversion algorithm. 


3.1.4 Routine IV 

If the error does not fall within the prescribed error bounds, the target 
and environment parameters are checked in order to find whether changes in these 
parameters created pseudo-errors in the sensor responses. If there has been a 
change in either the target or environmental parameters, they are re-defined 
and re-input to Routine I. In spite of repeated checks, if the error cannot be 
minimized, then the thorough review of instruments will be required to determine 
whether the error bounds are to be reassessed. In this routine, when the target 
and environment are checked and relatively free from the erroneous performance 
then we would postulate that the instrument's measurement capability might be 
deteriorated due partially to the thermal and environmental contamination or by 
other means such as aging of the parts and so on. 


70 


A. DISCUSSION OF PARAMETER ESTIMATION 


The simulation code was developed to study the sensor instrument responses 
and behavior during the ground and inflight calibration procedures as described 
earlier in Section 2.1.6. This simulation code uses the irradiance at the FOV 
limiter of the instrument as the input data. All the constants representing 
the instrument geometry, material, thermal and radiative properties and the 
electronic thermal control systems are parameterized to be used as coefficients 
in the system of equations. Then the simulation code generates the counts and 
the nodal temperatures of the model corresponding to the input irradiance. 

When the simulation code uses the input irradiance which is the same as the 
calibration source, the results of counts and housekeeping temperaures from the 
simulation are compared to the calibation results. If the differences from these 
comparisons for various levels of irradiance are wide, then the simulation model 
parameters are adjusted to minimize those differences. This procedure for esti- 
mating the model parameters is repeated until the differences can no longer be 
reduced. These estimations also continue for various calibration sources until 
systematic response patterns for each calibration source are established. All 
procedures mentioned above are performed for each individual radiometer sensor. 
Even though two instruments of the same type are used in the calibration pro- 
cedure, for example, the results may differ due to the inherent differences in 
the many factors involved, e.g. slight differences in the sensor assembly process 
of uniform surface roughness, coatings, glue or contact areas. 

The following is an example that shows how the parameters for the nonscanner 
total channel are selected. Assume that the instrument output derived from 
calibration input irradiances are linear; then the factors which determine the 
slope and offset of the curve for simulated results would be judged to be the 
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offset temperature between the active and reference cavities, and the radiative 
constant for the term of irradiance in Eq. 34, i.e. To determine these two 
factors using the measurement data, the linear form of Eq. 34 can be written 


ax = 3 


(109) 


„here * Is a column vector consisting o£ the offset temperature and radiative 
constant, and a and 6 are defined later. The solution to Eq. 108 is 

_1 (in 

x = (a* a) A a' 3 


where a' is the transpose of a. 

The constant a is a matrix defined by 

o ± j " A 12 + A D 11 (T 2 ) 1’ E 1 

A 12 + * E 2 

. 

. 

A 12 + A D ll^ T 2^i’ E i 

where i is the calibration sequence. 

The constant 3 is a column vector 
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$ as an explicit form is 


3 = - B K V - <f> - B Q 

M 1 a a 1 i J i J i 


(111) 


<j> ■ J D., T?, , only when k = 1, replace T _ by T 

1 if lk ki 


2i 


( 112 ) 


The nominal values are used for unmeasured and except + 5T. 

The temperature of the reference cavity is a little higher than the heat sink 
temperature due to the Joule heating by electric current flowing probe wire 
of reference cavity. The temperature difference, <5T, due to the Joule heating 
cannot be directly measured but can be observed from the simulation results. 

The plots of the calibration and simulation results show that they are 
not actually linear but could be considered to be linear for estimation (see 
for instance, the steady-state response comparison of simulation of calibration 
data for the NS WFOVSW sensor on page 100 in Appendix A). Therefore, the above 
approach can still adjust the overall slope and offset of the curve by regarding 
them as a linear. The linearity and adjustment can be made by adjusting other 
parameters. But since almost all parameters are interrelated, the selection 
of a parameter that mostly influences the linearity of the curve in a plot of 
simulation results must be cautiously taken by testing the linearity with 
various parameters. 

Since the simulation code was developed, the calibration data for the 
nonscanner radiometers using MRBB as a calibration source was the only source 
available. Hence, the nonscanner instrument models have been simulated using 
MRBB data which was used in the calibration. Figure 22 shows the calibration 
data for the MRBB temperatures in the 1st column and the count outputs for 
the four nonscanner channels in the 6th through 9th columns. The 2nd column 
shows the blackbody radiation corresponding to the temperatures in the 1st 
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column. The 3rd and 4th columns show the percentages for the shortwave and 
longwave portions of the blackbody radiation described in the 5th column, 
respectively. These shortwave and longwave portions were calculated by the 
integration of Planck's equation for blackbody radiation with respect to the 
wavelength from 0 to 5 ym for shortwave and from 5 ym to 1000 ym for longwave. 
These portions were used for the MFOV and WFOV shortwave instrument simulations. 
The shortwave range selection of 0 to 5 ym was simply made by the shortwave 
filter (suprasil) characteristics. In the 6th through 9th columns, the numbers 
in the bottom part of the blaock signify the output (in counts) corresponding 
to MRBB radiation from the simulation. The 10th and 11th columns show the 
differences in count between calibration and simulation tests. 
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APPENDIX A 


PRELIMINARY RESULTS ON ERBE NON-SCANNER INSTRUMENT 
MODELING, COUNT CONVERSION AND SENSITIVITY ANALYSIS 


preceding Page blank m*r fumed 

bl ank not mMm 


GOALS 


DEVELOP A DYNAMIC SENSOR MODEL FOR EACH INSTRUMENT 

DEVELOP AN ADAPTIVE COUNT CONVERSION ALGORITHM THAT 
CAN ACCOMMODATE CHANGING CONDITIONS 

DETERMINE THE ACCURACY OF THE ERBE INSTRUMENTS 

DETERMINE THE ACCURACY OF THE ERBE MEASUREMENTS 


emujiNG PAGE BLANK NOT tfILWED 


APPROACH TO ASSESSMENT OF INSTRUMENT ACCURACY 


• DEVELOP A MATHEMATICAL MODEL AND COMPUTER SIMULATION FOR EACH INSTRUMENT 

• COMPARE GROUND CALIBRATION DATA TO CORRESPONDING SIMULATION DATA 

• MINIMIZE DIFFERENCE BY ADJUSTING MODEL PARAMETERS 

• PERFORM A SENSITIVITY ANALYSIS 

• ESTIMATE CONFIDENCE INTERVALS FROM GROUND CALIBRATION DATA 

• REPEAT SAME PROCEDURE FOR FLIGHT CALIBRATION DATA (ERROR ANALYSIS) 

. ADJUST ESTIMATE OF CONFIDENCE INTERVALS AT EACH FLIGHT CALIBRATION 
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APPROACH TO ADAPTIVE COUNT CONVERSION 


DEVELOP A MATHEMATICAL MODEL AND COMPUTER SIMULATION FOR EACH INSTRUMENT 

COMPARE SROUND CALIBRATION DATA TO CORRESPONDING SIMULATION DATA 
MINIMIZE DIFFERENCE BY FINE TUNING MODEL PARAMETERS 

DEVELOP A COUNT CONVERSION ALGORITHM (CCA) BASED ON THE MODEL 
SIMULATE THE INSTRUMENT AND CCA 
ESTIMATE ERROR DUE TO CCA 

REPEAT PROCEDURE FOR EACH FLIGHT CALIBRATION 
FINE TUNE MODEL AND CCA PARAMETERS 
OBTAIN NEW ERROR ESTIMATES 


NON-SCANNER MFOV TOTAL 



7 A1 HEAT SINK 

8 APERTURE 

H HEATER 

P TEMPERATURE PROBE 
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NON-SCANNER WFOVSW 



7 A1 HEAT SINK 


8 APERTURE 

9 SW FILTER DOME 


H HEATER 

p TEMPERATURE PROBE 
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C ‘ ^ 
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OFFSET TEMPERATURE # C 


NON-SCANNER DYNAMIC RESPONSE 


WFOV TOTAL SIMULATION 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 



87 






JURE °C 


NON-SCANNER DYNAMIC RESPONSE 


WFOV FILTERED SIMULATION 


ACTIVE CAVITY 
REFERENCE CAVITY 


APERTURE 





OFFSET TEMPERATURE °C 


NON-SCANNER DYNAMIC RESPONSE 

WFOV FILTERED SIMULATION 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SNK 



0 















OFFSET TEMPERATURE # C 


NON-SCANNER DYNAMIC RESPONSE 


MFOV TOTAL SIMULATION 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 
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RADIATION ERROR W/U 5 SENSOR OUTPUT (COUNTS) RADIATION W/M 
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TEMPERATURE # C TEMPERATURE °C TEMPERATURE 


NON-SCANNER DYNAMIC RESPONSE 


MFOV FILTERED SIMULATION 


ACTIVE CAVITY 
REFERENCE CAVITY 
APERTURE 



- CU HEAT SINK 

- TOP AL HEAT SINK 

- BOT AL HEAT SINK 



26 . 10 | — 


25.88 — 
25.67 — 
25.45 — 
25.23 — 

2£&2r — 

24.80 


1 


50 


100 


150 

TIME (SEC) 


200 


250 


300 
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RADIATION ERROR W/M SENSOR OUTPUT (COUNTS) 


NON-SCANNER DYNAMIC RESPONSE 


MFOV FILTERED SIMULATION 


IRRAD1ANCE 
CONVERTED COUNTS 




RADIATION ERROR 
ZERO LINE 
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STEADY STATE OUTPUT (COUNTS) 


NON-SCANNER STEADY STATE RESPONSE 
COMPARISON OF SIMULATION TO CALIBRATION DATA 

WFOV TOTAL 

O CALIBRATION 
O MODEL 





NON-SCANNER STEADY STATE RESPONSE 
COMPARISON OF SIMULATION OF CALIBRATION DATA 

WFOV FILTERED 

O CALIBRATION 
□ MODEL 
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STEADY STATE OUTPUT (COUNTS) 


NON-SCANNER STEADY STATE RESPONSE 
COMPARISON OF SIMULATION TO CALIBRATION DATA 


MFOV TOTAL 

O CALIBRATION 
□ MODEL 




CALIBRATION - MODEL (COUNTS) 


NON-SCANNER STEADY STATE RESPONSE 


COMPARISON OF SIMULATION TO CALIBRATION DATA 

MFOV TOTAL 
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STEADY STATE OUTPUT (COUNTS) 


NON-SCANNER STEADY STATE RESPONSE 
COMPARISON OF SIMULATION TO CALIBRATION DATA 


MFOV FILTERED 

O CALIBRATION 
O MODEL 




CALIBRATION - MODEL (COUNTS) 


NON-SCANNER STEADY STATE RESPONSE 


COMPARISON OF SIMULATION TO CALIBRATION DATA 

MFOV FILTERED 
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SENSOR OUTPUT (COUNTS) RAWATION W/M 


NON-SCANNER DYNAMIC RESPONSE 


WFOV TOTAL SIMULATION 
TIME DELAY (LAG) ESTIMATION 


RRADIANCE 
CONVERTED COUNTS 



RADIATION ERROR 
UNE 






NON-SCANNER DYNAMIC RESPONSE 

WFOV TOTAL SIMULATION 
TIME DELAY (LAG) ESTIMATION 


ACTIVE CAVITY 
REFERENCE CAVITY 
APERTURE 



BOT AL HEAT SINK 





OFFSET TEMPERATURE 


NON-SCANNER DYNAMIC RESPONSE 


WFOV TOTAL SIMULATION 
TIME DELAY (LAG) ESTIMATION 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 




RADIATION ERROR W/M 2 SENSOR OUTPUT (COUNTS) RADIATION W/M 


NON-SCANNER DYNAMIC RESPONSE 


MFOV TOTAL SIMULATION 
TIME DELAY (LAG) ESTIMATION 


RRADIANCE 
CONVERTED COUNTS 



RADtAHON ERROR 
ZERO LINE 









OFFSET TEMPERATURE °C 


NON-SCANNER DYNAMIC RESPONSE 


MFOV TOTAL SIMULATION 
TIME DELAY (LAG) ESTIMATION 


ACT CAV - RET CAV 
ACT CAV-CU HEAT SMK 
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SENSOR OUTPUT (COUNTS) RADIATION W/M 


NON-SCANNER DYNAMIC RESPONSE 

WFOV FILTERED SIMULATION 
TIME DELAY (LAG) ESTIMATION 


^RADIANCE 
CONVERTED COUNTS 



RADIATION ERROR 
ZERO LINE 











OFFSET TEMPERATURE °C 


NON-SCANNER DYNAMIC RESPONSE 


WFOV FILTERED SIMULATION 
TIME DELAY (LAG) ESTIMATION 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 




RADIATION ERROR W/M 2 SENSOR OUTPUT (COUNTS) RADIATION W/M 


NON-SCANNER DYNAMIC RESPONSE 


MFOV FILTERED SIMULATION 
TIME DELAY (LAG) ESTIMATION 


RRADIANCE 
CONVERTED COUNTS 
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NON-SCANNER DYNAMIC RESPONSE 


mfov filtered simulation 

TIME DELAY (LAG) ESTIMATION 


— ACTIVE CAVITY 

REFERENCE CA 

APERTURE 

SW FILTER 







OFFSET TEMPERATURE *C 


NON-SCANNER DYNAMIC RESPONSE 


MFOV FILTERED SIMULATION 
TIME DELAY (LAG) ESTIMATION 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SMK 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY EMISSIVITY CHANGED BY .02 (NOMINAL .92) 


^RADIANCE 
CONVERTED COUNTS 



8200 
8050 — 



RADIATION ERROR 
ZERO UNE 
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SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

ACTIVE CAVITY EMISSIVITY CHANGED BY .02 (NOMINAL .92) 


ACTIVE CAVITY 
REFERENCE CAVITf 



BOT AL HEAT SNK 
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SENSITIVITY ANALYSIS USIN6 DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY EMISSIVITY CHANGED BY .02 (NOMINAL .92) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SNK 



SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

REFERENCE CAVITY EMISSIVITY CHANGED BY .02 (NOMINAL .95) 


IRRADIANCE 
CONVERTED COUNTS 




RADIATION ERROR 
ZERO UNE 
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TEMPERATURE °C TEMPERATURE °C TEMPERATURE ®C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

REFERENCE CAVITY EMISSIVITY CHANGED BY .02 (NOMINAL .95) 


ACTIVE CAVITY 
REFERENCE CAVITY 
APERTURE 



- BOT AL HEAT SINK 


10.05 


K >.00 


59.95 


59.90 


59.85 


59.80 


39.75 


BASE PLATE 
FOV UMITER 
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OFFSET TEMPERATURE *C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

REFERENCE CAVITY EMISSIVITY CHANGED BY .02 (NOMINAL .95) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

DOME FILTER LONG WAVE TRANSMITTANCE (LEAK) CHANGED BY ,005 (NOMINAL .01) 


■ IRRADIANCE 

CONVERTED COUNTS 

950.00 

CM 

791 .67 

g 635.33 

F 

< 475.00 
O 

< 

O: 316.67 

158.33 
0 

8200 
8050 
7900 
7750 
7600 
7450 
7300 




RADIATION ERROR 
ZERO UNE 
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TEMPERATURE # C TEMPERATURE °C __ TEMPERATURE 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

DOME FILTER LONG WAVE TRANSMITTANCE (LEAK) CHANGED BY .005 (NOMINAL .01) 


ACTIVE CAVITY 
REFERENCE CAVITY 



Ctl HEAT SINK 
TOP AL HEAT SINK 



BASE PLATE 
FOV UMITER 
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OFFSET TEMPERATURE 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV FILTERED SIMULATION 

DOME FILTER LONG WAVE TRANSMITTANCE (LEAK) CHANGED BY .005 (NOMINAL . 01 ) 


KRMMANCE 
CONVERTED COUNTS 




RADIATION ERROR 
ZERO UNE 




SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV FILTERED SIMULATION 

FILTER LONG WAVE TRANSMITTANCE (LEAK) CHANGED BY .005 (NOMINAL .01) 


ACTIVE CAVITY 

REFERENCE CAVITY 

APERTURE 

SW FILTER 







OFFSET TEMPERATURE *C 


SENSITIVITY ANALYSIS USIN6 DYNAMIC MODEL 


MFOV FILTERED SIMULATION 

DOME FILTER LONG WAVE TRANSMITTANCE (LEAK) CHANGED BY .005 (NOMINAL .01) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SMK 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

DOME FILTER SHORT WAVE TRANSMITTANCE CHANGED BY .005 (NOMINAL .946) 


950.00 

.791.67 


5 633-33 

O 

I— 

< *475.00 
O 

< 

cc 316.67 


158.33 

0 


It 

I 

S 


IR RADIANCE 
CONVERTED COUNTS 



RADIATION ERROR 
ZERO LINE 
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TEMPERATURE # C TEMPERATURE °C TEMPERATURE °C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

DOME FILTER SHORT WAVE TRANSMITTANCE CHANGED BY .005 (NOMINAL .946) 


active CAvrrr 

REFERENCE CAVTTV 
APERTURE 



CU HEAT SNK 
TOP AL HEAT SNK 
SOT AL HEAT SNK 


BASE PLATE 

" v limiter 



160 

TIME (SEC) 
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OFFSET TEMPERATURE 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

DOME FILTER SHORT WAVE TRANSMITTANCE CHANGED BY .005 (NOMINAL .946) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 
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TEMPERATURE °C TEMPERATURE °C TEMPERATURE °C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV FILTERED SIMULATION 

DOME FILTER SHORT WAVE TRANSMITTANCE CHANGED BY .005 (NOMINAL .946) 


ACTIVE CAVITY 

REFERENCE CAVITY 

APERTURE 



CU HEAT SNK 
TOP AL HEAT 9NK 
BOT AL HEAT SNK 




133 



SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV FILTERED SIMULATION 

DOME FILTER SHORT WAVE TRANSMITTANCE CHANGED BY .005 (NOMINAL .946) 


IRRABANCE 
CONVERTED COUNTS 




RADIATION ERROR 
ZERO UNC 
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OFFSET TEMPERATURE °C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV FILTERED SIMULATION 

DOME FILTER SHORT WAVE TRANSMITTANCE CHANGED BY .005 (NOMINAL .946) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 
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SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV FILTERED 
DOME FILTER LONG WAVE EMISSIVITY 


SIMULATION 

CHANGED BY .005 (NOMINAL .69) 


iRRADtANCE 
CONVERTED COUNTS 




RADIATION ERROR 
ZERO JNE 



J 

300 
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TEMPERATURE *C TEMPERATURE °C TEMPERATURE e C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV FILTERED SIMULATION 

DOME FILTER LONG WAVE EMISSIVITY CHANGED BY .005 (NOMINAL .69) 


ACTIVE CAVITY 
REFERENCE CAVITY 
APERTURE 



CU HEAT SNK 
TOP AL HEAT SINK 
BOT AL HEAT SINK 
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OFFSET TEMPERATURE *C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV FILTERED SIMULATION 

DOME FILTER LONG WAVE EMISSIVITY CHANGED BY .005 (NOMINAL .69) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SNK 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


mfov total simulation 

SENSOR BRIDGE VOLTAGE CHANGED BY .01 v (NOMINAL 10 v) 


RRADIANCE 
CONVERTED COUNTS 



8200 

8050 

7900 

7750 

7800 

7450 

7300 


I: 


RADIATION ERROR 
ZERO UNE 
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TEMPERATURE # C . TEMPERATURE °C TEMPERATURE 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

SENSOR BRIDGE VOLTAGE CHANGED BY .01 v (NOMINAL 10 v) 


ACTIVE CAVITY 
REFERENCE CAVITY 
APERTURE 


O 



40.4 


40.3 


40.2 


40.1 


- CU HEAT 9NK 

- TOP AL HEAT SNK 

- BOT AL HEAT SNK 


10.05 

10.00 


19.95 


19.90 


59.85 


59.80 


59.75 


BASE PLATE 
FOV UMITER 
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OFFSET TEMPERATURE 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

SENSOR BRIDGE VOLTAGE CHANGED BY .01 v (NOMINAL 10 v) 


ACT CAV - REP CAV 
ACT CAV-CU HEAT SINK 
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TEMPERATURE °C TEMPERATURE °C TEMPERATURE *C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY HEATER RESISTANCE CHANGED BY .1 fl (NOMINAL 1492.3 Q) 


- ACTIVE CAVITY 

- REFERENCE CAVITY 

- APERTURE 


40.6 

405 

404 

405 
402 



GU HEAT 9NK 
TOP At HEAT SINK 
BOT AL HEAT SINK 


1005 

♦OOO 

30.05 


30.00 


5|23| 



BASE PLATE 

FOY UMITER 
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OFFSET TEMPERATURE °C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY HEATER RESISTANCE CHANGED BY .1 ft (NOMINAL 1492.3 fl) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .01 0 (NOMINAL 2319.1 fl 8 40*C) 


950.00 

M 

791.67 
z <1 33 ,33 
<47M0 
2 316.67 

166*35 

0 


f 




IMUOIANCC 
CONVERTED COUNTS 



RADIATION ERROR 
ZERO EMC 



J 

300 
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SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

ACTIVE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .01 fl (NOMINAL 2319.1 ft 8 40*C) 


40.6 


O 

° 40.6 

tl 

UJ 

or 

P40.4 


— 

< 


£]4Q.3 

— 

CL 


jjj 40.2 

— 


40.1 


- ACTIVE CAVITY 

- REFERENCE CAVIT 

- APERTURE 


40.0 


CU HEAT ShK 
TOP AL HEAT SNK 



BASE PLATE 
FOV A4TCR 
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OFFSET TEMPERATURE a C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .01 ti (NOMINAL 2319.1 fl 9 AO*C) 


ACT CAV - REF CAv 
ACT CAV-CU HEAT SNK 
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SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

REFERENCE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .01 0 (NOMINAL 2327.1 fl 8 40°C) 


IR RADIANCE 
CONVERTED COUNTS 




RADIATION ERROR 
ZERO UNE 




SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

REFERENCE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .01 0 (NOMINAL 2327.1 (1 S 40*C) 


40.6 

— 

O 


° 4a5 
u 

r_ 

S 

p 40.4 

— 

< 


“ 40.3 

— • 



^4o.a 

ATS « 


W>1 

40.0 



ACTIVE CAVITY 
REFERENCE CAVITY 
APERTURE 


CU HEAT SINK 


TCP AL HEAT 9N< 

BUT AL HEAT 9N< 



BASE PLATE 
FOV UMfTER 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

SENSOR BRIDGE RESISTANCE CHANGED BY .01 ft (NOMINAL 2210 Q) 


930.00 

M 

\ 791 .67 
3E 

z 433.33 
O 

< -473X10 

$ 

cz 316.67 

15833 

0 







WRABANCE 
CONVERTED COUNTS 



RADIATION ERROR 
ZERO UNE 
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TEMPER A TU RE °C TEMPERATURE °C TEMPERATURE 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

SENSOR BRIDGE RESISTANCE CHANGED BY .01 fl (NOMINAL 2210 0) 


ACTIVE CAVITY 

REFERENCE CAVITY 

APERTURE 



TOP At HEAT 9NK 
BOT AL HEAT SNK 



BASE PLATE 
FOV LMTER 
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OFFSET TEMPERATURE *C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

SENSOR BRIDGE RESISTANCE CHANGED BY .01 fl (NOMINAL 2210 fl) 


ACT CA V - REF CA V 
ACT CAV-CU HEAT SINK 
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SENSOR OUTPUT (COUNTS) 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

ACTIVE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .5 0 (NOMINAL 2319.1 0) 


moo 

f 

791.67 

3s 

Z (33J] 
O 

< 47SM) 



2 i 

0£ 316.67 
136J3 
0 



CONVERTED COUNTS 



RADIATION ERROR 

ZERO LINE 
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SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY SENSOR UIRE RESISTANCE CHANGED BY .5 £2 (NOMINAL 2319.1 (2) 


40.6 

O 

° 403 
Ul 

®540l4 

£403 

CL 

m402 




ACTIVE CAVITY 
REFERENCE CAVITY 
APERTURE 


40.1 


4ao 


CU HEAT 9NK 
TOP AL HEAT SINK 
BOT AL HEAT SthK 



BASE PLATE 
FOV UMfTER 
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SENSITIVITY ANALYSIS USING DYNAMIC MODEL 

MFOV TOTAL SIMULATION 

ACTIVE CAVITY SENSOR WIRE RESISTANCE CHANGED BY . 5 £1 (NOMINAL 2319.1 8) 


ACT CAY - REF CAV 
ACT CAV-CU HEAT 3NK 





200 


TIME (SEC) 







OFFSET TEMPERATURE # C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 
REDUCED CONDUCTANCE BETWEEN NODES 1 AND 2 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 
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SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


MFOV TOTAL SIMULATION 

REDUCED CONDUCTANCE BETWEEN NODES 1 AND 2 
» ! 

ACTIVE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .5 ft (NOMINAL 2319.1 ft) 


§ 

§ 

2 

CC 


950.00 

791.67 
633J3 

475.00 

316.67 
158.33 


H 


i\ 


\ 


/ ^ 




tRRADIANCE 
CONVERTED COUNTS 


W 821 

5 

ZD do 

o 

o 


3 

3 

° 7 

on 

o 

00 7 

z 

ill 


RADIATION ERROR 
ZERO UNE 



C 
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TEMPERATURE ®C TEMPERATURE *C TEMPERATURE *C 





OFFSET TEMPERATURE *C 


SENSITIVITY ANALYSIS USING DYNAMIC MODEL 


* 


MFOV TOTAL SIMULATION 
REDUCED CONDUCTANCE BETWEEN NODES 1 AND 2 
ACTIVE CAVITY SENSOR WIRE RESISTANCE CHANGED BY .5 tl (NOMINAL 2319.1 fi) 


ACT CAV - REF CAV 
ACT CAV-CU HEAT SINK 
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ORIGINAL page b 

OE POOR QUALITY 
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CONCLUSIONS 


DYNAMIC MODELS HELP THE UNDERSTANDING OF THE ERBE INSTRUMENTS 

SIMPLE COUNT CONVERSION ALGORITHMS HAVE BEEN DEVELOPED 

STEADY STATE SIMULATION AND GROUND CALIBRATION RESULTS ARE IN 
GENERAL AGREEMENT 

THE ERRORS DUE TO COUNT CONVERSION ARE BEING ANALYZED BY SIMULATION 

INITIAL SENSITIVITY AND ERROR ANALYSES PROVIDE CONSTRUCTIVE RESULTS 

OVERALL APPROACH SEEMS WELL-SUITED TO ASSESS INSTRUMENT PERFORMANCE, 
ACCURACY AND VALIDATION 


PftSClDINO PAGE BLANK NOT fK4K 
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APPENDIX B 


ERBE SCANNER SIMULATIONS 




prscedinq Page 


PLANT 



These graphical displays were produced during development of the simula- 
tion model. Please note that the plot of "OUTPUT (CNTS)" versus "TIME (SEC)" 
in the graph on pages 169, 173, 174, 177 and 180 does not display quantita- 
tively absolute values for counts generated from the simulation model due to 
uncertain gain/offset parameters in the electronics. 


ERBS NFOV TOTAL CHANNEL 


£ 

</> 


5 


346.982906 
247.644933 
1 46.70698 0 
49.568987 
-49.566987 L 



COUNTS 


</> 


z 

o 


3 

CL 


3 

O 


3326.132711 
2721.381309 
2116.629907 
1511.875505 
907.127103 
302.375701 
-302.375701 L_ 




ACT1\€ ruKC T1 
RCFEREHCt FIAKC T3 


O 

o 

Cl 

2 

u 


47.148068 
47.136412 
47.124756 
47.113100 
47.101 444 
47.089788 
47.078132 




ACT1V€ TWICE Tt 


47.124045 
47.117364 
O 47.1 10683 
° 47.104002 
^ 47.097321 
U 47.090640 f- 
F 47.083960 L 



RcracNcc flakc ta 


47.144667 
47.139*13 
a 47.1 3435 i 
° 47.129105 
47.1 23851 
U 47.1 18597 f- 
F 47.11 3343 L 



T1 - T3 


O 


CL 

2 

U 


-.020823 

-.022249 

-.023676 

-.025103 

-.026530 

-.027957 

-.029383 



0 55.01 


110.02 



165.03 220.03 275.04 330.05 

TIME (SEC) 
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38.010796 
38.007507 
0 38.004217 
° 38.000927 
^ 37.997637 
UJ 37.994347 
p 37.991057 



DC TEC TO# ASSDfiLY T2 


36.065236 
38.081709 
<j 38.0581 83 
° 38.054657 
^ 38.051 131 

U 38.047605 
F- 36 044078 



PRCQ& APERTURE H 


36.078190 
38.072219 
<j 38.066248 
° 38.060277 
^ 38.054306 

U 38.048335 
r- 38.042365 



PRIMARY NSOH T5 


31.849668 
31.078157 
0 30.306647 
° 29,535136 
^ 28.763626 

LJ 27.9921 15 
r* 57 750605 



PRIMARY MPROR TS 


31.728329 
34.622611 
0 34.517293 
° 34.411775 
^ 34.306257 
U 34.200740 

r" 34 0Q5555 



SECONDARY MRftOR TO 


36.855956 
36.777971 
0 36.699986 
° 38.622002 
^ 36.54401 7 

U 36.460032 
F- 36.388046 



165.03 220.03 

TIME (SEC) 







ER3S NFOV TOTAL 


a: 

i — 

on 


% 545.258853 
446.1 20879 
£ 346.952906 
2 247.644933 
2 148. 70S 960 
< 49.56S987 

^ -49.568967 

a: 



1031.968569 

00 1580.701557 
2 1229,434544 
§ 878.167532 
a: 526.900519 
CJ 175.633506 
y -175.633506 
UJ 




.251737 
(/y .205967 
2 .160196 
§ .114426 
.068556 
O .022885 
y -.022885 



UJ 


1440.705236 
^ 982.203458 

y 523.611680 
§ 65.019902 

£ -393.571876 
O -852.163653 
gjUJ -1310.755431 

< 

.251738 
.205907 
y .160197 
O .1 1442*“ 
<£ .09o056 
U .022685 
y -.022885 
Ld 



0 55.01 110.02 165.03 


TIME (5EC) 


MT»1 T1« 


153.944516 ^LlO* 7 

33.430035 
g -57.054447 
-207.598929 
U -328.113410 
X -448.627892 K 
-569.1 42374 L_ 



HIM Tii 



in 


z 

o 


3 

£L 


3 

O 


3326.132711 

2721.381309 

2116.629907 

1511.878505 

907.127103 

302.375701 

-302.375701 




CCA RAO 



545.2588531 
446.1 20F79 
346.982906 
g 247.844933 
o: 1 48.706950 
g 49.568987 
-49.568987 



55.01 


± 


1 10.02 


± 


_L 


165.03 220.03 

TIME (SEC) 


275.04 


330.05 
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RADIATION W 






ELECTRONICS 



fftMMANCE 



to 6*5271 
H -6.5271 



active rv*t Ti 

KCFERCNCC rVMCC T3 


DCTCCTOA ASSEMBLY T2 

pRcas A«imj»E ta 


O 

d 

It 

Ui 


< 

5 

5 ! 

H* 

O 


O 

U 

Z 

V) 

m 

o: 





177 



HEATER 











• 89.1737 

* -89.1737 
t -287.521 1 
5 -445.8684 
q. -624.2158 
5-6025632 
y -980.9105 


q; 39.0526 
3 36.9474 
< 34.8423 
£ 327372 
^ 30.6321 


q; 38.1857 
3 37.8143 
5 37.4429 
£37.0716 
0- 36.7002 





HTftl TM 
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n 0521 
sy -.052 
0-156 1 
^ -.260 

-J-.466 
w -.572 






TIME (SEC) 
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APPENDIX C 


ERBE NON-SCANNER THERMAL MODEL 
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NON- SCANNER MODEL 


The energy balance equation 


n T. - T. 

= y _j k 

k j-i 


**k T k , L ^ + \ + ^source 


( 1 ) 


The second term in the right-hand side of the above equation is the 
radiation interaction of the k^ h element to the other elements in an en- 
closure. This variables are defined by the following approach. 

For the radiant interchange in an enclosure having both specular and 
diffuse surfaces, the surface properties are defined as the followings: 

The reflectivity is 

d S 

P 3 = P 8 + P S 

d s 

p s = P l + P l 

where subscripts, s and l, signify short and long waves and superscripts, 
d and s, denote diffuse and specular. 

For diffuse and gray approximation, the following realtionships are 

obtained 

p + e + X =1 
^s s s 

p £ € .e T e = 1# 

From now, the subscripts for short and longwaves are omitted for brevity. 


The radiosoty at the k h surface is 

M? = e,aT? + p. d E. + X. E. . 
k k k k k kk 

The incident radiant flux at the k ^ surface is 

E, = l F .M d 

k h k-j j 


( 2 ) 


( 3 ) 
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Combining thu above two equations yields 


d- ,-l 


M j ' | (S ji " P j F j-i ) <l k fi i + T i E i ) 


( 4 ) 


From (3) and (4), 


where 


E. = T(f .fl. + T.E.) y F, , P. 
k 4 i i 11“ k-j J-1 

i 3 


P j-i " (6 ji ‘ P i ¥ 3-± ) 


(5) 


n 


a T 4 . 


The radiant engergy balance equation is the difference between the 
incident energy and the radio so ty. That is. 




" E k ' 'k n k " p k E k p k E k T k \ 


Substituting (5) into (6) we have 


(6) 


♦k - \ la - p k> I Vj p j-i - % V n i 


Let 


H. 


+ l K 1 - P k ) T i 


= y f p 

k-i j k-j j-i 


V 5 


4 r k-j P j-i T k 5 ki-* E i 


(7) 


Then, for the cavity a, the equation (5) can be 
expanded as 

E 10 ■ H lo-l 'l n i 

+ [H 10-8 e 8 + H 10-9 *9 ] Q 8 


+ H 10-10 € io n io 



+ H 10-10 T 10 E 10 


(A) 
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For the cavity, b 


E 10 " H 10-5 £ 5 S + H 10-6 6 6 P “6 + H 10-10 £ 10 R 10 


+ H 10-10 T 10 E 10 + H 10-ll T n e fov 


From (A) and (B) , 


E 10 " G £ 1 H 10-l p ‘i 


+ G < 5 T 10 » 10 _ 5 H 10 _ 10 


+ G ( 6 T 10 H 1q _ 6 H 10 _ 1q fi 6 


+ G ^ H iO-8 ( 8 + H 10-9 *9^ ^8 


+ (3 t, n [H' 


10 lH 10-10 + T 10 H 10-10 n 10-10 J “10 


H 


J n. 


+ G T 10 T 11 H 10-10 H 10-ll e fov 


2 a \ 1 

inhere G » (1 t 10 H 10-10 h 10-l(r * 


E 10 - G € 1 T 10 H 10-10 H 10-1 P 1 


+ G £ 5 H 10-5 P 1 


+ G £ 6 H 10-6 P 6 


+ G T 10 H 10-10 [H 10-8 € 8 + H 10-9 ‘ 9 1 P “8 


+ G £ 10 H io-io [1 + T io H 10-10 ] P 10 


+ G T H 10 -n e fov 


♦t ■ ? [(l - Pi> H iJ ‘j - 4 U 'j J "3 + ^ t(1 ' p i ) "ij T j - 4 ij T j 1E j 


In cavity a: 
.a 


<1>1 = _ + ^ p i^ T 10 G H l-10 H 10-10 H 10-l 


h,« , - 1] n. 


+ (1 - P x ) < 5 T 10 G H 1 _ 10 « 10 _ 5 

+ (i - p i> * 6 t 1q G ^^0 ^0-6 fi 6 

+ (1 - P 1 ) {((g H 1 _ 8 + e 9 H^g) + t 1q G Hj_ 10 H 10 -10 8 H 10-8 


+ H 10-9 € 9 )J 


+ (1 - p^ « 10 ® 1 _ 10 1 1 + T !o G H 10-10 (1 + T 10 H 10-10 )] 


10 


+ (1 - P x ) T 1q T n G H 1 _ 1() H 10 _ 11 E fqv 


.a A 8 . . A 9 . 

^8 = IT" “ ^8 + *9 


(8+9) 


(8+9) 


^8 = (1 " P 8 )tf l H 8-l + € 1 T 10 G H 8-10 H 10-10 H 10-l^ °1 

+ (1 - p 8 ) « 5 T 10 G H g _ 10 H 1q _ 5 fi 5 

+ d - Pg) « 6 x 10 g h 8 _ 10 h 10 _ 6 n 6 

+ {< s t(l - P a ) H r _ 8 1] + (1 p g ) t 1q G Hg_ 10 H 10 _ 10 ^ H 10-8 '8 


8 u '* K 8 ' 8-8 

+ H 10-9 *9* + (1 " 


8) € 9 H 8-9* ^8 


+ (1 - pg) e 1Q H g _ 10 [1 + T 1f) G H 1q _ 10 (1 + T 10 H 10 _ 10 ^ 


+ (1- Pg) t 1q T U G Hg_ 10 H iQ _ 1;l E fqv 
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(Jjg = (i - Pg) + T 10 G H 9-10 H 10-10 H 10-l ^ ^1 

+ (1 - P 9 > ( 5 t K) G H 9-10 H 10-5 S 
+ (1 - p g ) « 6 T 10 G H 9-10 H 10-6 

+ {(1 - Pg) €g Hg_g + £g[(l - Pg) Hg_g ~ l] + (l ~ P g ) G Q ^Q-lO 

(H 10-8 *8 + H 10-9 e 9)}n g 

+ (i - pg) ( 1Q Hg_ 1() [l + T io G H lO-10 (1 + T 10 H 10-10^ fi 10 

+ (1 - pg) t 1q T u G Hg_10 H 10 _ n E fqv 

In Cavity b: 

^5 " (1 P 5 J *1 T 10 G H 5-10 H 10-l 

+ € 5 “ ^ 5 ^ H 5-5 1] + (1 P 5 ) T 10 G H 5 _iq H iO-5 H 10-10* S 

+ (1 - P 5 ) € 6 t H 5 _ 6 + T 10 G H 5-10 H 10-6 H 10 - 10 ] ^6 

+ (1 - P 5 ) t 1q G H 5 _ 1q (Cg h 10 -8 + 6 9 H 10-9^ ^8 

+ (1 - P 5 ) e 10 H 5 _io ^ + T 10 G H 10-10 ^ + T 10 H 10-10 )] “ Q 10 

+ (1 - P 5 ) t iQ x u G H 5 _ 1q H 1q _ 10 U l0 _ n E pov 

^6 = f l (1 " P 6 ) T 10 G H 5-10 H 10-l Q 1 

+ ^ P 6^ f 6^ H 6-5 + T 10 G H 6-10 H 10-5 H 10-10-* S 

+ e 6 {[(1 - p 6> H 6-6 - 1] + (1 " p 6 } T l 2 0 G VlO H 10-5 H 10-10 } ^6 

+ (i - P 6 ) T l0 G h 6 _ 10 ( ( 8 H l0-8 + € 9 H 10-9^ 9 '8 

+ (i - P 6 ) « 1Q h 5 _ 10 [1 + T io G H 10-10 (1 + T 10 H 10-10 )] fi 10 

+ (i - p 6 ) T l0 T 11 G H 5 _io H 10-10 H 10-ll E FOV 
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♦Jo - «1 T 10 Gt(1 - Pio> H io-io - 11 H 10-l 

+ f 5 H 10-5 " P ]V + T 10 G ^ 1 ■ H 10-10 ~ ^ H 10-10 } S 

+ 4 6 H 10-6 ~ P 10^ + T 10 G ^ 1 " P 10^ H 10-10 ” H 10-10^ ^6 

+ T 10 G ^ 1 “ p lO^ H 10-10 ~ 1 '* ^8 **10-8 + ( 9 **10-9^ fl 8 
+ 1 10 [ (1 - P 1Q ) H 10 _ 10 " 1] U + t 1q G lH 10 _ 10 + T iQ H 10 _ 10 H io-10^ ^10 

+ T 10 T 11 G [ (1 " P 10 ) H 10-10 " H 10-10 H 10-11 e fov 

. ^Oa .a A 10b .h 

ho filter - A 10(a+b) *10 A 1Q(a+b) *10 

Cavity C: 

<j[(i - p x ) h 1 _ 1 - i] o 1 + (i - Pj) s 2 h 1 _ 2 n 2 

+ w - P x ) *4 H !^ n 4 

*2- » - P 2 > ‘l H 2-l R 1 + ‘2 [(1 - P 2 ) H 2-2 ' ^ “ 2 

+ (1 - P 2 ) « 4 H 2 _ 4 C! 4 

(1 - P 4 > <! Vl % + (1 - P*> <2 H 4- 2 °2 
+ * 4 I<i - p 4 ) h 4 _ 4 - i] n 4 

c c 

where e and p are e and p . 

Cavity dr 

*2 - 4 2 [(1 - p 2 > h 2 . 2 - U a 2 + (1 - P 2 ) *3 h 2 _ 3 n 3 
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<t>3 = (1 _ P3) ( 2 ^3-2 ^2 + ~ ^3^ ^3-3 ~ ^ ^3 

+ (l - p 3 ) « 7 h 3 _ ? 

4> 7 = (1 - p 7 > 1 2 ®7_2 ^2 + ^ ^7^ ( 3 ^7-3 ^3 

+ c ? [ Cl - p 7 ) h ? _ 7 - 1] 

Cavity e: 

4* - ^[(l - P x ) « 1 _ 1 - 1] ftj + (1 - P x ) * 2 H x _ 2 fi 2 
+ (1 - p^) ®i-3 ^3 

<l >2 - <1 - P 2 ) «! Vi n i + € 2 [(1 " p 2 } H 2-2 ' 1] n 2 

+ (1 - P 2 ) 6 3 H 2_3 ^3 

<J >3 = C 1 “ P 3 ) £ i H 3 _i + ^ p 3^ * 2 H 3-2 ^2 

+ « 3 [(1 - P 3 ) H 3 _ 3 -1] ^ 


Cavity f : 

d>3 = € 3 [(1 - P 3 ) h 3 _ 3 - 1] fl 3 + (1 - P 3 ) f 1 h 3 _ 7 ft 7 

<J > 7 = (1 - P 7 > « 3 H 7 _ 3 ^ 3 + e 7 [(l “ P 7 ) H 7_7 “ ^ ^7 

The energy balance equation (1) is written for each node as the 
followings : 
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T 4 = 


T 5 = 


, T -T. T,-T- T.-T 0 T_-T, 

1 r 1 2, 4 2 , 3 2, 7 2 IA e A e,„i 

m“ + “5 + ~n + “r + A 4> 2 + Q } 

M 2 R 2-l R 2-4 K 2-3 K 2-7 1 L 


1 r T 2 T 3 o. A e A e j. * d A d 4- A f * f x A 1 

Mj { "r^ + A 3 *3 + A 3 *3 + A 3 ^3 + Q H 3 } 


i T 0 -T. T.-T. T_-T. 

1 / 2 4, 5 4, 8 4 , A c A c , ^ \ 

\ { ^ + ^ + ^7 + a ^ 4 + 94) 


, T. -T_ T,-T T --T,. , , 

1 r 4 5 , 6 5, 12 5 , A b A b-, 

^ + ^ + A ^’ 


T -T 

1 T 5 6 4 . A b A b l 

+ A 6 v 


T r^ 

7 A 2-7 

T,-T„ T,-T n T, rt -T„ 




1 r 4 *8 . ‘5 ‘8 , *10 ‘8 , A a A a , A a ^ 

8 ' mT {_ rTT “CT “5777 8 *8 + a 9 V 


x 

1 /"k 


10 M, 


•8-4 

R 8-5 

R 8-10 

T -T 
4 10 

^ V T 10 A a 

Ri rk / 

R-. n 

o T "10 


Where M = V y C. In addition there are three more equations which were 
derived from the electronics model. These equations are already discussed 
in Section II. 

In the above equations, the <J>s can be splitted as 


^k shortwave + ^k^ longwave " 


Substituting into the above equations and rearranging the constants 


yeild a following type of equation 

3 ~ ~ 3 


T. = a A , • T . + A D , • T 4 + B 0. + B, E. 
k ‘ k] j . kj j Qk k k 


FOV 
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where 


^ ■ Wi 


D * constants including all radiative parameters 

kj 


\ 


Qk ^ 


B * constants related with the irradiance at FOV limiter of 
k 

radiometers . 

The above constants are described in the following pages. 


{A. 


i-l M lRl _ 2 • 


1-2 


Vl-2 


■2-i M 2 R 2 . 1 


‘2-3 M 2 R 2 _3 


1 _ . 1 , _ 1 _ . 1 + 1 % 
2 -2 ' ' « 2 R 2-l R 2-4 R 2-3 R -27 


2 -‘ «2 R 2-4 


’ S-7 


2-7 M 2 R 2 . 7 


{A 


3-2 M 3 R 3 _ 2 ’ 


'3-3 M 3 R 3 _ 2 


4-2 M 4 R 4 _ 2 


A - 1_ ( l + _J_ + __±_) 
A 4-4 - - M 4 \_ 2 R 4 _ 5 R 4 _3 


4 - 5 M 4 R 4-5 


M 4 R 4-8 


'5-4 M 5 R 5 _ 4 


v = _i_ ( _I_ + - 1 - + -1_) 

55 M 5 R 5 _ 4 R 5 _ 6 R 5_i2 


5“6 M 5 R 5 _ 6 


5-12 M 5 R 5_ 12 
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{a, 


6-5 M,R, , 
6 6-5 


6-6 


M 6 R 6-5 


{A. 


7 - 2 - M 7 R 7 _ 2 ’ 


7-7 


M 7 R 7-2 


‘8-4 M 8 R 8 _ 4 


8 - 10 M 8 R 8-10 


8 ‘ 5 M 8 R 8-5 


A = - 1 ( 1 + l + _ JL _) 
8-5 M g Rg_ 4 R g _ 5 Rg_ 10 


{A 10-4 " M,„rL , » A 10-8 M, rt R,„ „ ’ A 10-10 M,«'R 


10 10-4 


1 + 


‘ 10 * 10-8 


10 10-4 10-8 


D 1-1 -■ ^ {A i € l t(1 - p i> h !-i + (1 - p i> T 10 G H l-10 H 10-10 H ?0-l 


- 1 ] 


+ aJ «°[(1 - pj) - 1] + a l < J[(l - p*) H*_ 1 - 1]} 


e e r 


ev „e 


^ o r 4 c c x c TT c , A e e* e „e 

D l-2 ' Mj A 1 ^ " p l* *2 *1-2 + A 1 (1 “ P 1 J *2 H l-2 


d i -3 • Iq (1 


e \ 6 TT 6 

P l ) £ 3 H l-3 


°l-5 " M l * A 1 (1 p 1^ e 5 T 10 G H l-10 H 10-5* 


1-4 


O . C * i C \ C ..C 

— A x (1 - p ^ f 4 H x _ 4 


1-6 


O A a /i a r a b 

M x A 1 (1 “ P 1^ 6 T 10 G H l-10 H 10-6 
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°l-8 " A 1 (1 " P 1 J * (e 8 H l-8 + € 9 H l-9 5 + T 10 G H l-10 H 10-10 
^ e 8 H 10-8 + 6 9 H 10-9^ 


D l-10 " A 1 6 10 H l-10 [1 + T 10 G H 10-10 (1 + T 10 H 10-10 )] 


B i " 5 ^ A i (1 ” P i ) T 10 T n G H 1-10 H io-n 


Vl “^ A 2 ^- p 2 > ^ H 2 - l +A 2 ^- p 2 > £ 1 H 2 - l ] 


"2-2 M„ 2 


S -^ ta - pJ ) h 2_ 2 - 1 ] + ta - 


il d 


p 2 } H 2-2 ' 1] 


+ A* €*[(1 - p®) H®_ 2 - 1]} 


°2-3 = U 2 (1 * p 2 } 6 3 H 2-3 + A 2 (1 ‘ p 2> *3 H 2-3 } 


6 \ © it ® 


2-4 


- <4 (i - p c 2 > * c 4 v 4 ) 


M 2 ‘ 2 


„ a .d d. d u d 

} 2-7 ' M 2 A 2 (1 P 2 } € 7 H 2-7 


3-1 


° A 6 


es e „e 


-3 a - »,) «x -3.3 


>3-2 ' t 3 tA 3 (1 - <4> '2 H 3-2 + *3 (1 


e \ e ti e 1 

P 3 ) € 2 H 3-2^ 


° (4 6^[(1 - P 3 ) H ^_ 3 - 1] + A® €®[<1 


3-3 M 3 3 3 


3 3' 


p 3> *3-3 - 1] 


+ A 3 £ 3 [ (1 - P3 ) H 3 _ 3 - 1 ]} 
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> 3-7 - 5 ^* 3 '? (1 ' p 3 > > 3-7 + *3 a - p*) Hj_ ? } 


>4-1 ■ 57 A 4 O - »V ‘i >4-1 
4 


TV 0 ^ A C / -| Cv C -rjC 

°4-2 " M. A 4 (1 P 4 € 2 H 4-2 
4 


4-4 


Cf A C C r /- C\ yC 1 «| 

M~ A 4 € 4 1 ” P 4^ H 4-4 “ 

4 


5-1 


M 5 ^ A 5 ^ P 5^ C 1 T 10 G H 5-10 H 10-l^ 


°5-5 " M 5 ^ A 5^ € 5^ ^ " p 5 ) H 5-5 “ ^ + ^ p 5^ T 10 G H 5-10 H 10-5 H 10-10^ 


_ a « b /- bv b r «»b § 2 p «.b „b i 

°5-6 = A 5 1 P 5^ e 5 H 5-6 + T 10 G H 5-10 H 10-6 H lO-lO^ 


°5-8 = M 5 A 5 (1 " P 5 ) T 10 G H 5-10 (< 8 H 10-8 + *9 


°5-10 " M 5 " p 5^ *10 H 5-10^ + T 10 G H 10-10 (1 + T 10 H 10-10^ 

B 5 * A 5 (1 P 5> T H tT 10 G H 5-10 H !o- 10 H 10-ll + W 


Vl^ A 6 f l »“pS> ^10 G H 5-10 H 10-l 


n O . b /i bv b r.-b . 2 n tjb »jb i 

°6-5 _ m7 A 6 ^ “ P 6^ 6 6 H 6-5 + T 10 G H 6-10 H 10-5 H 10-10-* 

0 


°6-6 “ M A 6 f 6 ^ (1 ~ p 6^ H 6-6 “ ^ + (1 P 6 ) T 10 G H 10-5 H 10-10* 

D 
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°6-8 = M, A 6 (1 ‘ P 6 ) T 10 G H 6-10 8 H 10-8 + 4 9 H 10-9 ) 

6 


O *b 


bs b „b 


D 6-10 " m 7 A 6 (1 “ P 6 } 4 10 H 6 - 10 fl + T 10 G H 10-10 (1 “ T 10 H 10 - 10 );1 
6 


B 6 - 57 A 6 (1 - p 6> Wll + T w G "6-10 H 10-10 "lO-lO 1 
6 


a ,d /. d N d TT d 
°7-2 " A 7 a * P 7 } c 2 H 7-2 


°7-3 3 {A 7 (1 " P ^ 4 3 H 7-3 + A 7 (1 " P 7 ) 4 3 H 7-3 3 


D 


7-7 " ^ {A 7 4 7 [(1 - p 7> H 7-7 ‘ 1] + A 7 4 7 [(1 ’ p 7> H 7-7 “ l3> 


f f, 


fv „£ 


o r 8 u a u. r, w a vfi H 3 ] 

u i ri i J 


8-1 M ‘a 

a A (8+9) 


<-r— a - p i^ 4 I h 8-i + 4 1 T io G Vio H io-io ^io-i 


+ _ (1 - P 9 ) + T 10 G H 9-10 H 10-10 H 10-l 33 A 


(8+9) 


(8+9) 


°8-5 = t {A 8 (1 * p 8> 4 5 T 10 G H 8-10 *lo-5 + A 9 (1 ‘ p 9 } 4 5 T 10 G *9-10 H ^0-5 } 
8 


°8-6 “ M - {A 8 (1 ' P 8 } 4 6 T 10 G H 8-10 H l0-6 + A g (1 * p 9 } 4 6 T 10 H 9-10 H 10-6 } 


a N b 


D 3-8 - k {A 8 [ '8 <a ' P 8 ) H 8-8 * 1J + (1 " p 8> ho C H 8-10 “w-lO (H 10-8 


‘l + “? 0-9 '? + (l - P? 4 9 H 8 “ 9 ^ + A ; i(1 - » 9 > ‘9 " 9-8 + *5 

(( 1 - P 9 ) H 9 _ 9 - 1 ) + (X - p g ) T 10 G H g _ 10 H ip_lq^ H io -8 * 8 + H 10-9 * 9 ^ 
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5 8-10 “ M fl * A 8 (1 " p 8* ‘ 10 H 8-10 11 + T 10 G H 10-10 (1 + T 10 H 10-10 )] 
o 

+ A 9 (1 P 9^ 4 10 H 9-10^ + T 10 G H 10-10 (1 + T 10 H 10-10*^ 


B 8 " M fl * A 8 (1 “ P 8^ T 10 T 11 G H 8-10 H 10-ll + A 9 (1 " p 9^ T 10 T 11 G 9-10 H 10-ll^ 

o 


>10-1 ■ 4 Ko •; Vl 1 * 1 - + T io G((1 - "Io> H 10- 10 

+ ^o e i T io G[(1 " p io ) H io-io “ 1] H io-i* 


h io - iq 1 


°10-5 " M 10 * A 10 ‘5 T 10 G[(1 " p 10* H 10-10 " 1 H 10-5 

+ A 10 € 5 H 10-5^ (1 ■ p 10* + T 10 G (1 ■ P 10 J H 10-10 ■ H 10-10^ 


° 10-6 ” M ,« tA 10 4 6 T 10 G[(1 p 10 ) H 10-10 1 H 10-6 


10 

+ A 10 ( 6 ^O-S^ 1 ~ P 10^ + T 10 G ^ 1 “ P 1 0^ H 10-l 


' 10 ' 10-10 


1 H 10 - 10 jJ 


■*10-8 " M 1q ^ A 10 (t 8 H 10-8 + 4 9 ^0-9* [ (1 ’ p 10* + T 10 G((1 p 10^ H 10-10 -1 ^ 

^O-IO 1 + A 10 T 10 G[<1 ~ P 10 ) H 10-10 " 1](< 8 H 10-8 + ‘ 9 H 10-9^ 


D io-io m 10 * A 10 £ 10 [(1 p 10 ) H 10-10 " 11 [1 + T 10 G H 10-10 (1 + T 10 H 10-10 )] 

+ A io 4 10 f(1 ' p 10 ) H 10-10 ” 11 11 + T io G H io-io (1 + T 10 H 10-10 )J * 
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1 r 




Bin * u {A, n n ^[(1 — P^q) ®]_0— 10 "*" J “10— 11 


10 M 10 1 10 10 11 


- 1 ] h: 


+ A 10 T 10 T 11 G[(1 ■ P 10 ) H 10-10 ■ 11 H 10-10 H 10-ll 


b ,, b <. b „b i 

+ A 1q (1 - p iQ ) T u H 10 _ 1;L } 


B 10 = M^. T 1 H 10-ll {A 10 T 10 Gt(1 " P 10 ) H 10-10 “ 1] 

+ A 10 (T 10 G H 10-10 “ p 10 ^ H 10-10 1 -* + 1 P 10 )} 
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APPENDIX D 


ERBE SCANNER THERMAL MODEL EQUATIONS 


PKEUED1NU PAGE BLANK NOT FILMED 



SCANNER MODEL 


The scanner model was developed by the same methods which were used 
for the non-scanner model development. The Equations (1) through (7) in 
the Appendix C are used to generate the constants of the thermal model 
Equation (24) in the Section 2.1.3. 

To define the coefficients in the themal model, the enclosure coupling 
technique was used as another alternative approach. The following pages 
contain those coefficients defined by using the enclosure coupling method. 

The attached figure shows the scanner model which has node description 
and enclosure descriptions. The enclosure has the number which is circled 
around. The numbers describing enclosure surfaces are underlined. In the 
illustration of the coefficients, the word that explains the procedure of 
derivation was omitted because they were already shown in Appendix C. 
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Heat Sink & Detector Assy 
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SCANNER MODEL 





" 1 * 1-2 


1,2 M 1 R 1 _ 2 


l l,3 ’ *1,4 " *1,5 ‘ *1,6 ’ *1,7 ' *1,8 = *1,9 ’ *1,10 * *1,11 = *1,: 


A = A =0 
1,13 1,14 


2,! ' 

, = 1 — 

2,3 M 2 R 2-3 


A 2 » 2 = ( M 2 Hr 2-1 + R 2-3 + R 2-4 + R 2-7 + R 2-8 + R 2-13 } 


, A, 


2,4 - M 2 R 2 _ 4 


A 2,5 " A 2,6 ° 


2,7 M 2 R 2 _ 7 


2,8 M_R, 


2 2-8 


2,13 M R 


2 2-13 


» A 2,9 A 2,10 A 2,ll A 2,12 A 2,14 ° 


A =0 
i 3,1 


3,4 3 , 


*4,1 * ° 


(— ) ( — + — + — 

M 4 R 4-2 R 4-6 R 4-8 


‘ 4 > 6 " 4 * 4-6 ’ A4 ’ 7 ' ° 

^4,9 = A 4,10 = A 4,ll = A 4,12 ~ A 4,13 “ \,14 


"3*3-2 

’ A 3,3 

M 3 R 3-2 

= A 

= A ^ = 

: A„ = A, 

3,9 

3,10 

3,11 : 

1 

A, „ 

= 0 

'4*4-2 

’ 4,3 


) 

. A ’ 

= 0 


’ 4,5 

1 


00 

*5 

e* 

‘ "4*4-8 


,A, „ - A. - A, A..* A. _ A, 0 
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A 5,l = A 5,2 " A 5,3 _ A 5,4 ° 




5 , 5 M 5 R 5 _ 6 


'5,6 M 5 R 5 _ 6 


A 5,7 " A 5,8 “ A 5,9 " A 5,10 A 5,ll " A 5,12 " A 5,13 A 5,14 = ° 


^ A 6,l " A 6,2 " A 6,3 ° 


6,4 


6,6 


6,7 


M 6 R 6-4 


6 ’ 5 M 6 R 6-5 


' + R 6-5 + R 6-8 ) 


6 ' 8 « 6 R 6-8 


V A 6,9 ” A 6,10 A 6,ll A 6,12 " A 6,13 " A 6,14 ° 


A 7,l * ° 


’ 7,2 M ? R 7 _ 2 


A = A_ - = A _ » A , = 0 
7,3 7,4 7,5 7,6 


*7.7--<s?<cT + ir-> 


*7 7-2 ‘7-8 


7,8 M 7 R 7-8 


V A 7,9 “ A 7,10 = A 7,ll A 7,12 “ A 7,13 A 7,14 " ° 
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f A 


8,1 


0 , A 

0 , A 


8,2 ' M 8 R 8 _ 2 • A 8,3 ‘ ° • V ’ M 8 R 8 _ 4 


8,5 ’ “8,6 M 8 R 8 . 6 


, A, 


8 , 7 M 8 R 8 _ 7 


L= - (^-)(^ — + + ~~ + i~-) 

8,8 M g R g _ 2 R g_4 R 8-6 R 8-7 K 8-10 K 8-ll 


a 8 , 9 - 0 , a. 


. K 


8,10 M 8 R 8-10 8,11 M 8 R 8-11 


^ A 8,12 A 8,13 = A 8,14 ~ ° 


A 9,l A 9,2 ” A 9,3 " A 9,4 A 9,5 ~ A 9,6 A 9,7 = A 9,8 " ° 


<A 


9,9 M 9 R 9-12 


» A 9,10 ° ’ A 9,ll » A 9,12 M 9 R 9 _ 12 


A 9,13 “ A 9,14 ° 


10,1 ” A 10,2 = A 10,3 A 10,4 A 10,5 ' A 10,6 A 10,7 


’ A 10-9 ° ’ 


10,8 M io R io-8 


^o-io = - (rr-x^— + 

M 10 i0-8 K 10-ll k 10-14 


’ A 10-12 ° ’ A 10,13 


v 10-11 “ M 10 R 10-11 


= 0 , A 


10 ’ 14 M 10 R 10-14 
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A ll,l ■ A ll,2 ' A ll,3 A ll,4 A ll,5 = A ll,6 " A ll,7 = ° 


n ’ 8 * 11 * 11-8 ’ 


^11*9 ’ All » 10 M ii R n-iO 


11,12 


, A 11,13 = A ll,14 ” ° * A ll,ll ' " < M 11 )( R 11 _g + R 11 _ 10 ) 


A —A =A = A s A = A = A * A = 0 

*12,1 A 12,2 A 12,3 *12,4 *12,5 *12,6 *12,7 12,8 


L 12,9 " M 12 R 12 _ 9 ’ *12, 10 = A 12,ll 


= 0 , A. 


12,12 M 12 R i 2-9 


A 12,13 " ° ’ A 12,14 = ° 


13,1 = 

° ’ A 13,2 ~ M 

13*13-2 

• *13,3 = A 13,4 


A — A 

= A 

a A ^ = A ^ 

13,7 = 

13,8 13,9 

13,10 

13,11 13,12 


1 

A 

= 0 

13,13 

M 13 R 13-2 

13,14 



A = A 

= A 

A = A , , = 

14-1 = 

14,2 14,3 

14,4 

14,5 14,6 


1 

A = 

0 , A = 0 

14,10 

M 14 R 14-10 

14,11 

’ 14,12 


1 



14ml4 

M 14 R 14-10 




14,13 
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For the total channel, the filter glasses are removed. Accordingly, 

the (— - — ) related with node number of filters have to set zero. 

R . . 
i-J 


B n = [0, 0, 0, 0, 0, 0, 0, i-, 0, 0, 0, 0, 0, 0] 

x 8 


Scanner : 


Let 


(H kji } A " ( ^k-j P J-i*A 


where subscript A denotes the enclosure A. 
Then let 


G 1 = t 1 T 13 ^ H 13jl3^D ^ H 13ml3^F^ 


S 1 " ^ H 17ml7^F + T 13 G 1 ^nd^F (H 13jl3 ) D (H 13ml7 ) F 


°2 ' 11 - h7 S l (H 1 7jl7 ) c rl 


S 2 " + T 17 S 1 G 2 (H 17jl6 ) C (H 16jl7 ) C 


g 3 - u - t 2 16 s 2 (H 16 ml 6 > B r 1 


S 3 " ^ H 14ml4^B + X 16 S 2 G 3 (H 14ml6 ) B (H 16ml4 ) B 


G 4 X 14 S 3 (H 14jl4 ) A ] 
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X 1 - F 1,A + F l,2 n 2 + P 1,A + F l,5 n 5 + F lA + F l,8 n 8 + F l,9 n 9 


+ F l,12 !i 13 


where 


F l,l T 13 f l G 1 ^ H 13jl^D 


F l,2 “ € 2 G 2 t(H 17m2 ) F + T 13 (H 13j2 ) D + T 13 M^m^F 


^ H 13m2^F^ 


F, . = e. G {(H ,)_ + T.,, G (H (H (H., .)_} 

1,4 4 2 17m4 F 13 1 17ml3 F 13ml3 D 13m4 F 


F l,5 = T 17 ( 5 S 1 G 2 (H 17j5 ) C 


F l,6 = T 17 € 6 S 1 G 2 (H 17j6 ) C 


F l,8 “ T 17 ( 8 S 1 G 2 (H 17j8 ) C 


F l,9 T 17 C 9 S 1 G 2 (H 17j9 ) C 


F l,13 " ‘l3 G 2 t(H 17ml3 )F + T 13 G l (H 17ml3 ) F (H 13jl3*D + T 13 G 1 (H 17ml3 ) F 
(H 13jl3*D ^ H 13ml3^F^ 


X 2 ' F 2,A + F 2,2 Q 2 + F 2,A^4 + F 2,5" Q 5 + F 2,6^6 + F 2,8 n 8 + F 2,9 fi 9 + F 2,lAl 


+ - F 2, 12 fl 12 + F 2,13 P ‘13 + F 2,1^14 
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where 


F 2,l “ T 13 T 17 G 1 °2 G 3^ H 16jl7^C (H 17ml3 ) F (H 13jl ) D 


F 2,2 “ X 17 f 2 G 2 G 3 (H 16jl7 ) C * (H 17m2*F + X 13 G l fH 17iiil3*F (H 13j2 ) D + T 13 G 1 


^ H 17ml3^F * H 13jl3^D 


F 2,4 = X 17 £ 4 G 2 G 3 (H 16jl7 ) C ^ (H 17m4 ) F + T 13 G l (H 17ml3 ) F (H 13m4 ) F^ 

F 2,5 “ £ 5 G 3 ^ (H 16j5 > C + T 17 S 1 G 2 (H 16jl7 ) C (H 17j5 ) C } 

F 2,6 “ 6 6 G 3 ^^6^ + T 17 S 1 G 2 (H 16jl7 ) C (H 17j6 ) C J 

F 2,8 = *8 G 3 *^638*0 + X 17 S 1 G 2 (H 16jl7 ) C (H 17J8*C* 

F 2,9 = 4 9 G 3 ^ (H 16j9 ) C + T 16 V^n^B + T 17 S 1 G 2 (H 16jl7 ) C (H 17j9 } C* 


F 2,ll " X 16 'll S 2 S^iemll^B 


F 2,12 “ T 16 € 12 S 2 G 3 (H 16ml2 ) B 


F 2,13 " T 17 f 13 G 2 (H 17tnl3 ) F {l + X 13 + X 13 G 1 

^ H 13jl3^D ^ H 13ml3^F* 


F 2,14 X 16 6 14 S 2 G 3 (H 16ml4 ) B 
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E 14 ■ °1,A + U l,2 n 2 + “l,A + U l,5 n 5 + "l A + U 1 A + “l A + U l,l<Ao 
+ "l.lAl + U l,12 n i2 + U l,13 n X3 + U I,14 a i4 + U 1,15 E 


U - S T . G/ (H- m F 

1,1 16 4 14ml6 B 2,1 


U l,2 " X 16 F 2,2 


U . * T , G. (H, . . F- , 
1,4 16 4 14ml6 B 2,4 


U l,5 T 16 F 2,5 


U, * - T., G. (H,. .,)„ F- - 
1,6 16 4 14ml 6 B 2,6 


U l,8 " T 16 G 4* H 14ml6*B F 2,8 


U l,9 ■ G 4 {T 16 (1 W6>B F 2,9 + VWb 1 


U 1,10 " T 10 S 3 G 4 < 10 (H 14jl0 ) A 


U. - , • G. {« (H . n ) + T.-.-CH... iJ 

1,11 4 11 14mll B 16 14ml6 B 2,11 


U l,12 " G 4 t€ 12 (H 14ml2 ) B + T 14 *12 S 3 (H 14jl2 ) A + T 16 (H 14ml6 ) B *2,12* 


U l,13 " T 16 G 4 (H 14ml6 ) B F 2,13 


U l,14 = G 4 tc l4 (H 14ml4 ) B + X 14 *14 S 3 (H 14jl4 ) A + T l6 (H 14ml6 ) B F 2,14* 


U l,15 ” X 15 °4 ^14jl5*A 
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U 2,l = T 14 T 16 G 4 (H 14jl4 } A F 2,l 

U 2,2 = T 14 T 16 G 4 (H 14hl4 ) A (H 14ml6 ) B *2,2 

U 2,4 = T 14 T 16 ^ H 14ml6^B F 2,4 

U 2,5 = T 14 X 16 °4 ^ H 14 jl4^ A ^14*16^ *2,5 

U 2,6 = T 14 X 16 G 4 (H 14jl4 ) A F 2,6 

U 2,8 = T 14 T 16 V^JiVa (H 14ml6 ) B ? 2,8 

U 2,9 = V^Wa {T 14 T 16 (H 14ml6 ) B F 2,9 + T 14 t 9 (H 14m9 ) B } 

U 2,10 = * 10^ H 14jl0^A {1 + T 14 S 3 °4 (H 14jl4 ) A } 

U 2,ll " G 4 (H 14jl4 ) A {T 14 T 16 (H 14ml6 ) B *2,11 + T 14 € ll (H 14mll ) B } 

U 2,12 = T 14 T 16 V^j^A (H 14ml6 ) B *2,12 + € 12 (H 14jl2 ) A tl + 

+ X 14 ( 12 G 4 (H 14jl4 ) A (H 14ml2 ) B 

U 2,13 = T 14 T 16 G 4 (H 14jl4 ) A ^lAm^B F 2,13 


U 2,14 " {t 14 *14 G 4 (H 14ml4 ) B + T 14 X 16 G 4 (H 14ml6*B F 2,14 + *14 


t 1+X 14 S 3 VVjU^ 


U 2,15 " T 15 (H 14jl5 ) A ^ + T 14 G 4 (H 14jl4 ) A } 


E 16 " U 3,l^l + U 3,2^2 + U 3,4 a 4 + U 3,5 fi 5 + U 3,6^6 + U 3,8^8 + U 3,9S + °3, 10 fl 10 
+ U 3, 11^11 + U 3, 12^12 + U 3, 13^13 + U 3,14 fi 14 + U 3,15 E 


U 3,l “ T 14 T 16 S 2 G 3 ^61114^ B U 2,l + X 13 X 17 G 1 G 2 °3 *l (H 16jl7^C ^ H 17ml3*F 

U 3,2 * X 14 X 16 S 2 G 3^ H 16nl4 ) B U 2,2 + X 13 X 17 *2 G 1 G 2 G 3 (H 16jl7 i C (H 17ml3 ) F 
(H 13jl3 J D + X 17 6 2 G 2 G 3 (H 16jl7 ) C (H 17m2 ) F + X 13 X 17 *2 G 1 °2 °3 
(H 16jl7*C (H 17ml3 ) F (H 13jl3*D * (H 13m2 ) F 


U 3,4 “ X 14 X 16 S 2 G 3^ H 16ml4^B U 2,4 + T 17 *4 °2 G 3^ H 16.il7^C ^ H 17m4^F 
+ X 13 X 17 *4 G 1 G 2 G 3 (H 16jl7 ) C (H 1 71113^ ^lSjU^D (H 13m4 ) F 


U 3,5 * X 14 X 16 S 2 G 3 (H lftnl4 ) B U 2,5 + °3 € 5 (H 16j5*C + X 17 *5 S 1 ®2 ^*16^7*0 
(H 17j5 ) C 

U 3,6 = X 14 X 16 G 3 (H 16ml4 ) B U 2,6 + *6 + X 17 f 6 S 1 °2 °3 (H 16jl7 ) C 


(H 17j6 ) C 
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U 3,8 " X 14 T 16 S 2 G 3 (H 16ml4 ) B U 2,8 + *8 G 3 (H 1638*C + X 17 *8 S 1 °2 G 3 (H 16jl7 ) C 


^ H 17j8^C 


U 3,9 = T 14 X 16 S 2 G 3 (H l6ml4 ) B U 2,9 + *9 G 3 (H 1639*C + X 17 *8 S 1 °2 


(H 17j9 ) C + X 16 * 9 S 2 G 3 (H 17j9 ) C 


U 3,10 T 14 T 10 S 2 G 3 (H 16ml4 ) B U 2,10 


U 3,ll = T 14 T 16 S 2 G 3 (B 16iii14 ) B U 2,ll + T 16 *11 S 2 G 3 (H 17jll ) C 
U 3,12 = T 14 T 16 S 2 G 3 (H 16ml4 ) B U 2,12 + X 16 *12 S 2 G 3 (H 17jl2 ) C 

U 3,13 = X 14 X 16 S 2 G 3^ H 16ml4 ) B U 2,13 + X 17 *13 G 3 ^*16317*0 (H 1 7ml3 ) F 

+ T 13 X 17 *13 G 1 G 2 G 3 (H l6jl7 5 C (H 13jl3*D (H 17ml3 ) F [1 + X 13 (H 13ml3 ) F ] 


U 3,14 = X 14 X 16 S 2 G 3 (H 16nl4 ) B U 2,14 + X 16 *14 S 2 G 3 ^17314*0 


U 3,15 ’ X 14 X 16 S 2 G 3 (H 16ml4 ) B U 2,15 


E 16 " \l*l + U 4,2“ Q 2 + U 4,4°4 + U 4,5 B 5 + U 4,6°6 + + + \l<PlO 

+ U 4, 11 Q 11 + U 4,12 n i2 + U 4, 13^13 + U 4, 14^14 + U 4,15 E 

U 4,l = X l4 (H l6ml4 ) B [1 + X 16 S 2 G 3 (H 16ml6V U 2,l + X 16 X 17 G 3 (H 16ml6 ) B 


^6317^ F l,l 
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U 4,2 ' T 14 (H 16ml4 ) B [1 + X 16 S 2 G 3 (H 16ml6V U 2,2 + X 16 
F l,2 

U 4,4 = X 14 (H 16ml4 ) B [1 + X l6 S 2 G 3 (H 16ml6 ) B ] U 2,4 + X 16 
^jl^C F l,4 

U 4,5 = X l4 (H 16ml4 ) B [1 + X 16 S 2 G 3 (H 16ml6V U 2,5 + X 16 
F l,5 + X 16 f 5 G 3 (H l6ml6 ) B ( ®L6J5*C 

U 4,8 = X 14 (H 16ml4 ) B tl + X 16 S 2 G 3 (H 16ml6V U 2,8 + X 16 
F l,8 + X 16 € 8 G 3 (H 16ml6 ) B ^jB^C 

U 4,9 = X 14 (H 16ml4 ) B tl + X 16 S 2 °3 (H 16ml6V U 2,9 + X 16 
^ H 16jl7^C F l,9 + + X 16 *9 G 3 (H l6ml6 ) 

+ X 16 € 9 S 2 G 3 (H 16ml6 ) B (H 16m9*B 

U 4,10 = X 14 (H l6ml4 ) B [1 + X l6 S 2 G 3 (H 16ml6 ) B 1 U 2,10 

U 4,ll = T 14 (H 16ml4 ) B 11 + X 16 S 2 G 3 (H 16ml6 ) B 1 U 2,ll + 6 
+ X 16 £ 11 S 2 G 3 (H l6ml6 ) B ^Bmll^B 

U 4,12 = T 14 (H 16ml4 ) B [1 + X 16 S 2 G 3 (H 16ml6 ) B ] U 2,12 + € 
+ X 16 € 12 S 2 S^iend^B (H 16ml2^B 


X 17 

X 17 G 3 (H l6ml6 ) B 

X 17 G 3 (H l6ml6^B 

X 17 G 3 (H 16ml6 ) B 

X 17 G 3 (H l6ml6^B 
B (H 16j9 ) C 

ll* H l6mll*B 

12 (H 16ml2 ) B 



"4,13 = T 14 (H 16ml4 ) B fl + T 16 S 2 G 3 (H 16ml6 ) B 1 U 2,13 + T 16 X 17 C 3 (H 16ml6 ) B 
(H l6jl7 ) C F l,13 

U 4,14 = T l4 (H 16ml4 ) B tl + X 16 S 2 G 3 (H 16ml6 ) B 1 U 2,14 + X 16 *14 S 2 °3 <H 16ml6 ) B 
^ H 16ml4^B + < 14 (I, 16ml4 ) B 


U 4,15 T 14^ H 16ml4^B^ 1 + T 16 S 2 G 3^ H 16ml6^ U 2,15 


g 17 ’ U 5,l“l + ”5,2 n 2 + U 5,4 n 4 + U 5,5 n 5 + "sA + “sA + U 5,9 n 9 + U 5,10 n i0 
+ “s.lAl U 5,12 n i2 + U 5,13 n i3 + U 5,14 n i4 + ° 5,15 E 


u 


S 1 - ( 1 T !3 G i G 2^ H 17ml3 ) F ^ H 13jl^D + T 16 T 17 S 1 G 2 (H 17jl6^C U 4,l 


U 5,2 ~ C 2 G 2 (H 17m2 ) F + *2 T 13 G 1 G 2 (H 17ml3 ) F [ (H 13j2 ) D + (H 13m2 ) F 

T 13 ] + T 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,2 

U 5,4 = € 4 G 2 (H 17m4 ) F + *4 X 13 G 1 G 2 (H 17ml3 ) F ^UjlS^D (H 13m4 ) F + T 16 X 17 
S 1 G 2 < ' K 17jl6 ) C U 4,4 

U 5,5 = ‘5 X 17 S 1 G 2 (H 17j3 ) C + X 16 X 17 S 1 G 2 (H 17jl6 } C U 4,5 

U 5,6 = ‘6 X 17 S 1 G 2 (H 17j6 } C + X 16 X 17 S 1 G 2 (K 17jl6 ) C U 4,6 

U 5,8 = ( 8 X 17 S 1 G 2 (H 17j8 } C + X l6 X 17 S 1 G 2 (H 17jl6 ) C U 4,8 
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U 5,9 = € 9 T 17 S 1 G 2 (H 17j9 ) C + T 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,9 


U 5,10 = T 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,10 


U 5,ll X 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,ll 


U 5,12 ~ T 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,12 


i,13 " € 13 G 2 (H 17ml3 ) F + 'l3 T 13 G 1 G 2 (H 17ml3 ) F (H 13jl3 ) D [1 + T 13 (H 13ml3 ) F ] 


+ X 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,13 


U 5,14 " T 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,14 


U 5,15 " T 16 T 17 S 1 G 2 (H 17jl6 ) C U 4,15 


E 17 " U 6,A + U 6,2^2 + U 6,4 fi 4 + U 6,5S + U 6 y 6^6 + U 6,8 fi 8 + U 6,9 P "9 + V 10^10 
+ U 6,lAl + U 6, 12^12 + U 6,13 n i3 + V 14^14 + U 6,15 E 

U 6,l = f l T 13 T 17 G 1 G 2 (H 17jl7 ) C (H 17ml3 ) F + T 16 G 2 (H 17jl6 ) C U 4,l 

U 6,2 = f 2 T 17 G 2 (H 17jl7 ) C (H 17m2 } F + £ 2 T 13 T 17 G 1 G 2 (H 17jl7 ) C (H 17ml3 ) F 
^ ( H 13 j 2>D + T 13^ H 13jl3^D * H 13m2^ + X 16 G 2^ H 17jl6^C U 4,2 

U 6,4 = € 4 T 17 G 2 (H 17jl7 ) C (H 17m4 ( F + % X 13 X 17 °1 G 2 (H 17jl7 ) C (H 17ml3 ) F 


(H 13jl3 J D (H 13m4 ) F + X 16 G 2 (H 17jl6 ) C °4,4 
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U 6,5 " *5 G 2 (H 17j5 ) C + T 16 G 2 (H 17jl6 ) C U 4,5 


U 6,6 ' *6 G 2 (H 17j6 } C + T 16 G 2 (H 17jl6 ) C U 4,6 


U 6,8 f 8 G 2 (H 17j8 ) C + T 16 G 2 (H 17jl6 ) C U 4,8 


U 6,9 f 9 ( V H 17j9^C + T 16 G 2 ( ‘ H 17jl6^C U 4,9 


U 6 , 10 T 16 G 2 (H 17jl6 ) C U 4,10 


U 6,ll ” T 16 G 2 (H 17jl6 ) C U 4,U 


U 6,12 X 16 G 2 (H 17jl6 ) C U 4,12 


U 6,13 = f 13 T 13 T 17 G 1 G 2 (H 17jl7 ) C (H 17ml3 ) F + T 13 (H 13ml3 ) F 1 

+ £ 13 T 17 G 2 (H 17jl7 } C (H 17nl3 ) F + T 16 G 2 (H 17jl6 ) C U 4,13 


U 6,14 T 16 G 2 (H 17jl6 ) C U 4,14 


U 6,15 T 16 G 2 (H 17jl6 ) C U 4,15 


E 13 = U 7, A + U 7,2“ Q 2 + U 7,4^4 + YsS + U 7,6^6 + U 7 ^8 + U 7,A + U 7,10 n 


10 


+ U 7,lAl + U 7,12%2 + U 7, 13 Y + U 7, 14^14 + U 7,15 E 


U 7,l “ ( i G l (H 13jl ) D + T 13 T 17 G l < ‘ H 13jl3 ) D < ' H 13ml7 ) F U 6,l 
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U 7,2 ( 2 G l (H 13j2 ) D + C 2 T 13 G l (H 13jl3 ) D (K 13m2 ) F + T 13 T 17 G l (K 13jl3 ) D 


(H 13ml7 ) F U 6,2 


U 7,4 * 4 T 13 G l (H 13jl3 ) D * H 13m4*F + X 13 X 17 G l ( ' H 13jl3 ) D (H 13ml7 ) F U 6,4 


U 7,5 X 13 X 17 G l (H 13jl3 ) D (H 13ml7 ) F U 6,5 


U 7,6 X 13 X 17 G l* H 13jl3*D (H 13ml7 ) F U 6,6 


U 7,8 = X 13 X 17 G l (H 13jl3 ) D (H 13ml7 ) F U 6,8 


U 7,9 X 13 X 17 G 1 ^lSj^D ^SndJ^F U 6,9 


U 7,10 X 13 X 17 G l (H 13.il3 ) D (H 13ml7*F U 6,10 


U 7,ll X 13 X 17 G l (H 13jl3 ) D (H k3ml7 ) F U 6,ll 


U 7,12 X 13 X 17 G 1 ^ H 13jl3^D (H 13ml7 ) F U 6,12 


U 7,13 £ 13 + ( 13 X 13 G l* H 13jl3*D * H 13ml3*F + X 13 X 17 G 1 

(H 13ml7*F U 6,13 


U 7,14 X 13 X 17 G l (H 13jl3 ) D (H 13ml7 ) F U 6,14 


U_ = T _ T._ G.(H. (H 17 )_ U _ 
7,15 13 17 1 13jl3 D 13ml7 F 6,15 
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\ 3 = U 8,A + U 8,2 Q 2 + U 8,4 Q 4 + °8 A + U 8 A + U 8 A + U 8 A + YAo 
+ U 8,lAl + V 12^12 + U 8, 13°13 + U 8?14 14 + U 8,15 E 


U 8,l = e i T 13 + T 17 G l (H 13ml7 ) F U 6,l 

U 8,2 = f 2 X 13 G l lH 13inl3 ) F (H 13i2 ) D + € 2 SAAf + T 17 G l (H 13ml 7 } F U 6,2 

U 8,4 =£ 4 G l (H 13m4 ) F +T 17 V^Af U 6,4 


U 8,5 = X 17 G l (H 13ml7 ) F U 6,5 


U 8,6 " T 17 G l (H 13ml7 ) F U 6,6 


U 8,8 = X 17 G l (H 13tnl7 ) F U 6,8 


U 8,9 ' T 17 G l (H 13ml7 ) F U 6,9 


U 8,10 X 17 G l (H 13ml7 ) F U 6,10 


U 8,H " T 17 G l (H 13ml7 ) F U 6,ll 


U 8,12 X 17 G l (H 13ml7 ) F U 6,12 


U 8,13 ‘ 13 


T 13 G l (H 13ml3 ( F (H 13jl3 } D + ‘ 13 G l (H 13ml3 ) F + T 17 G l (H 13ml7 ) F U 6,13 


U 8,14 = T 17 G 1< H 13»17 ) F U 6,14 


U 8 , 15 T 17 G l (H 13ml7 ) F U 6,15 
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0A 1 

D l,l " 6 1 ^ljl^D “ € 1 + (1_p l ) T 13 (H ljl3 ) D U 8,1 ] S.L 

°1,2 = (- M^ )l(1_p l ) *2^2^ + (1-p l ) T 13 (H ljl3 ) D U 8,2 ] S,L 


°1,3 " ° 


D 


1,4 


1,5 


1,6 


( M ^) Id-Pi) T 13 ^ H 1 ji3^ D U 8,4^S,L 
aA 

(- M^ )[(1 “ p l ) T 13 (H ljl3 ) D U 8,5 ] S,L 
aA 

( ~M^ )t(1 ~ P l ) T 13 (H ljl3 ) D U 8,6 J S,L 


° 1,7 " ° 

°1,8 = (_ M^ )[(1 “ P 1 ) T 13 (H ljl3 ) D U 8,8 ] S,L 

D l,9 = (_ M^ )[(1 ‘ P 1 ) T 13 (H ljl3 ) D U 8,9 J S,L 
aA 

°1,10 ~ (_ M^ [ ^“Pi^ T 13 (H ljl3 ) D U 8,10 J S,L 
aA 

D l,ll = ( ^’ )t(1_P l ) T 13 (H ljl3 ) D U 8,11^S,L 
aA l 

D l,12 “ ( T 13 (H ljl3 ) D U 8,12^S,L 

aA l 

°1,13 = ( M 1 )[(1-p l ) T 13 (H ljl3 ) D U l,13 + (1 P l ) 4 13 (H ljl3 Vs,L 
aA 

D i i /, = <— rr") t d-pT > t..(h . 1? ) n u ] 

1,14 1 13 ljl3 D 8,14 S,L 
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°2,1 = {A 2D (1 - p 2 )[£ l (H 2 j l ) D + T 13 (H 2 j 13 ) D V 1 + A 2F (1 - p 2 )lT 13 

(H 2ml3 ) F U 7,l + T 17 (H 2ml7 ) F U 6,1 ]} S,L 


°2,2 * ^ A 2E t(1 ' P 2 ) € 2 (H 2j2^E ' £ 2 ] + £ 2 (H 2 j 2 ) “ £ 2 

+ (l-p 2 ) T 13 ( H 2 jl3 ) D U 8,2 ] + A 2F f(1_P 2 ) £ 2 (H 2m2 ) F ' £ 2 + (1_p 2 ) 

+ (1-P 2 ) t 13 *fH 2ml3 ) F U ? ^ 2 + (1-P 2 ) T 17 ( H 2 m i 7 )F U 6,2 ]} S,L 

^M^ A 2E^ 1-p 2^ £ 3 (H 2j3^E-*S,L 

^ A 2D (j - p 2 ) T 13 (H 2jl3 ) D U 8,4 + A 2F t(1 - p 2 ) £ 4 (H 2m4 ) F + (1 ‘ P 2 ) T 13 
(H 2ml3*F U 7,4 + (1_p 2 ) T 17 (H 2ml7 ) F U 6,4 ]} S,L 


2,3 


2,4 


°2,5 = ^ {A 2D (1 " P 2 ) T 13 (H 2jl3 ) D U 8,5 + A 2F (1 " P 2 ) [T 13 (H 2ml3 ) F U 7,5 + T 17 


(H 2ml7 ) F U 6,3 ]} S,L 


D 2,6 " <§"> {A 2D (1 - p 2 ) T 13 (H 2 j 13 ) D U 8,6 + A 2F (1 - p 2 )l [ t 13 (H 2«13>F U 7,6 + T 17 


(H 2ml7 ) F U 6,6 ]} S,L 


°2,7 " ( M 0 )[A 2E (1_P 2 ) £ 7 (H 2j7 ) E ] S,L 


2,8 


<*£> {A 2D (1 ' P 2 ) T 13 (H 2jl3 } D U 8,8 + A 2F (1_P 2 ) [t 13 (H 2«13 ) F U 7,8 + T 17 
(H 2ml7 ) F U 6,8^S,L. 
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°2,9 " ( M^ {A 2D (1_P 2 ) T 13 (H 2jl3 ) D °8,9 + A 2F (1_P 2 ) tT 13 (H 2ml3 ) F U 7,9 
+ T 17 (H 2ml7 ) F U 6,9^S,L 

°2,10 = ( M^ {A 2D (1_P 2 ) T 13 (H 2jl3 ) D U 8,10 + A 2F (1_P 2 ) tx i3 (H 2ml3 ) F U 7,10 
+ T 17 (H 2ml7 ) F U 6,10^S,L 

D 2,ll = tA 2D (1_P 2 ) X 13 (H 2jl3 ) D U 8,ll + A 2F (1_P 2 ) [x l3 (H 2ral3 ) F U 7,H 

+ X 17 (H 2ml7 ) F U 6,11 ]} S,L 


°2,12 “ ( M 2 } ^ A 2D (1 P 2^ X 13 (H 2jl3 ) D U 8,12 + A 2F (1 " P 2 ) [x i3 (H 2ml3^F U 7,12 
+ X 17 (H 2ml7 ) F U 6,12^S,L 

°2,13 = ^ A 2D (1_P 2 )t€ 13 (H 2jl3 ) D + X 13 (H 2jl3 ) D- U 8,13 ] + A 2F (1 " P 2 ) te i3 

(H 2ml3 ) F + X 13 (H 2ml3 ) F U 7,13 + 17 (H 2ml7 ) F U 6,13^S,L 

°2 ,14 = < ’ A 2D (1_P 2 ) X 13 (H 2jl3 ) D U 8,14 + A 2F (1 * P 2 ) ^ x 13^ H 2ml3^F U 7,14 

+ T 17 (H 2ml7 ) F U 6,14 ]} S,L 


°3,1 = ° 


°3,2 ( M 3 ) A 3 t(1 p 3 ) 6 2 (H 3j2 ) E^S,L 
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D 3,3 " ( M^ A 3 [(1 P 3 } *3 (H 3j3 ) E ( 3 ] S,L 


°3,4 = °3,5 = °3,6 " ° 


°3,7 " ( M 3 } A 3 [(1 P 3 } VSjyVs.L 

°3,8 = °3,9 = D 3,10 = °3,11 = °3,12 = °3,13 = D 3,14 = ° 


c a _ 

- (-JT-) {(l-p 4 )tT 13 ( H 4 m i 3 ) F U 7>1 + T 17 (.H 4ml7 ) F U 6>1 ]> S)L 

9 II 


UA, « 

D 4,2 " t(1 "' > 4 )t '2 (H 4m2 ) F + T 13 <H 4ml3 , F U 7,2 + T 17 (H 4iia7 ) F U 6,2 ll S,L 


D 4,3 ■ ° 


Q A. 

D 4 4 ‘ ( TT ) '<W , 4 <H 4m4 ) F ' <4 + ‘W T 13 (H 4ml3 ) F C 7,4 + (1_P 4 ) T 17 
’ 4 

(H 4ml7*F U 6,4^S,L 

cA 

°4 5 = ( "m 7 ) {(1_p 4 )[T 13 (H 4ml3 ) F U 7,5 + T 17 (H 4ml7 ) F U 6,5 ]} S,L 
* 4 

ok 

°4,6 ‘ ( -H7 )t(1 - p 4 )[T 13 (H 4 m 13>F U 7,6 + T 17 (H 4»17 ) F Vs'Vl 


D 4.7 - ° 


°A 

°4,8 = ( lT ){(1 ' P 4 )[T 13 (H 4ml3 ) F U 7,8 + T 17 (H 4,17 ) F U 6,8 ]} S,L 
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c A, 


°4,9 “ ( M 4 ^ (1 P 4 )tT 13 (H 4ml3 ) F U 7,9 + T 17 (H 4ml7 ) F U 6,9^S,L 


ok 

D 4,10 = (_ M^' >{(1 ' P 4 )tT 13 (H 4ml3 ) F U 7,10 + T 17 (H 4ml7 ) F U 6,10^S,L 

Ok 

D 4,ll = (1_P 4 ) * T 13* H 4na3*F U 7,ll + T 17 (H 4ml7 ) F U 6,11^S,L 

ok 

°4,12 = (1_p 4 ) [X 13 (H 4ml3? F U 7,12 + T 17 (H 4ml7 ) F U 6,12^S,L 

oA a 

°4,13 = + T 13 (H 4ml3 ) F U 7,13 + T 17 (H 4ml7*F U 6,13^ 1_P 4^S,L 

Ok 

°4,14 = ( 'll^ ){(1 ' p 4 )[T 13 (H 4ml3 ) F U 7,14 + T 17 (H 4ml7 ) F U 6,14 ]J S,L 


ok 

°5,1 = ( "Mj^ (1 “ P 5^ T 16 (H 5jl6^c U 4,l = T 17 (H 5jl7 ) c U 5,1^S,L 

Ok 

°5,2 " ( -^ ){(1 - p 5 )tT 16 (H 5jl6>c U 4,2 + T 17 (H 5J17 > c U 5,2 ]} S,L 


°5,3 ‘ ° 


5,4 


PA 5 

( V U/ / +T n ( H q . 17 ) U. ]} 

M,. 5 16 5jl6 c 4,4 17 5jl7 c 5,4 S,L 


D 


5,5 


aA 

(—j i){(l-p )[t. a (H_..,) U +t 17 (H ) U • ] + [(1-p ) € 

5 16 5jl6 c 4,5 17 5jl7 c 5,5 5 5 


(H,.,) - €,]}_ 

5j5 c 5 S,L 


a A r 


D , , = (— ^>{(1- p )[c,(H..,) + T 1A (H ) U + t 17 (H ) U fi ]} 

5,6 M_ 5 6 5j6 c 16 5jl6 c 4,6 17 5jl7 c 5,6 S,L 
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v ■ 0 


q A 

°5,8 = ( lJ' ){(1 ' p 5 )[e 8 (H 5j8 ) c + T 16 (H 5jl6 ) c U 4,8 + T 17 (H 5jl7 ) c U 5,8 ]> S,L 
aA 

D 5,9 ' (_ M i){(1 ' p 5 ,le 9 (H 5j9 ) c + T 16 (H 5J16 ) c U 4,9 + T 17 (H 5jl7 ) c D 5,9 ], S,L 


D 5,10 ’ ( -5; ){(1 - C '5 >[T 16 (H 5 j 16 ) c U 4,10 + T 17 (H 5J17>c U 5,10 1} S,I- 

aA . 

°5,11 = (_ M^ ){(1_P 5 )[T l6 (H 5jl6 ) c U 4,ll + T 17 (H 5jl7 ) c U 5,H } S,L 

aA 

D 5,12 = (- M^ ){(1_P 5 )[T 16 (H 5jl6 ) c U 4,12 + T 17 (H 5jl7 ) c U 5,12 ]) S,L 
aA 

D 5,13 = (_ ^ ){(1_p 5 )tT l6 (H 5jl6 ) c U 4,13 + T 17 (H 5J17 ) c U 5,13 ]) S,L 

aA,. .. 

°5,14 = ( 1i^^ (1 “ P 5 ) ^ T 16^ H 5jl6 ) c U 4,14 + T 17^ H 5jl7' ) c U 5,14 S,L 


aA 

°6 1 = ( !T ){(1 “ P 6 )[T 16 (H 6jl6 ) c U 4,l + T 17 (H 6jl7 } c U 5,l ]} S,L 
’ 6 

aA 

S 6 2 = (- M~ ){(1 ' p 6 )[T l6 (H 6jl6 ) c °4,2 + T 17 (H 6jl6 } c U 5,2 ]} S,L 

* A 


D 6,3 ’ ° 


aA, 


°6 4 = (_ M i){(1_P 6 )[T 16 (H 6jl6 ) c U 4,4 +T 17 (H 6jl7 ) c U 5,4 ]) S,L 

* A 
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0 A 

D 6,5 ‘ ( 1T ){(1 -‘ ) 6 )l, 5 (H 6J5 ) c +T 16 (H 6 jl 6 ) c U 4,5 + T 17 (H 6jl7>c U 5,S I} S,L 


OA^ 


°6 6 = ( “M^ ){(1 ' P 6 )[t 6 (H 6j6 ) c + T l6 (H 6jl6 ) c U 4,6 + T 17 (H 6jl7 ) c U 5,6 “ Vs.L 


V ' 0 


ok. 


°6,8 = ( ~M7 ){(1_p 6 )[t 8 (H 6j8 ) c + T l6 (H 6jl6 ) c U 4,8 + T 17 (H 6jl7 ) U 5,8 ]) S,L 


Ok 


D 6,9 = (- M7 ){(1 " P 6 )[€ 9 (H 6j9 ) c + T 16 (H 6jl6 ) c U 4,9 + T 17 (H 6jl7 ) c U 5,9 1} S,L 


Ok 

°6,10 = ( -M^H< 1 - p 6 >t T i6 (H 6hl6 ) c U 4,10 + T 17 (H 6jl7 } c U 5,10 ]} S,L 

ok 

D 6,ll = ( ~M^ ){(1_P 6 )tT 16 (H 6jl6 ) c U 4,ll + T 17 (H 6jl7 } c U 5,11 ]} S,L 

ok 

°6 ,12 = ( 1T ){(1 " p 6 )[T l6 (H 6jl6 ) c °4,12 + T 17 (H 6jl7 ) c U 5,12 ]} S,L 
6 


a A. 


°6,13 “ ( M, ){(1 p 6 )[T 16 (H 6jl6 ) c °4,13 + T 17 (H 6jl7 ) c U 5,13 ]} S,L 


aA, 


D 6,14 " ( M, ){(1 " p 6 )tT 16 (H 6jl6 ) c U 4,14 + T 17 (H 6jl7 ) c U 5,15 ]} S,L 

D 


“7,1 ' ° 


D, „ = ( 


7,2 M 


ok 

( 1T ) {(1_P 7 ) £ 2 (H 7j2 ) E } S,L 
7 E 
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OA 

D 7,3 " ( T5^ ) e (( 1 ' P 7 ) t 3 (H 7j3 > E J S,L 


“7,4 ’ “7,5 - ”7,6 ' ° 


CA 


°7,7 “ ( M^ ) „ {(1_P 7 ) ~ € 7 } S,L 

°7,8 = °7 ,9 = °7,10 = °7,11 = D 7,12 = °7,13 = °7 ,14 = ° 


a A 

°8 1 = ( 1T ){(1-p 8 )[T 16 (H 8jl6 ) c U 4,l + T 17 (H 8jl7 ) c U 5,1 ]) S,L 
8 

OA 

°8,2 = (_ lT^ (1_P 8 )tT 16 (H 8jl6 ) c U 4,2 + T 17 (H 8jl7 ) c U 5,2^S,L 
* 8 


“ 8,3 - ° 


D 


8,4 


D 


8,5 


D 


8,6 


( !T ){(1 ' P 8 )lT l6 (H 8jl6 ) c U 4,4 + T 17 ( , H 8jl7 ) c U 5,4 ]} S,L 
8 

CfA 

+ T 16 (H 8J16>c U 4,5 + T 17°Wc V’Vl 

8 

OA 

( 1T ){<1 - p 8 )t6 6 (H 8 j6 ) c + t 16°Wc V« + ’lAlll’c U 5,6 1} S.L 


“8,7 ' ° 


C?A 

“8,8 = + T 16 (H 8jl6>c “4,8 + T 17 CH 8jl7 : 'c “s. 8 1 ' ‘s’s.L 


I 
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a A, 
*8 
aA 


a a 

D 8,9 - O ((1 -V [e 9 (H 8J9 > = + T 16 (H 8J16>c U 4,9 + •'l7 (H 8J17 ) c U 5,9 1J S,L 


8,10 “ (_ MT ){(1_P 8 )[T 16 (H 8jl6 ) c U 4,10 + T 17 (H 8jl7 ) c U 5,10 ]} S,L 
8 


Ok n 


D 8,ll ” ( 'ML ){(1 " P 8 )lT l6 (H 8jl6 ) c U 4,ll +T 17 (H 8jl7 ) c U 5,11 ]> S,L 
o 


Ok 


°8 ,12 " ( H fl ){(1 ' P 8 )[T 16 (H 8jl6 ) c U 4,12 + T 17 (H 8jl7 ) c U 5,12 1} S,L 


aA o 

D 8,13 " (- M7 ){(1_p 8 )[T 16 (H 8jl6 ) c U 4,13 + T 17 (H 8jl7 ) c U 5,13 3} S,L 

O 

ok 

D 8,14 = (_ M7 ){(1-p 8 )tT l6 (H 8jl6 ) c °4,14 + T 17 (H 8jl7 ) c U 5,14 ]} S,L 
8 


D 9,l ’ ^ ) <V 1 -‘>9 ) B 1 t 14 (I Wb U 2,l + t 16°Wb V 1 + Nc^c 

U 4,l + T 17 (H 9jl7 ) c U 5,1^S,L 

D 9,2 ■ ^'V^VVWb U 2,2 + 'l S <W B “b.B 1 + A 9c (1 - p 9>c 

' t 16 <h 9J16 > c U 4,2 + T 17 <H 9jl7 ) c U 5,2 1J S,L 


D 9,3 ’ ° 


°9,4 ' < 57> tA 9B (1 - p 9 B[T l 4 ( Vl4 ) B U 2,4 + t 16°Wb V + V (1 ' P 9>c 


U 4,4 + T 17 (H 9J17 ) C U 5,4 I5 S,L 
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D 9,5 " ^M 9 ^ A 9B^ 1 " C> 9^B^ T 14^ 11 9mI4^B U 2,5 + T 16^ H 9ml6^B U 3,5^ + A 9c ^ P 9^ 
l e 5 (H 9 J 5 > c + T 16 (H 9J16 >c U 4,5 + T 17°Wc U 5,5”s,L 

D 9,6 * <r )(A 9B (1 - p 9 ) B [T U (H 9ml4 ) B U 2,6 + T 16 (1 Wb U 3,6 ] + V^Vc 
9 9 

Ie 6 (H 9i6 ) c + T 16 (H 9.il6 ) c U 4,6 + T 17 (H 9jl7 ) c U 5,6 1} S,L 


°9,7 ' ° 


D 9,8 ' ( ^> {A 9B (1 - p 9 ) B [T 14 (H 9Ml4 ) B U 2,8 + t 16°Wb U 3,8 ! + Nc^^c 
+ T 16 <H 9jl6 ) c U 4,8 + T 17 (H 9jl7 ) c D 5,8 !) S,L 

D 9,9 ‘ <|-> (A 9B <1 - p 9>B I '9 (H 9„9 , B + T 14 °Wb U 4,9 + T 16 ( ^9jl6^c * A 9B ( 9 

+ VWc'Wc + T 16 (H 9jl6 ) c U 4,9 + T 17 (B 9jl7 ) c U 5,9 ] ‘ V *9^S,L 

D 9,10 ■ ( 5-> {A 9B (1 - p 9 ) B lT 14 <H 9»,14 ) B U 2,10 + T 16 °Wb Vo 1 + V T 16 
* 9 

(H 9jl6 ) c U 4,10 + T 17 (H 9jl7 ) c U 5,10 ]} S,L 
°9 11 ’ ( M" ){A 9B (1_p 9 ) B [e ll (H 9mll ) B = T 14 (H 9ml4 ) B U 2,ll + T l6 (H 9ml6 ) B U 3,ll ] 


+ A 9c (1 - p 9 ) c [T 16 (H 9 j l6 ) c U 4,ll + T 17 (H 9 j 17 ) c U 5,U 1} S.L 


°9,12 = ( M" ){A 9B (1_P 9 ) B [e i2 (H 9rnl2 ) B + T 14 (H 9ml4 ) B U 2,12 + T 16 (H 9ml6 ) B U 3,12 ] 
+ A 9 c (l-p 9 ) c [T 16 (H 9jl6 ) c U 4j12 + t 17 (H 9 . 17 ) c U 5>12 ]} S>L 
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D 9,13 = ( ^ ){A 9B (1 - p 9 ) B [T 14 (H 9ml4 ) B U 2,3 + VWm U 3,13 ] + 
lT 16 (H 9jl6 ) c U 4,13 + T 17 (H 9jl7 ) c U 5,13 ]) S,L 

°9,14 = ( ^" ){A 9B (1 ' P 9 ) B [t 14 (H 9mX4 ) B + T 14 (H 9ml4 ) B U 2,14 + T 16 (H 9ml6 ) B U 3,U ] 


+ A 9c a " p 9'c lT l6' n 9jl6 ; c °4,14 1 l 17 v “9jl7'c “5,14 JJ S,L 


(l‘Pn) U /. U + T 1 7^ H Q-n 7^r U S 14^ 


Ok. 


D 10,1 ■ ( Tl" )Ui ‘ p io ) Wi4>a “i.i’s.l 


10 


a A 


°10,2 ■ ( M^q ){( 1 " P 10 ) T l4 (H 10jl4 ) A U 1,2 } S,L 


D 10,3 " ° 


a A 


( M 10 ){(1 Pl0> (Hl °j U 1,4 } S,L 


10,4 


Ok. 


J 10,5 = ( M^ ){(1 “ P 10 ) T 14 (H 10jl4 ) A U 1,5 } S,L 


Ok. 


°10,6 * (- M^ ){(1 P 10 ) T 14 (H 10jl4 ) A U 1,6 } S,L 


°10,7 ° 


Ok 


Vs ■ ( T^ )!<1 -V WVa Ys’s.l 


ok 


°10,9 ■ ( P 10 ) T 14 (H 10jl4 ) A U 1,9 } S,L 
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10,10 


0 A 

■ ( lir ){( 1 - p io > V’Wo'a - e io + ( 1 ‘ p io ) Wiu’a “l.io's.L 


aA 1 

D 10,11 “ (_ H^ ){(1 " P 10 ) T 14 (H 10jl4 ) A U l,ll } 


10,12 


10,13 


10,14 


0 a 

< -5^ ){( 1 - p 10 )l£ 12 <H 10 j 12 ) A + V^OjuV U l,l 2 n s.L 

aA in 

T 14 (H 10jl4 ) A U 1,13*S,L 
OA 

^ - M^'^ < 14^ H 10jl4^A + T l4 ^ H 10jl4 A U l,14 


D l lf l “ ( -M^ ){(1_p ll )[T 14 (H llml4 ) B U 2,l + T 16 (H llml6 ) B U 3,1 ]} S,L 
oA 

O u , 2 - <-M ii)t<1 - p U )[T 14 <H Uml4 > B U 2,2 + T 16 (H llml6 ) B U 3,2 1} S,L 


D ll»3 ' ° 


aA, 


D ll,4 = ( — ){(1 ' P H )tT 14 (H H I nl4 ) B U 2,4 + T 16 ^ H llml6' ) B U 3,4 S,L 


OA - 

°11,5 = ^ m, 1 ^^ 1 ” P ll^ T 14^ H llml4^B U 2,5 + T 16^‘llml6 B U 3,5 S,L 


11 
aA. 
*11 
a A. 

*11 


OA . .. 

D ll,6 = ^ m,, ^ ^ 1-p il^ T 14^ H llml4^B U 2,6 + T 16^ H llml6 B U 3,6 S,L 


D U,7 = ° 


a A n 

D ll,8 “ ^^^^ (1 “ p ll )tT l4 (H llml4 ) B °2,8 + T l6 (H llml6 B U 3,8 S,L 


234 


a A. 


D 


11,9 


<_S 11 ){€ 9(Hllm9)B + Tl4<Hllinl4)B U2 ’ 9 + Tl6(Hllml6)B U 3,9 ]1 S,L 


aA 

D ll,10 = ^ M i;L ^ T 14^ H llml4^B U 2,10 + T l6^ H llml6^B U 3,10 S,L 
aA 

' ' 'U + T 14 (H nml4>B U 2,U 

+ T 16 <H llml6 ) B U 3,n's,L 


aA 


°11,12 = (- M^ ){(1_p ll )tt 12 (H llml2 ) B + T 14 (H llml4 ) B °2,12 


+ T l6 (H llml6 ) B U 3,12 1} S,L 


aA 


D ll,13 = (- M^ ){(1 ' p ll )[T 14 (H llml4 ) B U 2,13 + X l6 (H llml6 ) B U 3,13 ]} S,L 
aA 

D ll,14 = (_ M^ ){(1_P ll )tt l4 (H llinl4 ) B + T l4 (H llml4 ) B U 2,U 
+ T l6 (H llml6 ) B U 3,14^S,L 


°12,1 = ( ^ ){A 12A (1 - p 12 ) A T 14 (H 12 j 14 ) A U l,l + A 12B (1 - p 12 ) B [T 14 (H 12ml4 ) B U 2,l 
+ T 16 (H 12ml6 ) B U 3,1^S,L 

°12,2 = ( ^ ){A 12A (1 " P 12 ) A T 14 (H 12jl4 ) A U l,2 + A 12B (1 " P 12 ) B [T 14 (H 12ml4 ) B °2,2 
+ T l6 (H 12ml6 ) B U 3,2^S,L 


° 12 , 3 = ° 
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D 12,4 ' < M^ )tA 12A (1 - p 12 ) A T 14 (H 12J14 ! A “l,4 + W^IbVVW’b U 2,4 
+ T l6 (H 12ml6 ) B U 3,4 ]} S,L 

D 12,5 = ^ A 12A (1_P 12 ) A T l4 (H 12jl4 ) A U l,5 + A 12B (1_p 12 ) B [T l4 (H 12ml4 ) B U 2,5 
+ T l6^ H 12ml6^B U 3,5^S,L 


D 12,6 = ( ^ ){A 12A (1 ' P 12 ) A T 14 (H 12jl4 ) A U l,6 + A 12B (1 " P 12 ) B [T 14 (H 12»14 ) B U 2,6 
+ T l6 (H 12ml6 ) B U 3,6 ] ^S,L 

D 12,7 = ° 

°12,8 = ( ^ ){A 12A (1 ' p 12 ) A T 14 (H 12jl4 ) A U l,8 + A 12B (1 ‘ P 12 ) B [t 14 ' (H 12ml4 ) B U 2,8 
+ T l6 (H 12ml6 ) B U 3,8^S,L 

°12,9 “ ( ^ ){A l2A (1 ' Pl2) A Tl 4 (Hl2 j 1 4 )AUl , 9 + A 12B (1_P 12 ) B [( 9 (H 12m9 ) B 
+ T i4^ H l2ml4^B U 2,9 + T l6^ H 12ml6^B U 3,9^S,L 
°12 ,10 = ( ^ ){A 12A (1 - p 12 ) A [€ 10 (H 12 j 10 ) A + T 12 (H 12 j 14 ) A U l,10 ] + A 12B (1 - p i2 ) B 
t T l4 ^12ml4^B U 2,10 + T 16^ H 12ml6^B 3,10^ S,L 


D 12,ll " ( ^ ){A 12A (1 - p l 2) A T 14 (H !2 J1 4 ) A Yll + A 12B <1 “ P 1 2 ) B [4 1 1 (H I2^1 ) B 


+ T l4 (H 12ml4 ) B U 2,ll + T l6 (H l2ml6 ) B U 3,11 ]} S,L 
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D 12,12 = ( 4 ){A 12A (1 ^12 ) A t< 12 (H 12 j 12>A + T 14 (H l 2j 14>A U l,12 ] " A 12A <12 
+ A 12B (1 ~ P 12 ) B [€ 12 (H 12ml2 ) B + T 14 (H 12jl4 ) B U 2,12 + T 16 (H 12ml6 B 
U 3,12 ] " A 12B £ 12 } S,L 

°12 ,13 = ( M^' ){A 12A (1_P 12 ) A T l4 (H 12jl4 ) A U l,13 + A 12B (1 " P 12 ) B It 14 (H 12jd 14 ) B U 2,13 
+ T 16 (H 12ml6 ) B U 3,13^S,L 

D 12,14 = ( M^' ){A 12A (1_P l2 )[£ 12 (H 12jl4 ) A + T 14 (H 12jl4 ) A U l,14 + A 12B (1 " P 12 ) B 
t€ 14 (H 12ml4 ) B + T 14 (H 12ml4 ) B U 2,4 + T l6 (H 12ml6 ) B U 3,14 ]} S,L 

D 13 f l = ^M^'^ A l3F^ 1 ” P 13^F^ T 17^ H 13nil7^F °6,1 + T 13^ H 13ml3^F U 7,l^ A 13F °7,1 
+ A 13D ( 1 -P 1 3> D t f 1 ( H i3ji ) D + T 13 (H 13jl3 ) D U 8,l ] ' A 13D T 13 U 8,1 } S,L 

°13,2 = ^M^^ A 13F^ 1_P 13^F^ T 17^ H 13ml7^F °6,2 + T 13^ H 13ml3^F U 7,2 
+ *? (H 1 W } F ] " A 13F T 13 U 7,2 + A 13D (1_p 13 ) D [£ 2 (H 13j2 ) D 
+ T 13 (H 13jl3 ) D U 8,2 ] “ A 13D T 13 U 8,2 } S,L 

°13,3 = ° 

D 14,4 = < ‘M^^ A 13F^ 1_P 13^F^ T 17 < ' H 13ml7^F U 6,4 + T 13^ H 13ml3^F + f 4 < ‘ H 13m4^F^ A 13F 
T 13 U 7,4 + A 13D [ (1_P 13 ) D T 13 (H 13jl3 ) D ~ T 13 ] U 8,4^S,L 
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°13,5 = ( !^ ){A 13F <1 “ P 13 ) F lT 17 <H 13ml7 ) F U 6,5 + T 13 (H 13ml3 ) F U 7,5 ] " A 13F 
T 13 U 7,5 + A 13D [(1_P 13 ) D T 13 (H 13jl3 ) D “ X 13 ] U 8,5 J S,L 

°13,6 = ( S^ ){Al 3F (1 " Pl3)FtTl7(Hl3 m 17) F U 6,6 + T 13 (H 13ml 3) F U 7,6 ] " A 13F 
X 13 U 7,6 + A 13D [(1 " P 1 3) D T 13 (H 13jl 3) D “ T 13 ] U 8,6 J S,L 

D 13,7 = ° 

D 13,8 = ( ^ ){A 13F (1 - p 13 ) F [T 17 (H 13 m 17 ) F U 6,8 + U 7 ,8 ] ' A 13F 

T 13 U 7,8 + A 13D t(1-p l 3) D T 13 (H 13jl 3) D " X 13 ] U 8,8 J S,L 

°13,9 = ^M^^ A 13F^ 1_P 13^F^ T 17 ^ H 13ml7^F U 6,9 + T 13^ H 13ml3^F U 7,9-* A 13F 
T 13 U 7,9 + A 13 d [( 1 " P 13 ) D T 13 (H 13jl 3) D " X 13 ] U 8,9 } S,L 

°13,10 = ( ^ ){A 13F (1 " P 13 ) F [T 17 (H 13ml7 ) F U 6,10 + T 13 (H 13ml 3) F U 7,10 ] ' A 13F 
X 13 U 7,10 + A 13D [(1_P 1 3) D T 13 (H 13jl 3) D “ X 13 ] U 8,10 } S,L 

D 13,ll = ( M^ ){A 13F (1 ' p 13 ) F [X 17 (H 13rnl7 ) F Vll + X 13 (H 13ml 3) F U 7,ll 1 " A 13F 
X 13 U 7,ll + A 13D [(1 - p l 3) D X 13 (H 13jl 3) D “ X 13 ] U 8,H } S,L 

D 13,12 = ( i^ ){A 13F (1_p 13 ) F [x 17 (H 13ml7 ) F U 6,12 + X l3 (H 13ml 3) F U 7,12 ] " A 13F 
X 13 U 7,12 + A 13D [(1_P 1 3) D X 13 (H 13jl 3) D “ X 13 ] U 8,12 } S,L 
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°13 ,13 = ( M^'^ A 13F (1_p 13 ) F [T 17 (H 13ml7 ) F U 6,13 + X 13 ^Und^F U 7,13 + *13 
(H 13ml3 ) F 1 " A 13F (e 13 + T 13 U 7,13 ) + A 13D (1_p 13 ) D [e l3 (H 13jl3 ) D 
+ T 13 (H 13jl3 ) D U 8,13 ] " A 13D (€ 13 + T I3 U 8,13^S,L 

°13 ,14 = ( M^ ){A 13F (1 " P 13 ) F [T 17 (H 13ml7 ) F U 6,14 + T 13 (H 13ml3 ) F U 7,14‘ l " A 13F 
X 13 U 7,14 + A 13D I ' (1 " P 13 ) D 13 (H 13jl3 ) D ' X 13 ] U 8,14 } S,L 


°14,1 - ( m 7~ ){ A 14A X 14 [(1 - p 14 )(H 14 j 14 ) A ‘ 1] U l,l + A 14B [ ^ 1_P 14^ X 14 
14 

(H 14ml4*B ‘ X 14 1 U 2,l + A 14B X 16 (X_P 14 ) ^Und^B U 3,1^S,L 

°14,2 = ( M^ ){A 14A X 14 t(1_P 14 )(H 14jl4 ) A " ^ U l,2 + A 14B [(1_P 14 ) X 14 
(H 14«14 ) B ‘ X 14 ] U 2,2 + A 14B X 16 (1 " p 14 ) (H 14ml6 ) B U 3,2 } S,L 


14, 


14,5 


( M77 ){A 14A ^14t( 1 - p 14 )(H 14 J 14 ) A “ 11 U l,4 + A 14B X 14 [(1 - p 14 
14 

(H 14ml4 ) B " 1] U 2,4 + A 14B X 16 (1 " P 14 ) (H l4ml6 ) B U 3,4 } S,L 

^) {A 14A X 14 t(1 - p 14 )(H 14 j 14 ) A ’ 11 U l,5 + A 14B X 14 [(1 ' p 14 
(H 14ml4 ) B " 11 U 2,5 + A 14B X 16 (1-p 14 ) (H 14ml6 ) B U 3,5^S,L 


) 


) 
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D U 6 ■ ( ^ ){A 14A ' 11 U l,6 + A 14B T 14 1 

’ 14 

<Hi 4nl4 ) B - 1] U 2,6 + A 14B T 16 <1 ' P 14 )(H 1 4»16 )b Vb's.L 

D 14,7 ■ ° 

D 14 8 ■ ( fT7 ){A 14A T 14 [(1 - p 14 )(H 14hl4>A ‘ 1] U l,8 + Vb 'u 1 ^ 1 
5 14 

<H 14b 14>B - 11 U 2,8 + A 14B T 16 (1 - p 14>' 'W. V’s.L 

D 14 9 ■ < if _HA 14A • t 14 [(1 - p 14 >(H 14 3 14 ) A ‘ 11 °1,9 + 9 (I Wb 

’ 14 

+ T 14 (H 14ml4 ) B U 2,9 + T l6 (H 14ml6 ) B U 3,9 ] " Vb T 14 U 2,9 } S,L 

°14 ,10 = ^ {k UA^U )[( 10^Uil0\ + T 14 (K 14 j 14 ) A U l,10 ] ' A 14A T 14 U l,10 

+ A 14B (1_P 14 ) [T 14 (H 14ml4 ) B U 2,10 + T 16 (H 14ml6 ) B U 3,10-* " A 14B T 14 U 2,10*S,L 

°14 11 = ( r" ){A 14A T 14 [(1_p 14 )(H 14jl4 ) A " 1] U l,ll + A 14B (1 “ P 14 )[€ 11 (H 14«11 ) B 
* 14 

+ T 14 (H 14ml4 ) B U 2,ll + T 16 (H 14ml6 ) B U 3,ll ] " A 14B T 14 U 2,11 } S,L 

°14 12 = ( M 2-){A 14A (1 " P i4 )[e 12 (H l4jl2 ) A + T 14 (H 14jl4 ) A U l,12 ] " A 14A T 14 U l,2 
’ 14 

+ A 14B^ 1_P 14^ f 12^ H 14ml2^B + T 14^ H 14ml4^B U 2,12 + T l6^ H 14ml6^B U 3,12^ 

A 14B X 14 U 2,12^S,L 

D 14 13 = ( S^ ){A 14A T 14 [(1_p 14 )(H 14jl4 ) A “ 1] U l,13 + A 14B (1 “ p 14 ) It 14 (H I4«14 ) B 
5 14 

U 2,13 " T l6 (H l4ml6 ) B U 3,13 ] " A 14B T 14 U 2,13 } S,L 
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°14 , 14 = ( M^‘ ){A 14A [(1 ' P 14 )(H 14jl4 ) A " 11( ‘l4 + T 14 U l,l^ + A 14B (1 “ P 14* 
[t 14 (H 14ml4 ) B + T 14 (H 14ml4 ) B U 2,14 + T 16 (H 14ml6 ) B U 3,14 ] " A 14B 
(f 14 + T 14 U 2,14^S,L 


^ B 1^S,L T 13^ H ljl3^D U 8,15^S,L 

(B 2 ) S,L = ^ [A 2D (1_P 2 ) T 13 (H 2jl3 ) D U 8,15 + A 2F (1 " P 2 )tT 13 (H 2ml3 ) F U 7,15 
+ T 17 (H 2ml7 ) F U 6,15^S,L 

^Vs.L =( M^ )[A 4 (1 " P 4 ) ^ T 13 (H 4ml3 ) F °7,15 + T 17 (H 4ral7 ) F U 6,15^ S,L 
^ B 5^S,L = ( M7 )[A 5 (1 " P 5 ) ^16^16^ U 4,15 + T 17 (H 5jl7 ) c U 5,15* ] S,L 


( Vs,L = ^ [A 6 (1_P 6'’ {T 16 (H 6jl6 ) c U 4,15 + T 17 (H 6jl7 ) c U 5,15 }] S,L 

( Vs,L = ( M^ )tA 8 (1 ' P 8 ){T 16 (H 8.il6 ) c U 4,15 + X 17 (H 8jl7 ) c U 5,15^S,L 

( Vs,L = ( M^ tA 9B (1 " p 9 ){x l4 (H 9ml4 ) B U 2,15 + X 16 (H 9ml6 ) B U 3,15 } + A 9c (1_p 9 ) 

U 4,15 + X 17 (H 9jl7 ) c U 5,15 }] S,L 


10'S,L M lCr 1 14 v 10jl4 / A "l,15 1 15 v “l0jl5 / A J J S,L 


11' S,L ' ^/‘“ll^ M ll ,ll 14 v "ll«14'B "2,15 T l 16 v "llml6'B "3,15 JJ S,L 
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* B 12*S,L 


* b i 3 *s,l 


(B 14 5 S,L 


= ( M^' )tA 12A (1 ' P 12 ){T 14 (H 12jl4 ) A U l,15 + T 15 (H 12jl5 ) A } + A 12B (1_P 12 ) 
{t 14 (H l2ml4 ) B U 2,15 + T 16 (H 12ml6 ) B U 3,15 }] S,L 

= ( M^ )[A 13D T 13 {(1 - p 13 )(H 13 j 13 ) D " l}U 8,15 + A 13F T 13 ((1 - p 13 ) 
^ H 13ml3^F “ 1}U 7,15 + A 13F (1_P 13 ) T 17 (H 13ml7 ) F U 6,15 ] S,L 

- <r" )[A 14A T 14 {(1 - p 14 )(H 14jl4 ) A " 1>U 1,5 + T 15 

14 

fflujis 5 * + A 14B T 14 ta-»i4 , < H Wml4 ) B ‘ l)0 2,15 + A 14B T 16 
<1 " p 1 4 )<H 14ml6 ) B U 3,15 1 S,L 


For the total channel. 


( Vs,L 


T 13 * T 14 ' 1 > 


*13 “ ‘14 ‘ ° ' 


Scanner Model Equations 
Total Channel 


*i - Vj + Vj + V + B i E + b + e 


SW/LW Channels 


T. = A..T. +D..T 4 + B Q + (B ) E + (B ) E + b + e 
i ijj ljj Qi iSb ILL 


242 



Scanner (using Enclosure Coupling Method) 
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0 0 0 0 0 
0 0 0 0 0 
1 0 0 0 0 
0 0 0 0 0 


E 


2 


0 0 10 


0 0 0 0 0 0 

0 0 0 0 

E 1 + 

0 0 0 0 0 0 

0 0 0 0 


1 0 0 0 0 0 

0 0 0 0 


0 0 0 0 0 0 

0 0 0 0 


0 0 0 0 0 0 

_ _ 


L j 


- 

0 0 10 0 


0 0 0 0 

0 0 0 0 0 

E 2 + 

0 0 0 0 

0 0 0 0 0 


0 0 0 0 

0 0 0 0 0 


0 0 0 0 

0 0 0 0 0 


10 0 0 

0 0 0 0 0 


0 0 0 0 

_ _ 


- 




_ — 

0 0 0 0 1 0 


0 0 0 

0 0 0 0 0 0 

E 3 + 

0 0 0 

0 0 0 0 0 0 


0 0 0 

0 0 0 0 0 0 


10 0 


0 0 0 1 
0 0 0 0 
0 0 0 0 


= 0 


Thus D matrix is formed by the above Es. 
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APPENDIX E 


ERBE SOLAR MONITOR THERMAL MODEL 


FACE BLANK NOT FILMED 




SOLAR MONITOR MODEL 


As described in the non-scanner and scanner modeling efforts, the same 
notations are used for constants and variables in developing the energy 
balance equation of the solar monitor. 

To avoid the burden to write the same equations as used in developing 
the other models, the same variable name in the equation such as is 

here simply used. 

A general form of the equation for the solar monitor which is divided 
into 19 nodes is presented here as 

• 4 

T = AT + DT +BqQ + BE + b. 

Now the computational algorithm for the constants of the above equation 
has been developed by using the enclosure coupling technique. It is avail- 
able and can be used for all 8 non— scanner, scanner and solar monitor models. 

The following constants, A^ , represent the coefficients of themal 
conduction terms in the energy balance equations. 


An - - l/CR^V 


A i2 * ° 


A 13 1/(R 13 M 1> 


A 21 ' ° 


A 22 11 <R 23*y 


A 23 ' 1/<R 23 M 2> 
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c 


1 2 3 4 5 6 



C = (I - p d E) 1 E T • 

* ■ [ (I - P - TD) (I - CD) 1 A - + [ (1 - p- tD) (I - CD) 1_ C - x] F S 


253 





A 31 = 1 ^ R 31 M 3^ 


A 32 " 1/(R 32 M 3 ) 


A 33 = -<T- + R ~ + -T~ )/ll 3 
33 R 31 R 32 R 34 3 


A 34 = 1 ^ R 34 M 3^ 


A 43 = 1/1 (R 43 M 4 ) 


a 44 ’ ' ( f“ + + r ~ + r^— )/ m 4 

44 R 43 R 45 4 7 R 4-19 


\5 - ^^ R 45^4^ 


A 47 = l/CR^^M^) 

A 4-19 " 1/(R 4-19 M 4 ) 


A 54 - i/(r 34 m^) 


A 55 ■ - ( r~ + F^ )/M 5 
55 E 54 R 56 5 


A 56 ' 1/( W 


A 65 ~ ly ^ R 65*V 


A 66 - -<4 + 4 >/Ma 


A 67 " 1/(R 67 M 6 ) 
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SOLAR MONITOR MODEL 

1. Active Cavity 

2. Reference Cavity 

3. Button, Duplex Cone 
Cavities 

4. Heat Sink 

5. Aperture, Primary 

6. Shield, Aperture 

7. Body (Rear Barrel) 

8. Body (Front Barrel) 

9. Aperture, Secondary 
(or FOV Limiter) 

10. Motor Housing 

11 . Chopper Blade 

12. Chopper Pendulum 

13. Motor, Chopper 

14. FOV Limiter 

15. Front Top Housing Cover 

16 . Rear Top Housing Cover 

17. Front Housing 

18 . Rear Housing 

19. Holder, Heat Sink 

20. Fictitious Surface 
p. Temperature Probe 

E. Total Incident Radiation 
E. • Irradiance in enclosure a 
E. Irradiance in enclosure b 
E. Rate of Incident Radiation 
From a Neighboring Surface 
a-h. Cavities Within The 

Sensor, i.e.. Enclosures 
or Spaces Between Sensor 
Nodes 
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A 74 1/(R 74 M 7 ) 


A 76 ' 1/<R 76 M 7> 

A 77 ■ " <r ^ + ^7 +r ^7 )/M 7 
A 78 = 1/(R 78 M 7 ) 


A 87 * 1 ^ R 87 M 8^ 

A «l(-i- + -i- + — - — ^ M 8 
88 R 87 R 89 R 8-10 

A 89 “ 1/ ^ (R 87 M 8 ) 


A 8-10 “ 1 /( R 8 - 10 M 8 ) 


A 9_8 1 ^ R 9_8 M 9^ 


9-9 R 9-8 M 9 


A 10-8 " 1 /( R 10 - 8 M 10 ) 


A 10-10 1 /( R 10 - 8 M 10 ) 


1/(R 10-12 M 10 ) " 1/(R 10-13 M 10 ) 


A 10-12 1 /( R 10 - 12 M 10 ) 


A 10-13 1 1/(R 10-13 M 10 ) 


A ll-12 1/(R 11-12 M 11 ) 


A ll-ll 1/(R 11-12 M 11 ) 
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A 12-ll 

^ 2-12 

A 12-13 

A 13-12 

A 13-13 

^ 3-10 

A 14-14 

A 14-15 

A 14-17 

A 15-15 

^ 5-16 

A 15-14 

A 16-15 

^ 6-16 

A 16-19 


1/(R 12-llV 

■ <1/R 12-11 + 1/R 12-13 )/M 12 
1/(R 12-11 M 12 ) 


1/(R 13-12 M 12> 


■ (1/R 13-12 + 1/R 13-10 )/M 12 


1/(R 13-10 M 13 ) 


* (1/R 14-15 + 1/R 14-17 )/M 14 


1/(R 14-15 M 14 ) 

1/(R 14-17 M 14 ) 


- 1 /( r 15 _ 16 m 15 ) - 1/(R 15-14 R 15 ) 
1/(R 15-16 M 15 ) 

1/(R 15-14 M 15 ) 


1/(R 16-15 M 16 ) 

" (1/R 16-15 + 1/R 16-19 )/M 16 


1/(R 16-19 M 16 ) 
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A 17-14 1/(R 17-14 M 17 ) 


A 17-17 " - (1/R 17-14 + 1/R 17-18 )/M 17 


A 17-18 1/(R 17-18 M 17 ) 


A 18-17 1/(R 18-17 M 18 ) 


A 18-18 (1/R 18-17 + 1/R 18-19 )/M 18 


^8-19 1/(R 18-19 M 18 ) 


A 19-4 = 1/(R 19-4 M 19 ) 


A 19-16 ~ 1 ^19-16 M 19^ 


A 19-18 " 1/(R 19-18 M 19 ) 

A 19-19 = " (1/R 19-4 + 1/R 19-16 + 1/R 19-18 ) /M 19 


For the radiative interaction term, considering the specular reflection, 

( 1 ) 
( 2 ) 


d s 

P = P + P 


€ + p + T = 1 


< ■ ‘A + p k E k + \\ 


(3) 


E - l E. . M° 
k h k-j j 
J 


(4) 


The energy balance equation is now 


♦k * \ ” M k 


(5) 
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where 


»k 


k\ + < 


Thus 


K, = 


\ - ‘A 


d. s — 

p. E, - p, E, - T. E. 
Ik M kk k k 


where 


n k ■ ° T t • 


From (3) and (4), 
, n 


- c jVi rl( 'i n i + T iV 


From (7) through (4) 


E k - | V, l + T i E 1 ) 


n 

V 


- 7(6 .ft. + t.E.) ) E 
,11 1 1 . 
i J =1 


- - - / . r . . 

li i i . L , k-j j-i 


Equation (6) becomes 


* k - (1 - °k )E k - ‘A - W 


■ p (1 - V'i | E k-j p 3 -i - ‘At 1 fi i 

- ? [<1 - V T 1 1 - T k s ki) E i 

1 J J 


Let 


H. 


kji 


■l 


E. .P. . 
J k-J 3-1 


H 


, V 


E P 


nm£ ^ n-m m-Z 


m 


(6) 


(7) 


( 8 ) 


(9) 
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From the above drawing, the enclosures, a, b and c, are found. For 
each enclosure, the incoming radiance to a certain surface can be described 
by Equation (8). 

Choose the enclosure b and c and couple them together, then one can 
define the transmitted radiance through the filtering surface. 

First , 

enclosure b - elements 6,7,8,9,11,20 

enclosure c - elements 1,5,20 


Index 

i,j,k 

^ , m , n 


From Equation (8) 


e, = y e .n. y e. .p. . + y t.e. y e .p. . 

k L . l l h k-j j-i . ii“ k-j j-i 


where in this equation i, j and k are dummy index 


or 


E, = y e .Q. H . . + y t .E. H 

k h i i kii r ii kj i 

l i 


For the enclosure c. Equation (8) becomes 


E ?0 ~ e £ fi £ H 20m£ ) c + X 20 E 20 (H 20m20 ) c 


(10) 


where £, m and n = 1, 5 and 20. 
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Likewise, for the enclosure b 


E 20 ” 4 i^i H 20ji^b + T 20 E 20 ^20^20^ + T 11 E 11 * H 20.111* 


20 j 20 


20jll b 


( 11 ) 


where i and j =6, 7, 8, 9, 11 and 20. 


Combining (10) and (11) to obtain E 2 Q and E 2 Q as a function of E-^> we 


have 


E 20 < A H 20m£ ) c + T 20 (H 20m20 ) c € i fi i H 20ji*b 


+ T 11 T 20 (H 20m20 ) c * H 20jll*b E ll ] 


( 12 ) 


and 


E 20 G l t( ^ f i^i H 20j i^n + T 20 (H 20j20 ) b *A H 20irt'\ 


+ T ll (H 20jll ) b E ll ] 


(13) 


where 


G 1 " T 20 (H 20j20 ) b (H 20m20 ) c 1 


The radiation energy balances in the enclosures b and c are 


*k|. ■ ? [<1 - V 6 1 V - 'Ai 1 n i 

1 b 1 J 


+ ? l<1 - Pfc )T i H kji • 'Al 1 <E i=20 ° r 


(14) 
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and 


n 


- ltd 

i 


' p n> e Am l ' % 6 nl ] *1 


+ lid 

l 


" P n ) T £ H nm£ T n <S n l ] h 


(15) 


elements 11,14,21 
elements 6,7,8,9,11,20 


Index 

£,m,n 

i,j 


Second , 

enclosure d 
enclosure b 

From Equation (8) 


E 11 " € i^i H llj i^b + T ll (H lljll ) b E ll + T 20 (H llj20 ) b E 20 


(16) 


and 


J 11 


(| + T il (H ll m ll ) a E 11 + T 21 (H llm21 ) a E 


Substituting (12) into (16) , Equation (16) becomes 


(17) 


E 11 " € i R i H llji ) b + T 20 G l (H llj20 ) b * A H 20m£*c 
+ T 20 G l (H llj20 ) b (H 20m20 ) c € i fi i H 20ji*b 

+[T ll (H lljll ) b + T 11 T 20 G l (H llj20 ) b (H 20m20 ) c (H 2-jllV E 11 (18) 

Substituting (18) into (17) 

E 11 = G 2 [( | € A H nm£ ) a + T 21 (H llm21 ) a E + T ll (H llmll ) a ‘i^Alji^b 
+ T 11 X 20 G l^ H llmll‘ ) a ( ' H llj2o\ € l Q l H 20m£^c 


+ T 11 T 20 G l* H llmll*a (H llj20 ) b < ‘ H 20m20 ) c 6 i^i H 20j iV (19) 
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where 


G 2 [1 T 11 S l (H llmllV 


S 1 " AljlA + T 20 G l (H llj20^b (H 20m20 ) c (H 20jll ) b 


On the other hand, substitution (17) into (18) 

E 11 = G 2^^ € i n i H llji^b + T 20 G l (H llj20 ) b 6 A H 20mA 
+ T 20 G l^ H llj20 ) b (H 20m20^c ^ € A H 20ji*b 
"*11 S l ( | * A H llm^^a + T 11 T 21 ^AlstfA E] 

For the enclosure a: nodes: 11,14,21 


( 20 ) 


—0 for open position 


K 11 


[(! - p 1]L ) e 11 (* J i lmll ) a " ^lAl + (1 “ P 14 ) e 14 (H llml4 ) 


CL, 

a 14 


+ [1 - p u ) T 2i (H n m 2l ) a E + t(1 p ll ) ^l^llmlA " T 11 ]E 11 

0 for closed position 


p 14 


f. 


.0 for open position 


-(! - P i4 > /n^UmlA n n + [(1 P 14 J G 14 (H 14ml4 ) a ‘ *14 1 fi 14 


+ (1 p 14 ) T 11 ( K 14inll ) a E n + P 14 ) T 21 (H 14m21 ) a E 
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For the enclosure b: nodes: 6,7,8,9,11,20 


* K 1 - P 6 ) € 6^ H 6j6^b € 6 ] P 6 + (1 p 6 ) *7 (H 6j7 ) b P 7 


.0 for open position 


+ (1 " p 6 ) € 8^ H 6j8^b P 8 + (1 p 6 ) f 9 (H 6j9 5 b P 9 + (1 P 6 ) ^11 ^jll^b °11 

0 for closed position 

+ (1 - P 6 ) ^ii^gjU^b \l + (1 P 6 ) T 20 (H 6j20 ) b E 20 


= (i - p y ) € 6 ( H 7:j6 ) b Q 6 + t(l " P y ) € 7 ^ H 7j7^b " *7 ] Q 7 


+ (1 - p ? ) £ 8 ( H 7j8 ) b ^ 8 + (L ~ P 7 ) 6 9^ H 7j9^ b S + (1 ‘ P 7 ) ^ll^jll^b R 11 


)/ 


0 for open position 


+ (i - P 7 ) /ll ^ H 7 j 11 ^ b E 11 + ^ P 7^ T 20^ H 7j20^b E 20 


0 for closed position 


" (1 - p 8 } * 6 (H 8j6^b n 6 + (1 ‘ P 8 } < 7 CH 8j7 ) b R 7 + [ (1 “ P 8> 6 1 


^ H 8j8^b “ C 8^ P 8 


+ (1 pg) ‘9^ H 8j9^b P 9 


0 for open position 

+ (i - p 8 ) - Q n 


+ (1 p 8 ) /u^gjll^b E n + (1 " <V T 20 (H 7j20 ) b E 20 




0 for closed position 
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(1 p g ) £ 6 ( H 9 j 6 ) b + (1 - Pg) € 7 ^ H 9j7^b ^7 + ( L " Pg^ £ 8^ H 9j8^b fi 8 


0 for open position 

+ [(1 - P 9 ) { 9 ^ H 9 j9^b " { 9 ] S + (1 " P 9^ ^ll^jll^b Q 11 


+ P^y^Ll/Vjl^b E 11 + (1 p 9 ) T 20 (H 9j20 ) b E 20 


0 for closed position 


K llL ” (1 P ll^ € 6^ H llj6^b n 6 + (1 ” P ll ) € 7 ^lljy^b U 7 + (1 P ll ) *8 
^ H llj8^b ^8 


+a - p u > f 9 ( H n j 9 ) b + ( (1 p n ) ^{i (H iijii ) b " f n ] n n 

0 for open position 

+ ( L - P u ) T 20 (H llj20 ) b E 20 + f(1 ■ p ll ) yil (H lljll ) b ” T ll ] E 11 

0 for closed position 


^20 j b (1 “ P 20 ) € 6 (H 20j6^b ^6 + (1 P 20 ) e 7 (H 20j7 ) b ft 7 


+ (1- P 20 ) e 8^ H 20j8^b % + ^ P 20^ f 9 ^ H 20j9^b fi 9 + ^ P 20^ "11 


Ofor open position 



(H 20jll 5 b fl ll 


+ (1 P 20 ) T n (H 20jll ) b E 11 + [(1 P 20 ) T 20^ H 20j20^b ” X 20^ E 20 
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Enclosure c: Nodes: 1,5,20 


♦J - Id - P X > *!%!>. - ‘l 1 "l + (L - p l> E 5 

|c 

(1 - P-^ T 20^ H lj20^c E 20 

<J >5 = (1 - P 5 ) n ]_ + t(1 " p 5 ) € 5 (H 5h5*c " *5 J 

c 

A 

+ (1 - p 5 ) T 2 o (H 5j20 ) c E 20 

^20 = (1 “ p 20 ) € l (H 20jl ) c n i + (1 " P 20 ) f 5^ H 20h5^c S 

c 

/s 

+ td - p 2Q ) T 2 0 (H 20j20 ) ~ X 20 ] E 20 

For the enclosure d: Nodes: i,j,k = 1,3, 4, 5 

♦i = Id-pj.) e i (H iji ) d _£ i ] fi l + (1 ‘ p l ) fc 3* H lj3*d S 
d 

+ (1 - p x ) e 4 ^4 + (1 P l^ € 5 (H ij5^d S 

(J> 3 = (1 - p 3 ) n i + [(1 _ p 3* € 3 (H 3j3 ) d " C 3^ ^3 

d 

+ Cl - P 3 ) t 4 < H 3 :) 4> d n 4 + (1 ■ P 3 > ‘ 3 <li 3j5 > d ”5 

*4 ‘ (1 - p 4> W'd B 1 + U ' P 4 > Wd S 

d 

+ Id " P 4 ) £ 4^ H 4j4^d ~ f 4 ] ^4 + (1 " P 4 ^ *5^ H 4j5*d S 
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- (1 - p 5 ) ‘l <H 5Jl > d °1 + (1 ' 0 5 > f 3 "Vd fl 3 


+ <i-p 5 > 6 4 (H 5J4 ) d SJ 4 +U1-P5) < 5 <H 3j5 ) d - .3] n 5 


For the enclosure e: Nodes: 4,16,18,19 


[(! p 4 ) < 4 < H 4 j4> e f 4 ] n 4 + d P 4 ) * 16 r ‘l6 


+ (1-P 4 ) * 18 fl 18 + ^ 1-p 4^ *19 * H 4jl9*e Q 19 


'l6 


(1 P 16 ) *4 (H 16j4 ) e ^4 + t(1 " P l6 ) *16 (H 16jl6 ) e " *16 ] fl 16 


+ (1“P 16 ) £ lg ^gjis^e n i8 + (1_p 16 ) * 19 n i9 


18 


(1 " p 18 ) *4 (H 18 j 4 ) e fi 4 + *16 Vj^e °16 


+ t < 1_ P 18 > e 18 ^ H i8jl8^e _ *18 J fi 18 + (1 p 18 } *19 ^ H 18jl9*e S 19 


*19 1 e (1 P 19 ) *4 (H 19j4 ) e fi 4 + (1 * p 19 ) *16 (H 19jl6 ) e fl 16 

+ (1 “P 19 ) c i8 (n^j^e R 18 + ^ 1-p 19* *19 ” *19"* “19 


For the enclosure f: Nodes: 4,7,8,14,15,16,17,18 


= t(l-P 4 ) f 4 (H 4j4 ) f - e 4 ] ^ + (l-p 4 ) c ? (H 4j? ) f 


(! P 4 ) « 8 ( R 4j 8 )f fi 8 + * 1_p 4* *14 (H 4jl4 ) f fi 14 + (1 ' p 4 ) *15 * H 4jl5*f fi 15 


+ (1 P 4 ) * 16 ( H 4jie) f ^X6 + (1 p 4 ) *17 (H 4j 17^f R 17 + (1 p 4^ *18 * H 4jl8*f fl 18 
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(1 P 7 ) l 4 (Hy f + [(l“Py) « 7 (^7 j 7 > f € 7 ^ ^7 + ^ p 7 ^ £ g 


(H 7j8 ) j fi 8 


+ (l-p 7 ) « 14 ( H 7 j 14 ) f n !4 + ( 1_ P 7 ) e 15 ^jis^f n i5 + (1 ”<V e l6 


n i6 


+ (l-p ? ) « 17 (H 7 ji7 ) f n i7 + ( 1_ P 7 ) e 18 ^18 


(l- p 8 ) ^ H 8j4^f ^4 + ^ 1_P 8^ i 7 ^ H 8 j 7 ^ f ^7 + ^ 1 ~ p 8^ *8 


^ H 8j8^f “ f 8-* ^8 


+ (1 p 8 ) * 14 (H 8jl4 ) f n i4 + (1 ' p 8* * 15 (H 8jl5*f n i5 + ^ 1-p 8^ *16 
(H 8jl6 5 f n i6 

+ (l-Pg) f 17 ^ H 8ji7^f Q 17 + ^ 1-p 8^ e l8 ^ H 8jl8^f Q 18 


"141 


(1-p l4 ) e 4 (H 14jl4 ) f ^4 + (1 p 14 ) f 7 (H 14j7 J f S + (1_P 14 ) f 8 


* H 14j8*f S 


+ [(1"P 14 ) « 14 ( H !4ji4> f ( L ^ Q U + (1 p i4 ) ( 15 (H 14jl5*f ft 15 + 

(1 * p 14 ) f 16 (H 14jl6*f fi 16 


+ ( 1 “P 14 ) f 17 (H i4ji7 ) f ^17 + (1 ' P 14 ) e 18 (H 14J18*f fl 18 
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4>l 5 l f = ( 1 _ P 15 > ( 4 (H i5j4 ) f ^4 + ^" p 15^ i 7 (H 15j7 ) f ft 7 + (1_p 15 ) *8 

+ (1"P 15 ) e 14 fi i4 + t (1- P 15 ) e !5 ^iSjiS^f " € 15 3 fi 15 + 

( 1 -p l5> £ 16 ^16 

+ ( 1 *P 15 ) 4 17 (H i5jl7*f fi 17 + (1 " p 15 ) 6 18 (H 15jl8 ) f ^18 


*16 L ^” p 16^ *4 * H 16j4*f ^4 + (1-f W £ 7 (H 16j7 ) f fi 7 + (1 " P 16^ *8 


^ H 16j8^f ^8 


+ (1-P 16 ) e l4 (ni6jl4^f fi 14 + (1_p 16 ) € 15 (H 16jl5^f n i5 + ^ 1_p 16* *16 
* H 16jl6*f “ € 16 ] n i6 

+ ( 1 _ P 16 ) f 17 (H i6jl7^f fi 17 + (1-p 17 ) € 18 ^lejlS^f P “18 


k 17 


(1_P 17 ) ( 4 (H 17j4 ) f ^4 + (1 p 17 ) € 7 (H 17j7 ) f fi 7 + (1 * p 17 ) *8 


(H 17j8 ) f ^8 

+ (1”P 17 ) ( 14 (H i7ji4 ) f n i4 + ( 1 ~Pi7 ) 4 15 (H l7jl5 ) f fi 15 + (1_p 17 ) * 16 
(H 17jl6 ) f ®16 


+ [(1 p 1? ) e 17 (H 17 ji 7 ) f " € 17 1 fi 17 + (1 P 17 ) e 18 (H i7jl8 ) f P ‘18 
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For the enclosure g: Nodes: 8,9,14 

^8 = [(1 ' p 8 ) £ 8 ^ H 8j8^g ' f 8 ] U 8 + (1 " P 8 ) *9 ^9^ S 

g 

+ (1-Pg) f 14 ^ H 8jl4^g n i4 

4-9 - d-P 9 ) € a (H 9j8 ) g ^ + t(l-p 9 ) « 9 (H 9j9 ) g -* 9 ] « 9 

g 

+ (1-P 9 ) < 14 Q 14 

^14 = (1-P 14 ) € 8 (Il 14j8 ) g fi 8 + (1_p 14 ) € 9 ^l^^g R 9 

g 

+ [(l-p l4 ) e l4 (H i4ji4^ g " f 14^ R i4 

For the enclosure h: Nodes: 8,9,14,17 - i,j,k 

^8 ^ 1_p 8* f 8 ^ H 8j8^h ” f 8-* 0 8 + (1_p 8 ) *9 (H 8j9\ ^9 

h 

+ ( 1 ~P 8 ) £ 14 ^ H 8jl4^h a i4 + (1_P 8 ) *17 R 17 
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^ql “ (1~Pq) € q (Hoio)i, [(1 P Q ) { Q (^q^q)u ~ f cJ ^ 


|h 


V 8 9j8h 8 


9 ' 9 9j9 h 9 9 


+ (1 P 9 J *14 * H 9jl4*h ^14 +(1 P 9 J * 17 ^ H 9jl7^h fl 17 


^14 j h ~ (1 p i4 ) * 8 * H 14j8*h fl 8 + (1-p 14 ) 6 9 * H 14j9*h °9 

+ t(l-p l4 ) « 14 (H 1 4ji4 ) h ' € x4l fi i4 + (1_ Pi4 > c i7 * H l4J17*h fl 17 


17 


h " (1_P 17 ) *8 ^17:8^ % + (1 ’ p 17* € 9 (H 17j9 ) h fl 9 


+ (1“P 17 ) « 14 (H i7ji4^h ^14 + [(1_P 17 ) € 17 (H 17jl7 ) h ~* 17 ^ fi 17 


For the enclosure I: Nodes: 2,3,4 - i,j,k 


[(1-P 2 ) i 2 (H 2 j 2 ) i - e 2 ] fi 2 + (1 p 2 ) e 3 (H 2 j 3 ) i 
+ (1~P 2 ) « 4 ^ H 2j4^i ^4 


d- p 3 ) e 2 ( H 3 j 2 ^i S + ^ 1 -p 3 ^ e 3 ^ H 3j3^ f 3^ S 

+ < 1 - p 3 ) 6 4 (H 3j 4>i ^4 


= (1 ' p 4 } *2 (H 4 j 2 ) i ^2 + (1 - p 4> * 3 (H 4 j 3 ) i °3 


+ [(1-p^) f 4 ^4j4^i € 4^ ^4 
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The following are the constants for the radiation transfer terms in 


the energy balance equation. 


G 1 = ^ T 20 ^ H 20j20^b ^ H 20m20^c 3 

S 1 = ^lljll^b + T 20 G 1 <H llj20 ) b (H 20m20 ) c -^Ojll^b 
G 2 ^ “ T 11 S 1 ^llmll^a 3 

U k " G 2 £ k [ (H llhk^b + T 20 G 1 (H llj20 ) b (H 20m20 ) c (H 20jkV 

V k = G 1 € k T 20 [(H 20j20 } b (H 20mk ) c + T 11 G 1 °2 ^Ojll^b ^llmll^ 


^ H 20mk ) c 3 


2 2 


Y k = G 1 € k ^Onrf^c + T 11 T 20 G 1 G 2 (H 20m20 ) c ^Ojll^b (K llmll\ 


( Wc ] 


W k _ G 1 I(H 20jk J b £ k + T 11 (H Hmll*a (H 20jll } b V 


D 11 = ( r> {A lc [(1_P l ) £ 1 (H ljl ) c - £ 1 ] + A lc ^ T 20 (H 1 J 20 ) c V i 


+ A^Kl-p.^ < 1 ■ e i 33 


°12 = ° 


°13 " A id (1 - p l } £ 3 < H l J 3 ) d 


°14 - A ld (1 - p l> € 4 (H l j 4 ) d 


®llj20*b 


^ H llj2lPb 
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°15 ^ ^ A lc ^ ^ 1_p i^ ‘5 * H lj5*c + (1 P l ) T 20 (H lj20 ) c V + A ld (1 ' P 1 ) 

€ 5 (H lj5V 


°16 ~ A lc (1 - p l> T 20 (H X j2 0>c W 6 


°17 ( M 1 ) A 1c (1 P l ) T 20 (H lj20 ) c W 7 


°18 ~ ( M 1 ) A lc * 1-p l* T 20 (H lj20 ) c W 8 


°19 A lc (1_P 1^ T 20 (H lj20 ) c W 9 


° 1-10 = ° 


D l-ll ( M 1 ) A lc (1 p l* T 20 (H lj20 ) c [W 11 + °1 °2 T 11 € 11 ^Ojll^b (H llmllV 


°1-12 D l-13 ° 


°1-14 ( M 1 ) A 1c ^ 1-p l* T 20 (H lj20 ) c T 11 *14 G 1 °2 ^Ojll^b ^llndVc 

D l-15 = D l-16 = °1-17 = D l-18 = °1-19 = 0 


B 1 ( M 1 ) T 11 X 20 T 21 (1_p l ) G 1 G 2 (H lj20 ) c ^Ojll^b (H llm21^a 


°2-l = ° 


°2-2 ( M 2 ) A 2i [(1 ~ P 2 ) € 2 (H 2j2 ) i *2 ] 
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D 2-3 " A 2i (1_p 2 ) *3 <H 2j3*i 


D 2-4 ' <SJ> A 2 i (1 - p 2 > '4 <H 2j4 ) i 


J 2-S " D 2-6 " D 2-7 " D 2-8 “ D 2-9 °2-10 ~ D 2-ll D 2-12 2-13 2-14 

“ D 2-15 " D 2-16 ” D 2-17 ’ °2-18 " D 2-19 “ 0 


Yl - A 3d 0-»J> 6 1 <H 3Jl>d 

D 3-2 " <57> A 3d ( W € 2 «3 J 2 ) 1 


D 3 . 3 - ( 2 -) {A 3d [(l-p 3 ) e 3 (H 3j3 ) d - < 3 ] +A 31 ((1-P 3 ) r 3 (H^), 


* 3 » 


3-4 


( §^ f A 3d (1-P 3 > < 4 (H 3j4 ) d + A 3i (1-P 3 ) < 4 


D 3-5 - A 3d a - p 3> ‘5 V* 


D 3-6 ■ °3-7 ■ D 3-8 ’ °3-9 * u 3-10 ' “3-11 " “3-12 * “3-13 ‘ “3-14 ~ ”3-15 
" D 3-16 " °3-17 " D 3-18 " ®3-19 " 0 


■ 0 . 


- D. 
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